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> CHAPTER I 

DIGITAL COMPUTERS AND THE TEACHING 
OF COLLEGE GEOGRAPHY 

Although the large digital computer represents a relatively recent innovation in 
American universities, it is already established as a major tool for research 
geographers There is, in geographic research as in other disciplines, a definite 
interaction between tools and problems. Some problem areas may he dormant for a 
period of years or even decades because the necessary research tools are unavailable*. 
When a new research tool becomes available, either through primary innovation or 
through the transfer of knowledge from workers in other disciplines, there is 
frequently a surge of interest in topics which had previously been set aside, as well 
as development of new areas of investigation. 

An example of this phenomenon wrthin geography can be seen in the impact of the 
high-speed, digital computer on problems of multivariate regionalization/ The 
approach to an optimal multivariate classification of areas based upon a comparison 
of between region and within region variances had been developed, but empirical 
implementation was forced to wait until computing machinery capable of handling 
the calculations became available. When it did, a strong upsurge of empirical 
interest in this topic followed. These empirical investigations raised, in turn, 
additional theoretical questions and the cycle advanced through another phase. 

Parallel to the development of the computer as a research tool in many 
disciplines has been its use as an instructional aid. Early applications of the digital 
computer in instruction involved the use of the machine to permit students to 
handle problems of a size and scale of complexity that had hitherto been impossible 
within the standard classroom situation. In engineering, for example, students 
began to handle problems which* were considerably more realistic than the 
traditional, hand-computed examples. This proved quite fruitful, and the intro- 
duction o' the computer as a laboratory tool proceeded rapidly in many 
engineering disciplines and, in some schools, its use has been firmly established for 
well over a decade. An illuminating discussion of this approach can be found in the 
report of the Ford Foundation project at the College of Engineering of the 
University 'of Michigan [17). Instructional use of the computer has lagged 
considerably in most of the social sciences, though there have been some recent 
attempts witfiln geography and other social science disciplines to provide programs 
of this nature. For example, see the work by Rushton [II] as well as the more 
elaborate urban simulations such as that of Warncryd [ 19) . 
t Recent technological developments have permitted the central computer to 
handle a large number of users on a nearly simultaneous, real-time basis. Efficient 




time-sharing systems permit a number of users, equipped with remote terminals, to 
utilize a central computer through the intermediate efforts of a complex supervisory 
system which examines the priorities and needs of each user and allocates to him 
only that portion of the main computer's capacity required for his work, and for 
only that period of time which is necessary to process the next step in his problem. 
Much of the capaj^ty of these time-sharing systems arises out of the great disparity 
between the internal computation spee'ds of the central computer and the speed at 
which information can be transmitted and received from the remote terminals. This 
disparity can be seen in an examination of the relative printing rates of most of the 
commercially available remote terminals. These are very slow when compared with 
the main computer's dedicated peripheral printers and the speed factor is 
frequently on the order of 250 or 300:1, if not , more. Even with this speed 
differential, the main computer is frequently bound by the capacity of its "high 
speed" printers and will operate two or more of them simultaneously in order to 
attain greater output volumes. 

The introduction of time-sharing facilities m many universities has prompted the 
development of computer assisted instruction (CAI). There has been a great deal of 
interest in CAI on all educational levels ranging from preschool children to college 
graduates, and many of the CAI developments in college geography are discussed in 
previous technical papers of the Commission on College Geography. In the standard 
CAI environment, the time-sharing computer interacts with one or more students 
utilizing remote terminals. A series of queries are presented by the computer, based 
upon a program previously prepared by the instructor, with each query in the 
program designed to illuminate a particular point relating to the topic under 
discussion (e.g.. map projections, water balance, etc.). After each query the student 
presents a response which the computer checks against a series of stored answers 
which have also been provided by the instructor, and then, depending upon the 
nature of the student's response, branches to another section of the program. 
Internal loops are common in these programs so that the student can be given an 
opportunity to re-examine material that is giving him difficulty. This "match, 
branch and loop" format is basic to most CAI programs. In many cases, the CAI 
programs are not really computational ones since they utilize only a very small 
portion of the available facilities of the central computer and are normally 
restrained to "match, branch and loop" operations. Although many CAI programs 
are run on large, third generation digital computers which provide the necessary 
time-sharing environment, the "match, branch and loop" mode of these programs 
does not really require the facilities of a computer of such power and some 
universities have responded by installing smaller computers which handle only CAI 
work. An informal survey conducted in 1970 indicated that only a handful of U.S. 
geography departments were currently utilizing CAI on either an experimental or 
operational basis, but that many of the departments had a strong interest in future 
uses of CAI. 

It is possible to utilize the basic "match, branch and loop" format of the CAI 
program in combination with the major, computational capacity of the* machine. 
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This co.uld be done, for example,, in a CAI unit dealing with problems of spatial 
opti mizat ion which, after providing instruction in the basic elements of linear 
programming and its application to spatial problems, would require the student to 
set up a specific problem in a, hhear programming format. The CAI program would 
call upon a resident linear programming routine to solve the student s problem and 
then display the results of this computation* to the student together' with an 
appropnate commentary. Because of the multi-processing nature of large modern 
computers, the CAI program could continue the foreground conversation with the 
student while the linear programming problem was running in the background. This 
represents a rather sophisticated use of CAI, but one that will be increasingly 
utilized. There are, however, several basic disadvantages to current CAI approaches 
which arise largely out of the computing environment as it stands at the present 
time. These include: 

1) The lack of comparability of the programs. A great variety of CAI 
programming languages are available, some of which can be utilized only on 
specific machines, and sometimes only at specific installations. Other 
languages which appear to have greater comparability, particularly the BASIC 
language, suffer from a lack of standardization which inhibits the easy 
transfer of CAI programs from one installation to another. 

2) The slow printing speed ofmyst remote terminals. An inexpensive terminal 
tend^to print so slowly that 4here is difficulty in retaining student attention, 
wHile terminals with higher printing speeds are, at present , unduly expensive. 

3) The difficulty of producing any form of graphic output at a remote terminal 
This is an impediment which appears particularly critical foil geographic 
applications since many geographic CAi applications would benefit greatly 
from the output of simple graphics as well as maps. Available graphics 
terminals are currently both specialized and expensive, while pr nting a one 
page map at a Teletype terminal utilizing a printing area of, say 50 rows by 
70 columns would require the printing of 3,500 map elements/(characters) 
and would take five minutes or more. As terminals with higher printing 
speeds become available, this time will be cut to two minutes or undw, but 
even this is not adequate. At present, then, geographically based exercises 
which require more than exceedingly modest amounts of graphic output, or 
which require map output, can be carried out only on expensive remote 
terminal installations or by utilizing the batch processing capability of the 
central computer. 

It must bc^mphasized that these problems will undoubtedly be overcome in the 
next few t yeaVs by rapid increases in the level of computer and communications 
technology. It would appear that a profitable interim approach would be to develop 
a series of more complex programs based on highly standardized languages, such as 
FORTRAN, and oriented toward batch processing rather than time-sharing. 

As technology improves, these programs could be converted to the "foreground- 
background'* CAI mode mentioned earlier^ Experiments have shown that the 
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program described in tins paper, LANDUSE, can be run effectively in such a 
configuration, but to do so at present requires the use of a high-speed terminal, a 
specialized communications link and a 300 line/minute printer. This type of 
installation has a rental rate significantly in excess of Sl,000/month. Within a 
period of four to six years, it should be possible to incorporate these programs into 
a CAI environment in which such equipment configurations are less expensive and 
much more common than they are today. 

/ \ 

The LANDUSE Program ^ y 

The program described in this paper, LANDUSE, is designed as a laboratory tool 
for use in economic geography courses at both the beginning and intermediate 
levels. It is not a CAI program in the current sense of that term, since it does not 
utilize the basic "match, branch and loop" structure. It is, rather, directly related to 
the third type of instructional program termed "generative" by Huke, Fielding and 
Rummage (6). Programs of this nature permit students to enter parametric 
information into a process model and to note the fashion in which the model 
operates upon the information to produce the desired output. Students are 
expected to interpret the patterns produced with reference to the theoretical 
concepts contained in the process Vnodel. The impact of changes in the system can 
also be explored through alterations of the parameter values on successive runs of 
the program. In this manner, the program assists the students in attaining an 
understanding of some interrelationships present in complex models of spatial 
environments. A similar approach to industrial location is embodied in the 
computer model developed by Smith and Les [ 1 3 J . \ 

LANDUSE deals with a simple model of agricultural land use (the basic von 
Tliiinen model), which is relatively rich and capable of successive'complication. Jh^ 
program was developed after an examination of an early version of the program 
prepared by Rushton [II]. An entirely new program was prepared^ rather than 
modify the Rushton program, in order to provide a number of additional features 

in an efficient fashion. These include: \ 

\ 

1) replication of the basic case consisting of a single market center and simple 
concentric rings, 

2) introduction of one or more space distorting transportation routes into the 
system in a completely flexible form, and \ 

3) an ability to change the scale of analysis in a flexible btit simple fashion. 

T x 

Chapter II provides a summary pf the process theory which underlies the 
operation ol the program. The third chapter describes the structure of the program 
and some of the problems encountered in classroom use. The final chapter presents 
six selected cases which demonstrate the flexibility and scope of the<LANDUSE 
program. * s J 
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CHAPTER II 

/ 

AGRICULTURAL LAND USE: 
A BRIEF THEORETICAL STATEMENT 

The basic theory pertaining to the spatial structure of agricultural production 
originated a century and a half ago in the work of a German economist, Johann 
Heinrich von Thunen. Since then, other researchers have produced more complete 
and complex statements of this theory but we are concerned here largely with the 
simple case as originally presented in von Thiinen's classic book. The Isolated State 
[18]. In this work von Thunen attempted to discover those laws which influence 
the prices of agricultural products, and the way in which price variations are 
subsequently translated into patterns of agricultural land use. Basically, he was 
concerned with competition between various agricultural products and farming 
systems for the use of specific plots of land. A nrimary assumption of his model 
was that on each piece of land, the activity which yielded the highest net return 
from utilization of the parcel would be found .and that competing activities would 
be displaced to locations where tl^iey, in turn, would yield the highest net return. 

To facilitate discussion of thi^ problem, von ThUnen utilized a hypothetical case 
- his "isolated state" wherein he posited the existence of a market center located 
at the center of a large, fertile pjajn which containeujio transportation routes. The 
undifferentiated plain in which transfer costs are a simple function of distance, the 
"isotropic plain," appears as a starting point in many geographic theories. Von 
Thunen further assumed that throughout the plain the soil was equaj in quality and 
that the agricultural productivity of specific plots of land did not vary from one 
location to another. At some great distance from the town the plain turned into an 
uncultivated wilderness which effectively cut off all communication between the 
"isolated state" and the outside world. In this isolated community, the central 
town would supply the surrounding rural area with all manufactured products and, 
in return, would obtain all of its agricultural inputs from enterprises located on the 
plain. The problem as posed by von Thunen was. what patterji of cultivation would 
take shape, given the stated condition's, and how would farming systems in different 
portions of the plain be affected by their distance, from the central town? An 
explicit assumption in his work was that all agricultural entrepreneurs made 
decisions iri an economically rational fashion. That is, entrepreneurs responded only 
to profit maximization as an optimizing goal and were assumed to be in possession 
of complete information pertaining to all choices. 

Economic Rent and Location Rent 

As the basis of his analysis, von ThUnen introduced the concept of economic 
rent. While he developed this notion independently, it had also been introduced a 
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few years earlier by an fcnghsh economist, David Rieardo, who was interested in the 
way certain factors conditioned agricultural production decisions. He wafe, however, 
concerned primarily with variations in productivity rather than with the impact of 
puiely spatial factors. The Ricardian version of economic rent, as contested to that 
of von Thunen, is clearly set forth by Miclutel Chisholm in his book. Rural 
Setttement and Land Use (2 ) . ^ 

The concept of economic rent is* central to many questions of lapd use 
competition and it is utilised in the vori Thunen model to provide a means whereby 
this competition is resolved in a rational fashion. The observed gap bejween 
revenues an<i costs represents economic rent in a technical sense. Von Thunen 
observe^ that Ricardo had bas&d his discussion of the nature of economic rent upon 
dit^ferehces in t^ie inherent productivity of various parcels of, agricultural lanci, but 
that the same phenomena would arise'if productivity were held constant and tV 
parcels wefe subjected to varying transportation charges based upon t!*e > 
moving their production over different distances to the central market. In » « 
case, differential fertility of the plots influenced the level of gross returns,\wheieas 
in von Thunen s model, gross returns remained constant, but the cost of 
^transporttjiiton varied with distance to'thfijaarlcet. 

the two c,ases* are quite similar, bunrecenHiterature has tended to apply the 
term, 'location rents" to those economic rents arising solely out of differences in 
lot tion. This serves' to differentiate them from more general scarcity rents which 
may arise from factors such as simultaneous variations in productivity as well as 
* -location. It should be made clear, however, that location rents constitute a specially 
defined subclass of general economic or scarcity rents and are identical to the latter 
only in the case wheremjocation i$ the only factor differentiating the parcels of 
land. 

Given the assumption that we are dealing with a finite y/oup of activities (which 
/in themselves may consist of fixed bundles of subactwities, e.g., crop rotation 4 ' 
' schemes) it is then possible to view competition for the use of agricultural land as 
, equivalent to the following process, each 'entrepreneur wishing to engage in an 
activity evaluates the utility of a given parcel for carrying out a specific activity and 
determines the level of renf which can be expecteu. This figure represents the 
maximum amount that the operator will be willing to. bid for use of the parcel of 
land, since any bid higher than this estimate would require that the activity in 
question be operated at a loss. The bids for use t)f the parcel by various 
entrepreneurs desiring to undertake different activities may be ranked and the 
maximum identified. The activity associated with this maximurn bid will then. 
"Become the one which occi<pies that particular unit of land. This notion may be 
extended to*nor>agricultural activities, of course, and is roughly equivalent to the 
"highest andSest use" notion which is frequently encountered in the urban land 
economics literature. In fact, a significant proportion of the current theory dealing 
with the internal spatial structure of the city rests upon the extension of this notion 
of location rents to form a basis for land use competition in the city, e\.g., 
Alonso ( 1 ] . J 
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A Mathematical Statement^ the von Thiinen Model 

In his original work, von Thiinei^was concerned with specific cases, but he did 
produce a general statement relatingToftrtioo rent to distance from the market 
center. Typically the basic von Thiinen model can be stated as a simple linear 
equation: 

R, = YiP, - Y.C/- Y.T.D = Y.lPj - C, - t)d| Ji = I . .,n) (i) 

whf ;aMoii rent for ihc \ th activity n a specific paicel of land 

o!al ptoduction/unit of land for tlie \ th activity 
Pi - unit price at the market for-the i th activity 

unit production cost for the \ th activity 
T, = transport charge/unit/kilometer for the \ th activity 
D = airline distance of the parcel of land from the market in km. 

This equation states that the bid rent for activity / at a specific parcel is equal to 
the gross returns cn production /rom that parcel (Y,P,), lejss the total cost of 
production (Y,C,), le^s ihe cost of moving that production to the market (Y,P,D). 
The simplest situation^ obtained when we assume the existence of an isotropic 
transport plain in which T, is* linear function of distance and is assumed to have an, 
intercept of zero, i.e., there are no terminal charges in the system. It is also 
normally assumed that productivity is constant over the plain. That is, Y, does not 
vary locationally within ihe system and that a similar condition holds with respect 
to unit production costs thus yielding a constant value, C„ for every parcel. The, 
* unit price at the market is also assumed to be constant and the market will absorb, 
at that price, any amount which is offered. 

For a specific activity, these assumptions leave us with identical net returns (Y,P, 
- Y,C,) excluding transport costs for every parcel on our hypothetical plain. The 
only variations in rent are then introduced by variations in trajjsfer costs (YiT,D) 
which are L^sed entirely upon differences in the value of D. That is, when the 
parcel is located jnhe market site the t third term in equation (1) is zero and 
location rent is equal to gross returns minus total cost of production. As parcels 
~ further and further from the market are considered, the size of the third term 
XY^jD) in the equation increases. At some distance from the market transport 
charges will be equal to gross returns and the rent wilLfall to zero. This distance 
marks the limit of feasible cultivation for the activity since production beyond this 
point would be unprofitable. The term relating to unit production costs, C„ is 
defined in a somewhat unique fashion. It includes, not only th|£ normal cost of 
carrying out the activity but also includes a specific term representing a "normal" 
profit for th* operator. This means, theni that location rents can be viewed as a 
potential "excess" profit for the entrepreneur. 

Since equation (1) is linear and continuous, it is easily represented in a graphic 
foqn. In Figure 1, a bid rent curve for a specific activity has been plotted for 
parcels of land located at successively increasing distances from the market center 

er|c „ * ■ 7 . / 
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M distance — D 

Figure 1. The Bid Rent Curve Tor a Single Activity 

at M. The rent K is the net return derived from a parcel located at the market place 
and is equal to YP YC since transport costs at that location are zero. The distance 
MD represents the maximum distance at which a parcel can economically engage in 
the activity. Tor distances greater than MD either a subsidy v must be provided or it 
must be assumed that the operator is not behaving m an economically rational 
fashion. The :>lope of the bid rent ctiive is obviously equal to the unit transfer cost in 
the system and becomes increasingly steep as unit transfer costs increase. Figure 2 




Figure 2. Bid Rent Curves for Two Competing Activities 



displays a competitive situation where there are two activities, A and B, whose 
different bid rent curves result from differences in unit market price, unit 
'production costs, or unit transfer costs between the two activities. In each case, 
these values remain constant for the activity at any location, but they are permitted 
.'to vary between activities. As in the previous example. K A represents the rent at*' 
"the market for a parcel utilized by activity A, while K B represents the rent at the 
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market f|>r a parcel utilized by activity B The slopes of the two bid rent curves 
differ significantly reflecting higher unit transfer costs for A than for B If there 
were no competition present in the system and each activity were considered singly, 
then D A would represent the maximum distance at wfuch activity A could be 
undertaken with D B (D B > D A ) representing the maximum limit ol cultivation for 
activity B. However, the two activities are in competition for use of the parcels and 
within the distance MD* the bid rent of activity A always exceeds that of activity 
B, and hence activity A will then occupy all parcels of land between the market and 
location D*. However at distances greater than D* activity B produces a higher bid 
rent than activity A and activity B will occupy those parcels of land lying at 
distances from the market greater than D* and less than D B . An entrepreneur 
located on a parcel of land exactly D* kilometers from the market would be 
indifferent between the two activities since both produce identical bid rents at that 
point and no economic basis exists for differentiating between them. 

The graphic analysis of the competitive situation shown in Figure 2 is only 
one-dimensional and total rent in this case is equivalent to the area under the curve 
KaRDb. When translated to the two-dimensional case, two rings centered at M are 
formed. The first, with radius MD* enclosed the area of cultivation of activity A; 
the second, starting at distance D* and extending to radius MD B encloses the area 
of cultivation of crop B (see Figure 3). The resulting spatial pattern is one which is 




Figure 3 Simple Tinmen Rings 

commonly referred to as 'Tinmen Rings." The total rent in the two-dimensional 
case is then given by the volume of the <oiid of revolution formed by rotating 
KaRDb about K A M as an axis (see Appendix A of Dunn [3 | for a more detailed 
discussion). 

A great deal has been written on the von Tinmen model arid for discussions of its 
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many complkationsStfd extensions, see the works by Dunn [3), Isard [7], and 
Alonso [I]. The utility of even the simple model in explaining the observed 
patterns of land use has been noted by a number of researchers and much of this 
latter work is admirably reviewed in Chisholm [2] . 

Extensions of the Simple Model 

The basic von Thiinen model outlined above displays a number of important 
factors relating to agricultural production location decisions, and does so in an 
exceedingly simple fashion. While it may be characterized as one of the simplest 
models' of land use, it is exceedingly fertile since it is capable of Extension and 
modification in a variety of ways. Two of these, the expansion to'mfcltiple market 
centers and the introduction of space distorting transport routes will be discussed 
here. 

Multiple Market Centers 

Given the concept of land use competition, the extension of the basic von 
Thiinen model to the case of multiple urban centers is relatively straightforward. 
Basically, equatipn(l) is recast into the following form: 

Rij = Y.Pjj - YiC, - Y^jD, (i = 1 n; j = 1 , . . .,m) (2) 

where. R,j = rent produced by the t th activity when it is oriented toward the j //f 
"market center 
P,j unit price for activity i at the ) th market center 
T,j = transport rate /unit/km. for movement of output of activity i to 
market center j 

(the remaining terms arc as previously defined) 

Utilization of this formulation results in the establishment of multiple bi'd rent 
curves based upon the vanous market centers (see Figure 4). The structures arc as 
defined in thv previous example except that, in this case, an evaluation of the bid 
rent curves shows that at certain locations, for example, L, three "different" 
activities arc bidding for the use ot the land. At L we sec positive bid rents derived 
from activity A oriented toward market 1 and activity B oriented toward market 1, 
as well as a bid rent derived from activity B oriented toward market center 2. At L, 
activity A oriented toward market center I produces the highest bid rent and hence 
occupies that location. However, at distances greater than M,P we sec that 
orientations of the two activities, A and B, toward the second market center 
produce higher bid rent levels than a continued orientation toward the original 
market center. Thus, at distance M,L* from the original market center, it can be 
seen that while activity B oriented toward market center 1 is still making a positive 
bid for use ot the land, this bid ts exceeded by those made by activity A and activity 
B when oriented toward the second market center. In this particular case, activity 
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Figure 4, Bid Rent Curves for Competition Between Two Market Centers 



A, with an orientation toward the second market center, occupies the location. 
Thus, along the line from M| to M 2 , we see a succession of land uses, first A then B 
and then A again However, within the area devoted to the cultivation of land use B 
there is a breaking ooint, P, at which the land shifts from the tributary area of the 
first market center to the tributary area of the second center. 

An interesting problem involves the exact determination of the boundary 
between the tributary a\eas of the two centers. The boundary is clearly defined in 
the one-dimensional case v shown in Figure 4, and can be easily inferred in the 
two-dimensional case where the parameter values for the two market centers are 
identical and an overlap between the bid rent curves exists, the boundary in this 
latcer case wili be a straight line perpendicular to the line joining the two centers. 
However, as one of the centers becomes "dominant," the boundary between the 
tributary areas of the two centers will gradually shift from a straight line to a 
hyperbola to a parabola and finally to a closed elliptical form around the second 
area with the second (smaller) market center occupying the foci of the ellipse 
nearest to the larger market center. Demonstration of the development of the 
spatial structure in fhese cases is usually very difficult without resort to 
complicated mathematical derivations or an elaborate, carefully plotted graphical 
display. It is, nevertheless, an interesting extension of the basic model. 



System-Wide Transport Improvements 

The introduction of system-wide transport improvements, such as those that 
might occur when a new and more efficient type of vehicle is introduced into all 
parts of the transportation system can be handled within the model as given in 
equations (1) and (2). The net effect is simply to reduce the value of T, or T u , 
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which, in turn, leads to a reduction in the slope of the bid rent curves. The benefits 
of any system-wide improvement in transportation may be differentially distiibuted 
between activities since, for example, the introduction of a specialized and more 
efficient bulk carrier might favor certain commodities over others. Another 
example of differential impact would be the introduction of refrigerated, milk 
transports which would obviously influence T, for the activity of milk production, 
but would have little direct impact up&n the transport costs for commodities such 
as grain. 

Route-Specific Transport Improvements 

A far more interesting, and theoretically complex, case is encountered when 
transport improvements in the system differentially influence locations. This is 
typically caused in the real world by the coristruction or improvement of an 
existing transportation link. For example, a n^w highway linking two nodes in the 
system may be built or an existing road may be improved through resurfacing, thus 
leading to lower vehicle opeiating 'costs. Within the context of the basic von 
Thdnen model, the introduction of specific transport routes can be viewed as a 
distortion of the uniform plain since unit transfer costs over the improved routes 
will generally be lower than transfer costs over the remainder of the area. Within 
the context of equation (I), this change can be handled in one of two ways: 
through variations in the values of the T, depending upon the type of route over 
which the travel takes place or, somewhat more easily, by assuming that variations 
.occur in D itself. This latter approach arises out of a substitution of the notion of 
effective distance for that of airline distance which was previously utilized in the 
model. Effective distance is defined as follows: let D c be the minimum distance 
from the cell to the highway, D t the minimum distance from the market center to 
the highway, and D\ the distance travelled along the highway. Then the effective 
distance D is given by the formula: 



with a in this case representing the ratio of on-highway to off-highway movement 
charges within the range 0 < a < I . That is, for an a value of 0.5 travelling a unit 
, distance along the highway would cost exactly half as much as travelling the same 
distance off the highway. The use of a as a scaling factor in the equation for 
effective distance compensates for the real differences in unit movement costs over 
the different routes. This proves to be particularly useful 'when the definition of 
effective distance is expanded to include travel over a potential system of routes: 



De = D c + a D, + D t 



(3) 



(4) 



where there are a series of r. different highway segments in the route, each of which 
has associated with it an unique value of a. ' \ 
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The distortfng effect of low cost transportation routes was noted by von Thllnen 
in his original analysis in a discussion of the distorting effects of a river which 
would provide low cost transportation for agricultural commodities upon the limits 
of the area under cultivation. Derivation of the geometric impact of highway 
improvement upon land use patterns has been undertaken by several authors, but 
the discussion that follows is adapted from that of Werner [20). Werner was 
concerned with the distortion of the shape of the hinterland area that occurred 
after a new transportation route had been constructed. He defined hinterland as the 
area accessible from a given point with a specific level of transportation effort, 
where effort may be measured in terms of energy consumption, mileage, time or 
other criteria. The analysis assu) es as given the isotropic transport plain, a point P 
corresponding to a market cem?r and a straight transportation route running 
through the area witha cost per unit flow per kilometer of a, on the road and a 2 
off the road. This basic situation is illustrated in Figure 5. 

In Figure 5, for convenience," the new transportation route is assumed to be 
identical with the x-axis and P is a point on the y-axis with coordinates (0,b). The 
situation as shown is a symmetric one so that the analysis need only be completed 
for points lying to the right of the y-axis. The minimum co$t connection from an 
arbitrary point Q 2 located sufficiently far from P and lying off the transportation 
route is given by a direct movement from Q 2 towards the road entering with an 
angle equal to arc cos a 2 /a,, then along the road itself to a point R, and then 
directly to point P. The angle at which the minimum cost path leaves the road is 
again equal to arc cos a 2 /a! . . 

Werner demonstrates in his analysis that the half space under consideration can 
be divided into three regions, which are labeled I, II, and HI in Figure 5. Within 
region I all points have a minimum effective cost path which consists of a direct 
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Figure 5 Determination of the Sphere of Influence of a Transport Route 
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movement to point P with no movement over the road. For all points in region 11, 
such as the arbitrary point Q 2 , all minimum cost connections involve travel over the 
road in the fashion outlined above. The situation in region 111 is more complicated 
since for points arbitrarily close to the road, minimum effective cost paths will 
invojve a combination of direct travel and travel over the road. For points much 
further from the road, direct travel will be more efficient. The boundary separating 
the two parts of region 111, is shown to be a parabola with the point P lying in one 
of its foci. The resulting situation is summarized in Figure 6 which displays the 
areas of access via the new road versus the areas which retain direct access to point 
P. The shaded portion in Figure 6 represents those points in the hinterland of P tliat 




new road 
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Figure 6. The Sphere of Influence of a Transport Route 
Source: Werner [20) 



have had their distance relationship changed by the introduction of the new 
transport facility , while points lying within the area of direct access retain the same 
distance relationships as before. For those points in the shaded area, the effective 
distance has been reduced by the introduction of the highway and this means that 
they are, in a purely economic sense, closer to P than they were prior to the *' 
introduction ot the highway. By this transformation the areas of cultivation will be # 
extended along the road, and in the areas immediately adjacent to it, in the fashion 
shown in ^Figure 7 which displays a simple Thuncn ring structure for both the 
before and after cases. Werner goes on to show how the patterns of changed a cess* 
are influenced when a variety of pre-existing transport networks are assumed and 
the interested reader is referred to his original exposition for further details [20| . 
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Figure 7. Ring Distortion Produced by a Transport Route 
Source: Werner [20 1 



Summary 

Materia! presented in this chapter has provided a summary of basic aspects of the 
von Thunen model and some of its more direct extensions. The reader interested in 
a more detailed exploration of these problems is referred to the works by Alonso 
[1] , Chishobh (2] , Dunn (3] , Isard [7] and Losch [9] . 
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CHAPTER I II- 



THE LANDUSE PROGRAM 

i 

The objectives in using computer programs such as the one discussed here are set 
forth in the article by Huke, Fielding and Rushton [6] , and in the report of the 
Engineering Education Project of the University of Michigan [17], The selection of 
- this specific process model was based upon several factors: 

1) the von Thunen model is perhaps the simplest generative model for spatial 
patterns of land use, but is rich enough to permit a number of interesting 
extensions, 

2) the basic model is widely known and is discussed in most economic 
geography courses, and 

3) although the model is simple in structure, a number of operational 
difficulties frequently arise among students in visualizing jhe changes in 
spatial structure resulting from parameter shifts in the model. 

The program was developed to provide a means whereby the student could 
develop a stronger feeling for the' relationship between the abstract model and the 
spatial structure which it represents through manipulation of selected portions of 
the model and observation of the subsequent impact upon the spatial structure of 
land uses. The conventional "blackboard" approach utilizing one-dimensional 
cross-sections (such as those in Figures 1, 2 and 4) has not been proven to be an 
effective solution to this problem. Secondary development objectives included a 
desire to acquaint the student with the complex spatial structures that can arise out 
of even simple linear models (e.g., case V in the following chapter), and to provide a 
large group of undergraduates with an introduction to the use of the high-speed, 
digital computer as an analysis tool. 

Program Structure t><^ s 

The program presented here has been designed with flexibility of use in mind. It 
has the capability of handling exercises ranging from those that might be 
appropriate for an introductory course emphasizing simple concepts to those 
suitable for advanced courses dealing with more complex probjems. Examples of 
the former include generation of simple Thunen rings and the demonstration of the 
development of a competitive boundary between two market tenters. The latter 
could involve any of the many complex 'cases of transport improvement and 
. development. In all cases, the model relies upon the basic von Thunen model for 
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evaluation of competition among land uses. Competing market centers and crops 
are readily handled within this framework, and one or more distance distorting 
transport routes can be introduce^ through the concept of effective distance 
discussed in the previous chapter*. 

The program assumes the existence of an unbounded, isotropic plain upon which 
an arbitrary (x, y) coordinate systerii has been established. It then isolates a 
rectangular subregion from the plain and proceeds to develop the pattern of land 
use according to the basic rent equation. Unless otherwise instructed, the program 
assumes that the rectangular subregion is one whose lower left hand corner 
corresponds to the zero points of the coordinate system* For pragmatic reasons, 
related to both ease of computation and ease of output, the program assumes that 
this rectangular area is divided into a series of small areal units which form the basis 
for all computations. There are MO of these units along the horizontal axis and 90 
of the units along the vertical or y-axis, producing a total of slightly under 10,000 
unit areas. All coordinates which locate features on the plain, such as market 
centers*' are specified as integer numbers in a read right-and-up coordinate system* 
These integer numbers correspond to cell numbers. For example, the coordinates of 
a market center located approximately in the center of the rectangular region 
would be given as x = 55 and y = 45, Within the unbounded plain, the student Is 
permitted to locate up to a maximum of ten market centers through the simple 
process of specifying their x and y coordinates. However, only those towns lying 
within the rectangular area defined by I < x < 1 1 0 and 1 < y < 90 will appear on 
the output maps. 

The student also specifies a >enes of parameters which define up to twelve 
abstract crops, with each crop specified in terms of the following parameters: 

1) Yield per unit area (assumed to remain constant for that crop over the entire 
plain); 

2) Unit price at eaoh of the market centers; 

3) Unit cost of production (assumed also to be invariant over the plain); 

4) Unit cost of transport per unit distance (permitted to vary between market 
centers and between crops). 

With this information in hand, the program proceeds to evaluate the rent 
equation for each unit area in turn beginning with the unit in^he upper left hand 
corner. The distance from tins unit to each of the definctl market centers is 
calculated and the rent, according to the basic rent equation, that would be derived 
at tins parcel from each crop and market center combination is obtained- For 
example, in a system that had two market centers and two crops defined, each unit 
area would have four possible bid rents based upon, crop A at market center 1 % 
crop B at market center I , crop A at market center 2 and crop B at market center 
2, From this vector of four bid rent values, the program selects the maximum, and 
thus identifies the particular land use which is assigned to the parcel knd determines 
the assignment of the parcel to a specific tributary area. If all bid rents at that 
location are negative, then the use assigned is "idle" and the parcel is not assigned 
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to the tributary area of any urban center. For positive bid rents, the maximum 
value is also stored so that subsequent summaries of total rents and total areas 
under cultivation can be provided by crop and by crop-market center combinations. 

Transport routes are defined by the simple expedient of providing the x and y 
coordinates of their points of origin and termination together with a value for a, 
the ratio of on-highway to off-highway movement charges. Program operation 
proceeds as described before except that the minimum effective distance is 
established for the computational cell prior to evaluation of the rent equation. The 
information on minimum effective distance, when compared with direct access 
distances, permits definition of the spKere of influence of the transportation route. 

Subroutines 

The program is written entirely in FORTRAN IV and consists of a main routine 
and two specialized subroutines. The main routine, LANDUSE, is responsible for all 
input and computational phases. The only exception occurs when highways are 
present and the subroutine PROJCT is called to determine the shortest distance 
between the analysis cell urider consideration and" a specific hjghway. It will be 
recalled that Werner [20] demonstrated that entry into the transport route should 
take place at an angle which is a function of the ratio of the two costs. PROJCT, 
for reasons of computational simplicity, assumes that entry alwa>s takes place at 90 
degrees. This introduces some distortion, but it is felt that this can be safely 
disregarded in this particular case. 

The subroutine REGION is responsible for all map output- It is designed to store 
and output up to three one-page maps using the smallest possible amount of core 
memory. In order to fall into the most favorable level of the priority system, the 
program is written to load and go using only 14,000 words oi core. Storage of the 
three 54 x 110 maps would normally require nearly 18,000 words of computer 
memory, but REGION reduces this storage requirement to under 1,700 words 
through the use of masking statements. The amount of storage required b,y 
REGION is dependent upon the computer word length and since word length for 
the CDC 6000 series machines is significantly greater than for comparable 
computers, such as the IBM 360, the CDC version of REGION will not dpfej^Je on,^g£ 
other machines. Versions of REGION have been prepared, however, for use on the 
IBM 360 as well as on CDC 3000 series machines, both of which have different* 
word lengths than the CDC 6000 series. These modified routines produce output 
maps of slightly different sizes than those described here (the 110 units in the 
x-dimension of the analysis area represents eleven CDC words each holding ten 
characters), but the general structure is identical. / 

Production of maps on the high-speed printer requires that all output be scaled 
by a factor of 3/5 to produce identical vertical and horizontal scales. This results 
from the size of the characters on the pnnt chain (10 per inch)^and the spacing 
between lines (6 per inch). Tins is why LANDUSE uses a 90 x J 10 analysis region 
and produces 54 x 110 maps If this scaling were not clone, the Thiinen rings would 
be significantly distorted. 
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Program Limits ~ 

The program allows maxima of ten market centers, twelve abst racf crops and 
nine transport routes. The number of computations required in the maximum case 
can be quite large.\For example, m a problem with ten market centers and twelve 
crops, the program must evaluate 120 rent values for each of nearly 10,000 cells. 
The problem of determining patterns of land use in a system involving a relatively 
large nuiuher of market centers and crops is not nearly as severe, however, as the 
case where there are a large number of transportation routes with multiple 
intersections. Since each transportation route is permitted to have an individual 
value of a the determination of the minimum effective distance becomes quite 
complex. It is possible, therefore, within the defined maxima of the program to 
establish problems winch cannot be economically solved even with a 'large 
high-speed computer. These rather high limits have been established to provide the 
program with the maximum amount of flexibility and not with the view that all 
will be reached at the same time. From a pragmatic standpoint, it is desirable to 
have a lower set of limits for student^use. Within the program, these lower limits are 
established at four market-centers, four crops and two highways and this is the 
configuration which is utilized in the normal classroom situation., Access to the 
larger capabilities of the program may be obtained by punching a single decoding 
"password" on one of the problem control cards. If this password is absent, then 
the lower set of maxunas are in operation; if present, access is permitted to the full 
capabilities of the program. 

Data Input ' 

* Data input to the program takes place via a sequence of control cards, the first 
of which is a general problem control card specifying the number of market centers, 
the number of crops, the number of highways and the number of re-scale runs that 
are to be made. This card also carries basic title information which is reproduced to 
ider*ify the output, as well as the password if one is entered. The re-scale option is 
discussed in detail later, but basically a problem may be re-scaled through the 
specification of a new origin and a new scale factor for the map. If a problem is to 
be re-scaled, an indication of intent must be giveh on the problem control card. 
Following the problem control card, one or more datacards giving. the market 
center coordinates are read and are followed by a third group of cards containing 
information on the desired scaie of the map (given in terms of tl.e number of 
distance units per inch) a$ well as yield levels for each of the crops. The fourth 
group of cards is called crop parameter cards, and specifies for each crop and market 
combination: the unit market price, the unit production cost, and the unit 
transportation cost per unit of listance. If transportation routes are specified in the 
problem, then a fifth .group of cards follows specifying the starting and ending 
point of each transportation route as well as its associated a value If the pioblem is 
to be re-scaled one or jnore times, then the final group of cards consists of re-scale 
cards specifying x and y coordinates of the new origin as well as the new distance 

o 20 



^ < H3D026 



scale factor. *ny number of problems can be run in sequence by repeating the 
above deck structure K 

Detailed instructions on the preparation of input cards are given in Appendix A. 

* . 4 

Program Output 

* The output of the program consists of a series of maps as well as basic summary 
statistics. The first output stage consists of an echo check of data inputs and 
repeats, in tabular form, the information pertaining to market centers, crops and 
-"•""highways that has been read from the data cards (see Figures 8 and 9). This permits 
the student to check quickly for errors in his control cards. Following the echo 
<?heck, from one to four maps are produced depending upon the problem control 

J card parameters. The first jnap is always a general locatic s map which displays the 
location »of air market centers and highways in the problem, (see Figure 10). All 
maps are identical in format with market center locations shown as "X" *s, and all 
carry small cesses at the corners of the rectangular analysis area to serve as 
registration points should the student desire to overlay one or. more maps for 
comparative purposes. If highways have been defined in the system, then the 
second map is one showing the location of the market centers as well as those cells 
rf which fall within th£**spherer of influence" of each of the highways (see Figure 11). 
Highways <are identified by the program by reference to the sequence in which 
the highway control, cards are read. That is, the first card is identified by the 
progran} as being associated with highway number 1, the second card with highway 
number 2, and so forth. A cell whose effective distance to any of the market 

^ centers is reduced through the existence of highway 1 has a "1" printed in it. 
Similarly, any cell which is influenced by highway 2 ha's a "2" printed in it, etc. 
Those cells through which direct access still represents the minimuu effective 

9 distance are left blank. Thus, the area of the map which is filled with "1" *s 
represents the sphere of influence of highway 1, the area filled with "2" *s 
t represents the sphere of influence of highway 2, and so forth. 
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Figure 1 i. Output Sample-Sphere of Influence of Highway 




The next map produced is a general land use map (see* Figure 12). This map 
slxows the location of each of the market centeis as well as the pattern of land use. 
The first control card read is associated with crop A, the second with B, etc. Any 
cell for which crop A has the highest bid rent has an A printed in it on this map, 
any cell with crop B having the highest bid rent has a B printed in it, etc. Cells for 
which all bid rents are negative are left blank. For those cells for which two or more 
bid rents are identical an "= M sign is printed indicating that a tie exists between 
two or more land uses. The last map identifies the tributary areas of each market 
center (see Figure 13). The first set of coordinates read is assumed to define market 
center 1, the second set of coordinates market center 2, etc. A cell is assigned to 
center 1 if the highest bid rent for that cell is one which is associated with sales at 
market center 1, it is assigned to center^ 2 if the sales are associated with market 
center 2. ejc. Ties are again identified by printing of an sign in the cell. 

The final output section consists of a set of summary statistics relating to rent 
and area under cultivation for each market center and crop combination (see Figure 
14). Total rent and total area under cujtivation are also given for each crop and the 
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Figure 12 Output Sample- .Land Use ^ 



total area with zero rent level, that is, idle land, is also presented. The summary 
statistics are useful in providing quantitative measures of impacts due to parametric , 
changes, , * 

Intrcduciic 1 of Transportation Routes 

Transportation routes are defined within the study area through specification of 
the coordinates of their points of origm and termination as well as the associated a 
value. The a value represents the ratio of on-highway to off-highway movement 
charges and an a value of one implies that t'aese are identical. Thus, an a value^of 
0.5 implies that travel on the transportation route is onjy one-half as costly as travel 
off the highway. Based upon empirical investigations of vehicle operating costs, it is , 
possible to select a values which ,are' at least roughly in accordance with reality. A 
number of studies have examined the effect of transport improvements on vehicle 
operating costs and if the assumption is made that an a value of 1 is roughly 
equivalent to the crudest of dirt roads, then relative values of a may be established 
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Figure 13. Output Sample Tributary Areas 



for other common road surfaces as well as for different types of vehicles. Table i 
presents some suggested values for a which are based upon empirical investigations. 
These values are only suggestive since many other factors influence the level of 
operating costs. However the fact that they are based upon "reahworld" data tends 
to make them more viable in the mind of the student. 

The sphere of influence of a highway is defined as the set of all cells for which 
- the average effective distance from the cell to all N market centers is less than the 
average airline distance to the towns, or the average effective distance to the towns 
using all other highways. A cell which falls within the sphere of influence of a given 
• highway is identified on the map with that highway's number, while those cells 
whose minimum average distance is the average airline distance are left blank. If the 
average effective distances are equal for any two highways with respect to a given 
cell, an sign is then printed in the cell. 

A transport route may be of any length and may run in any direction, but care 
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TOTAL AQCA 



A 756363.0 1258.333 

A 990586.7 <.<.06.667 



AR£A »ITm ZfHO «£NT l.£V£L ■ *231.667 

Figure 14. Output Sample -Summary Statistics 

must be tanen that, either through initial specification or subsequent re-scaling, the 
effective length of the highway is not reduced to zero. If this should occur then an 
error message is issued and the program terminates. 

Changes in Origin and Scale 

All points on the isotropic plain are assumed to be specified in terms of integer 
coordinates so that a unit increment in either the horizontal or vertical direction 
corresponds to a distance of one-tenth of an inch on the printed map. The factor 
SCALE is a non-negative, floating point number which assigns a physical distance to 
this unit increment. The physical distance may be either abstractly defined or in 
feet, miles, kilometers, etc., and ten times the value of SCALE is the distance per 
unit' of the printed output maps. The distance scale which is in effect for a 
particular problem is printed at the bottom of each map. In order to obtain maps 
covering different regions within the plain, the coordinates of the defined market 
centers and highways must be modified and/or the value of SCALE must be 
changed leaving the remaining data the same. In this case, physical distances must 
1 be regarded as invariant, that is modifications of the value of SCALE or shifting of 
the origin must necessarily cause a consistent modification of the originally 
specified integer coordinates which wHl preserve the "actual" distances in a spatial 
system. Appendix A contains information on how these scale shifts and shifts of 
the origin are computed by the program. 
.The program contains an option which will automatically re-scale any specified 
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Table 1 

OPERATING COSTS AND ALPHA VALUES FOR ROAD SURFACE 
AND EQUIPMENT CLASSES IN A DEVELOPING NATION* 



Equipment 
Class 


Road Surface 


Earth 


Gravel 


Paved 


Small passenger car 


3.26 




1 QK 

1 .70 




(1 0) 




[v.Ol ) 


Large passenger car 


5.32 


4.13 


3.23 


- 


(1.0) 


(6.78) 


(0.61) 


1-ton truck 


6.79 


4.78 


3.56 




(1.0) 


(0.70) 


(0.52) 


3.5-ton truck 


11.79 


7.47 


5.13 




'(1.0) 


(0.63) 


(0.44) 


15-ton tractor 


20.53 


12.19 


7.92 


semi-trailer 










(1.0) 


(0.59) 


(0.38) 



♦Figures are cents/vehicle-kiiometer; bracketed numbers are alpha values assuming earth 
surface to be equivalent to a = 1.0 

Source Adapted from figures given in United Nations, Transport Modes and Technologies for 
Development. 1970. , 



problem. A problem which is specified by a data set immediately following a 
problem control card is called the reference case. Each re-scale data card which 
follows a reference case may contain the coordinates of a new origin as well as a 
new scale factor. The coordinates of a new origin on the re-scale card are assumed 
to be given irt terms of the original coordinate system. That is, the system may be 
viewed by the student as taking the original rectangular analysis area (which the 
program assumes to be located with its lower left hand corner at the zero point of 
the coordinate system) and moving its lower left hand corner to the new 
coordinates specified on the re-scale card. If a scale change is also specified, the 
relative size of the region and the analysis areas will be changed as well. The original 
coordinates of the market centers and highways, as specified in the reference case, 
are saved and new coordinates are computed relative to this reference problem for 
each subsequent re-scaling. Changes of origin may or may not be accompanied by 
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subsequent changes in scale. In many complex problems, it has been found that 
changes in origin together with simultaneous changes in scale produce useful results 
in terms of clarifying the analysis. For exarnple, a complex case may be developed 
and the student may desire to examine in greater detail some limited portion of it. 
He may then request that the Scale of the map be shifted to provide more detail and 
that the origin be translated to such a location that the area of greatest interest is 
centered in the new map. 

In the case of re-scale operations, the number of analysis cells defined in the 
rectangular region remains the same while their area is effectively changed. In the 
present version of the program, this means that the summary statistics in the 
reference case and the re-scaled case will not be comparable. This problem will be 
eliminated in a future version of the program. 

Using LANDUSE 

The LANDUSE program was developed at the Department of Geography, 
^Northwestern University, over a three year period. The program was initially used in 
an upper division course in economic geography having an enrollment of about 25 
students, most of whom were graduate students from other disciplines. Problems 
were assigned on an individual basis and were fairly complex in their structure since 
the student group possessed a substantial background in related areas such as 
transportation engineering. The use of the program with this small group of mature 
students posed a few problems. 

During the next year, with the cooperation of Dr. Ian Hamilton (Department of 
Geography, London School of Economics) who was a visiting professor at that 
time, the use of the program was extended to a larger, lower division course in 
introductory economic geography. This course had a registration of about 70 
students and because of the larger class si/A the students were "teamed" in groups 
of two or three in order to reduce operating costs and the load on rather limited 
keypunch facilities. Teaming in groups greater than three did not prove successful 
because of the increasing remoteness of the individual student.from the exercise. 

At the present time, the program is used in both the lower division and upper 
division economic geography courses. The types of problems assigned differ 
significantly with the more complex problems in the upper division course 
reflecting the more substantial background of these students, and their specific 
interest in questions of transportation and its relationship to the spatial structure of 
land use. 

General experience with the program indicates that \i has a high level of student 
acceptance and that it does increase their level of knowledge in the desired areas. 



Necessary Student Background ' 

Aside from a theoretical grounding in the notion of economic rents and location 
rents, and in the general structure of the simple von Thunen model, the student 
must know how to prepare the data car's for the program, how to submit the 

ERIC . . 28 

t)0034 



program to the computer and how to retrieve his output. In most cases, particularly 
when dealing with lower division students, the instiuctor cannot assume student 
familiarity with either computers or the operation of the keypunch. Our experience 
has shown that basic keypunch instruction can be given in about 15 or 20 minutes 
and that an introduction to the program, card layout and computer center 
operation take^ about an hour. In many universities, basic keypunch instruction is 
provided by the computer center through e/ther direct demonstration or by student 
access to video recordings, A sample of the keypunch reference material given to 
Northwestern students can be seen in Appendix C. 

The first exercise which is normally assigned at both class levels is a "set piece" 
wherein all parameter values are given. This is designed to test the student's 
familiarity with card preparation and system operation. Case Pin the following 
chapter is usually used for this purpose. 



Operating Costs 

The cost of operating the program is a function of the charge level established by , 
the local computer center, as wdfl as of the number and complexity of the problems 
which are assigned by the instructor. In the introductory course at Northwestern, 
where only a few simple problems are assigned to teams of students, operating costs 
are approximately $14.00 to S 15 00 per student. More complex problems, although 
not a greater number, are assigned in the upper division course on an individual 
basis and here the cost per student rises to between $50.00 and $70.00. In many 
cases these charges, particularly those for the lower division class, reflect a level 
which is related only to the number of job submissions made, since at Northwestern 
a minimum charge of Sl.OO per job is set by the computing center. In dealing with 
simple problems- one of the major factors influencing total class cost is the number 
of submissions per student. This can be minimized by making use of the provision 
within the program for stacking problems so that a number of different cases can be 
run on the same submission. When used, this tends to reduce significantly the cost 
per student. , 

The cost per student also increased when students substitute repetitive tests on 
the computer for thinking about the operation of the model. Prior to being given 
the exercises, the students are instructed to produce trial solutions utilizing 
standard graphic methods, such as those shown in Figures 2 and 4. However, the 
student finds that it is less work to change merely one or two numbers in his input 
deck and resubmit the problem. This substantially increases the number of runs and 
the total tost. Some computer installations have monitor routines available that will 
check the number of submissions made by an individual student and will either 
prohibit access after a certain number of tries or will call the situation to the 
attention of the instructor. Such a program is not available at Northwestern, but a 
comprehensive monitoring program is being developed winch should alleviate this 
problem and permit th* 1 instructor to identify those students who are having 
difficulties with the operation of the program, 
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While the dollar figures cited here refer to a $7.00/minute computer charge, it 
must be kept in mind that at many schools the assignment of dollar values to 
classroom use of the computer does not result in any actual transfer of funds. The 
"charge" in these cases provides only a means of allocating a fixed resource (which 
is paid for in the school's annual budget) among competing users. Charge policies 
vary widely amon^ schools and the prosepctive user should make local inquiries. 

Program Availability 

Because LANDUSE is a large, complex routine which is undergoing continuing 
development, a listing of the source deck is not included in this technical paper 
However, machine- readable versions of the current distribution version may be 
obtained from the Geography Program Exchange. The activities of the Exchange 
and methods of obtaining programs are described in Appendix The GPE will 
provrde either CDC or IBM versions of the program and persons requesting 
machine-readable copies should specify the machine which will be used. 

Persons using the program are requested to forward their comments to the 
author 
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CHAPTER IV 



SOME APPLICATIONS OF THE LANDUSE PROGRAM 



The following six cases have been chosen to demonstrate the versatility and 
flexibility of the program. Parameter values were selected to reflect this goal. 
Exercises featuring more "realistic" parameters can be developed easily within the 
context of a particular course. The cases presented here begin with a simple 
situation which reproduces the basic ThUnen ring structure and then proceed 
< through sucessive complications to a complex case of transport impact combined 
with a demonstration of the use of.the re-scale option to provide a finer grain 
: ; examination of the results. 

For ease of interpretation, boundaries and market center locations have been 
emphasized through the addition of hand inserted symbols in most illustrations in 
this chapter. 



Case I - Simple Thunen Rings 

This case demonstrates the ability of the model to replicate the formation of 
simple Thunen rings around a single market center. It is normally utilized as a "set 
piece" in order to familiarize the student with card preparation, program 
submission and interpretation of the output. The parameter values for generation of 
this case are shown in Table 2. and Figures 15 and 16 display the resulting land use 
pattern and market area structure. 
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Table 2 

CASE I - PARAMETER VALUES 



Market Center Locations 



Market Center 
1 



X Coordinate 
55~ 

Crop Yields 




Crop A 
CropB 



30.00 Units 
30.00 Units 



Market/Crop Parameters 



Market 
Center 

1 
1 



Cror) 

- A 
B 



Market 
Price 

35.00 
30.00 



Production 
Cost 

5.00 
5.00 



Transportation 
Cost 

1.50 
1.00 
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Figure 16. Case l~ Tributary Area 
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Case II - Competition Between Market CerHeW 

The case of competition between two market centers can be dealt with easily 
within the context 'of the program. Case II represents a somewhat arbitrary 
selection of parameters leading to a situation where the first center is "dominant" 
with respect to the second. The land use pattern which results is displayed in Figure 
17. An examination of this map shows that Thiinen rings develop around both 
centers and that land use B merges between them. A more interesting pattern is 
revealed when the tributary areas map, Figure 18, is examined. Here the dominance 
of the first center is clearly evident and the boundary between the trade areas of 
the two centers bends in the direction of the second center which occupies an 
off-center location within its trade area. Through suitable adjustment of the 
parameters, the boundary between the two trade areas may be bent even further to 
the right and in the Extreme, caused to close on the far right thus embedding the 
trade area of center 2 entirely within the trade area of center I. Table 3 gives the 
parameter values which will reproduce this case. 

' i 

Table 3 

CASE II - PARAMETER VALUES 
Market Center Locations 



Market Center X Coordinate Y Coordinate 

1 36 45 

2 74- 45 



Crop Yields 



Ciop A 30.00 Units 

Crop B 30 00 Units , 



Market/Crop Parameters 



Market 
Center 

I 

1 

2 

r> 



Crop 

A 
B 
A 
B 



Market 
Price 

45.00 
40.00 
35.00 
30.00 



Production 
Cost 

' 5.00 
5.00 
5.00 
5.00 



Transportation 
Cost 

1.00 
0.75 
1.50 
1.00 
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Figure 17. Case II Land Use 
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Figure 18. Case II -Tributary Areas 
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Case III - Transport Route Bypassing a Market Center 

This case is a replication of case 1 with the introduction of an^east-west 
transportation route running slightly to the south of the market center as shown in 
Figure 19. Figure 20 gives the sphere of influence of the transportation route with 
those cells containing a "1 M having a minimum effective distance which involves use 
of the transportation route. The thin vertical 1 line of Ts immediately below the 
market center results from an anomalous condition wherein crossing of the 
transport route is counted as use of the route. This will be removed in a future 
version of the program. Figure 21 demonstrates the resultant patterns of land use 
and should be compared directly with Figure 15. Although there is an observed 
extension of land use B along the transportation route, there is no impact upon the 
pattern of cultivation of land use A, such as one might have expected from an 
examination of Figure 7. The reason for the lack of distortion of the A rings can be 
seen through a comparison of Figures 20 and 21. Superimposition of Figure 20 over 
Figure 21 reveals that the sphere of influence of the highway does not extend into 
the A ring and hence no change is recorded for those cells. Figure 22 displays the 
resulting distortion of the trade area along the transport route and Table 4 gives the 
parameter values which will reproduce this case. 



Figure 19. Case III- Market Center and Highw. jcation 
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Figure 20. Case III Spnere of Influence of Highway 
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Figure 21. Case I II -Land Use 
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Figure 22. Case ill -Tributary Area 
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Table 4 

CASE III - PARAMETER VALUES 



Market Center 1 ocations 



Market Center 
1 



Market 
Center 

1 
1 



Crop 

A 
B 



X Coordinate 



55 

Crop Yields 



Crop A 
Crop B 



30.00 Units 
30.00 Units 



Market/Crop ?aiameters • 



Market 
Price 

35.00 
30.00 



Production 
/ Cost 



5.00 
5.00 



Highway Data 



Y Coordinate 
• 45 



Transportation 
Cost 

1.50 
1.00 



Starting Point 



EnJing Point 



Highway X Coordinate Y Coordinate X Coordinate Y Coordinate Alpha 



39 



108 



39 



0.45 
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Case IV - Market Center Competition and Transport Impact 

This case is a replication of case IT with the addition of a transport route 
bypassing the two competing centers a short distance to the south. The spatial 
positioning of the two market centers and the transport route is shown in Figure 
j 23. Because of the existence of twoicenters, a more complex pattern emerges when 
the sphere of influence of the road is examined (see Figure 24). The intact of the 
transport route produces an extension of cultivation of both crops A and B to a 
point where the new band of cultivation of crop A along the highway r arly joins 
the two centers (see, Figure 25). There is, however, still a smal 1 ,iand of B 
production remaining at the edge of the A-B ring around the secor enter. Again, 
many of the more interesting results of the insertion of the transr .ation route can 
be seen in the examination of Figure 26 w' n displays the .ifiguration of the 
two market areas. Some very significant changes can be „rved by comparing 
Figure 26 with Figure 18. The introduction of the transport route has extended the 
trade area of center 1 far to the east of center 2 and the lower boundary between 
the two trade areas has now b^en wrr/ped up and around center 2, so that its trade 
area has become a parabola oriented toward the northeast. The parameter inputs 
which replicate this case are given in Table 5., 



Figure 23. Case IV-Location of Market Centers and Highway 
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Figure 24. Case IV -Sphere of Influence of Highway 
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Table 5 

CASE IV - PARAMETER VALUES 



Market Center Locat ons 



Market Center 

1 , 
2 



X Coordinate 

36 
74 



Y Coordinate 

45 
45 



Crop Yields 



Crop A 
Crop B 



30.00 Units 
30.00 Units 



Market/Crop Parameters 



Market 
Center 

1 
I 

2 



Crop 

A 
B 
A 

B 



Market 
Price 

45.00 
40.00 
35.00 
30.00 



Production 
Cost 

5.00 
5.00 
5.00 
5.00 



Transportation 
Cost 

1.00 
0.75 
1.50 
1.00 



Highway Data 



Starting Point 



Ending Point 



Highway X Coordinate Y Coordinate X Coordinate Y Coordinate Alpha 



37 



108 



37 
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Figure 25. Case IV Land Use 
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Figure 26. Case IV-Tributary Areas 
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Case V - Complex Transport Impact 

This case is designed to demonstrate that the application of the simple linear 
model can produce relatively complex spatial patterns in the presence of both 
competing market centers ar\d off-cenier transportation routes. Figure 27 displays 
the spatial situation of market ^enters with respect to the two transport routes, and 
Figure 28 defines the relative spheres of influence of the two transport routes The 
operation of the model produces a very complex pattern of land use, as shown in 
Figure 29, which nevertheless bears a resemblance to that seen in cases II and IV. 
Three Thunen rings form around the first market center while an incipient ring 
structure seems to be forming ba:ed upon the second and third market centers The 
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Figure 27 C iso V Locution of Market Centers and Highways 
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border of cultivation is obviously distorted along each of the transport routes. 
Turning to Figure 30. the tributary areas map, we see that the particular parameter 
structure adopted here has generated a discontinuous market area uu center 2. A 
northeast extension of the market area of center 1 occurs as a result of\he effect of 
the second highway and divides the market area of the second center in^wo parts. 
The market area of the third center displays some of the characteristics s£en in the 
market area of the second center in the previous case. Table 6 gives the parameter 
inputs for replication of this example. 
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Figure 28. Case V Spheres of Influence of Highways 
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Figure 29. Case V-Land Use 
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Figure 30. Case V^Tributary Areas 




46 



\hW 00051 



Table 6 



CASE V - PARAMETER VALUES 







Market Center Locations 




Market Center 


X Coordinate 


Y Coordinate 




1 


46 


50 




2 


55 


50 




3 


64 


50 






Crop Yields 








Crop A 10 00 t Jniu 








Crop B 10 00 Units 








CropC 10 00 Units 








Market/Crop Parameters 




Market 




Market Production 


Transportation 


Center 


Crop 


Price Cost 


Cost 


1 


A 


47.00 19.60 


1.10 


1 


B 


47.00 19.80 


1.03 


1 


C 


47.00 20.00 


1.00 


2 


A 


39.00 19.60 


1.10 


* 2 


B 


39.00 19.80 


1.03 


2 


C 


39.00 20.00 


-1.00 


3 


A 


31.00 19.60 


1.10 


3 


B 


31.00 19.80 


1.03 


3 


C 


3100 20.00 


1.00 






Highway Data 






Starting Point Ending Point 


Highway 


X Coordinate 


Y Coordinate X Coordinate 


Y Coordinate Alpha 


1 


62 


1 50 


90 0.40 


1 


22 


1 94 


90 0.30 
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Case VI - Change of Scale and Origin 

This final case is included to illustrate the manner in which changes of scale and 
origin may be used to illuminate complex spatial patterns. Figure 29 displays a 
relatively complex spatial land use pattern whose interpretation is greatly hindered 
by the size of the printer elements. The insertion of a re-scale card at the end of 
previous example produces a more detailed map of a portion of the onginal arelT 
The origin was moved up and to the right, and the base scale of the map was 
changed from 10 units to the inch to 2.5 units to the inch, producing a fourfold 
enlargement of the subarea. The resultant map shows all three market centers as well 
as portions of the two roads, Figure 31, and the appropriate portions of the sphere 
of influence map, Figure 32. The re-scale operation produces, in Figure 33, a 
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Figure 31. Case VI -Location of Market Centers and Highways 
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significantly enlarged version of a portion of the original land use map. Here the 
relationship between the various crops and market centers can be seen more clearly, 
as can the impact of the two transport routes. Figure 34 shows the market area 
boundaries in the subarea in a more detailed fashion than did Figure 30. Thus, the 
change in scale and location of the origin leads to increased ease of interpretation 
by the student. Table 7 shows the factors which were entered on the re-scale cards 
to produce this case. 
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Figure 32. Case VI Spheres of lnlluerv.c of Highways 
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Figure 33 CiscVI Land Use 
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Figure 34 Case VI Tnbuia*v Areas 



Conclusion 



The examples presented here demonstrate the variety of Ways in which 
LANDUSH can provide a basis lor (he development of computer-based laboratory 
exercises lor eionomtc geography courses, The flexibility of the program permits it 
to be utilized in course* on nearly any level and the type oi exercisesdevc loped can 
be tailored to any group of students. 

LANDUSL is still undergoing development (e g . it will soon be possible to study 
the elicits ul systematic, ureal variations in productivity) and comments and 
suggestions fmin users aie welcomed. * 
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DETAILED OPERATING INSTRUCTIONS 

Version 2.0 / Department of Geography 

Northwestern University 

LANDUSE ^ 

Purpose: Development of land use and market area patterns according to 

the basic Thunen model. * 

Description: Within a rectangular area this program develops the pattern 

of land use according to the basic Thunen model: R = Y(P - C - TD) , 
where R is the economic rent, Y is the yield per unit area, P is 
the unit market price, C is the unit production cost, T is the 
unit transportation cost per unit distance and D is, the distance 
to the market. The use producing the highest rent level is 
assumed to occupy the unit of land. In addition to N market 
centers and M crops, NROAD highways m w y be defined for a given v 
problem. Crop yields (Y) are assumed constant for each crop over 
the total area. » 

In general, the addition of a highway will expand the amount of 
usable land by reducing the effective distance from a cell to a 
market center and, thereby, reducing the transportation cost. 
The restriction uf a, the ratio of on-highway to off-highway 
movement charges, to 0 < a < 1 implies that the effective distance 
; along the highway is reduced by the factor a below the same 
s distance off the highway. As it is possible for a problem to have 
1 multiple highways, it is necessary to consider all possible paths 
j including the direct path from a given cell to each market center 
I and to select from these the minimum effective distance. ^ Let 
J i^, i^* ... i f be subscripts which identify r highways to be 

traveled in sequence from a)cell to a market center, and assume 
/ further that this is the minimum effective distance path. Let 
D^ be the minimum distance from the cell to highway i^, the * 

minimum distance from the market center to highway i r , and 

D ij> •*•» D i ^ e ^ e di stances traveled along the r highways. 

Then the effective distance D is given by the formula: 

D = D_ ♦ a. D ♦ a. D. ♦ ... ♦ a. D. ♦ D~ 

c x i x \ x 2 x 2 x r *r ^ X 

The rectangular area relevant to program computations i^s between 
1 and 110 units horizontally and 1 arid 90 units vertically. This 
corresponds to the printed maps whic.h are 54"lines at 110 
characters per line. There are 10 characters per inch horizon- 
tally but only' 6 lines per inch vertically so that the computer 
output must be scaled by a factor of 3/5 to accomodate the printer 
and provide the same scale vertically and horizontally. All 
coordinates locating towns and roads are, specified as arbitrary 
integer numbers; however, only the towns and portions of highways 
lying-within' the rectangular area defined by 1 < X < U0 and 
1 < Y 1 90 will -appear on .the output maps. 
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The program currently allows maxima of 10 market centers, 12 
crops, .and 9 highways. Since the minimum distance must be 
computed from each cell to each town, problems could be devised 
using these maxima which might require a month or more of 
computer time. The lorfer set of maxima of 4 market- centers , 
4 crops, and 2 highways has been included in the program for 
student usage. Access to the expanded set of limits is pro- 
vided by punching a decoding "password" of up to 10 characters 
in columns 71 - 80 of the problem control^ card. 

Data input to the program is an arbitrary sequence of problems 
each one beginning with a problem control card. The problem 
control card contains, among other things, N, the number of market 
centers, M, the number of crops, NROAD, the number of highways, 
and NEWCASE, the number of re-scale runs to be made, A problem 
is re-scaled by specifying a new origin and a new scale factor 
for the maps. The original (reference) problem is then recomputed 
with updated coordinates for the market centers and highways. 
As N must be greater than zero, a data card giving the market 
center coordinates always follows the problem control card. If 
*M > 0, the scale/crop-yield cards followed by the crop parameter 
cards appear next in sequence. If NROAP > 0, the highway data 
cards (one per highway) are next. Finally, if NEWCASE > 0, 
the re-scaling cards (one per re-scale for a total of NEWCASE) 
are included. It should be noted that program executioli will 
be terminated and an error message will be printed out under the 
/ *' following conditions of the input data:< 

1) N < 0 or N > MAXN, where MAXN is either 4 or 10; 

2) M > MAXM, whew V AXM is eithe- 4 or 12; 

3) NR0AD > MaXR, where MAXR is either 2 or 9; 

4) a highway of zero length is specified, i.e. the initial 
' and final coordinates are identical. 



The program output consists of A] an echo check of data inputs, 
Bj frctt one to four maps (depending upon the proMe^ control 
card parameters), and C] simmary tables giving areas in cultiva- 
tion and levels of economic rent (if M > 0) . The maps are: 
1) a location of market centers map; 2) a sphere of influence 
of roads map; 3) a land use map; and 4) a market areas map. 
The first map always appears in the output, and the second 
appears if NROAD > 0. Maps 3) and 4) are printed if M > 0. 

The market center location map contains an X at the position 

of each town appearing within the previously defined rectangular 

area. Each highway or portion of a highway (if any) within this 



ERIC 



54 ' 

00058 



rectangular region is also displayed with identifying numbers 
at the initial and final coordinates. An option exists for 
selecting the highway symbol. It may he either a dot (•)» 
an asterisk (*), or a plus (♦)'. The default option is a iot. 

The sphere of influence of a highway is the set of cells for 
which the average effective distance from the cell to all N 
market centers is less than the average airline distance to 
the towns or the average effective distance to the towns using 
all other highways. A cell within the sphere of influence 
of a highway is identified on the map with the highway number. 
Those cells whose minimum average distance is the average 
airline distance arc left blank. If the average effective 
distances are equal for two highways with respect to a given 
cell, an equal (=) sign is printed in the cell. 

The land use nap. displays the regions of crop cultivations 
with a letter symbol printed in each cell associated with a 
crop. The market areas map indicates with a number symbol over 
the same region the market where the crops are sold. An equal («) 
sign is printed in those cells where it is equally profitable 
to grow either of two crops or sell the product at either of 
two markets, respectively. Regions for which the rent is <, 0 
are left blank.. 

If MsO, distances on the ,map may be scaled by a multiplying factor, 
SCALE. This factor has meaning only if M>0 as a result, of the 
effect of distance on the rent computations. All points on or 
off the< maps are specified in terms of integer coordinates so 
that a unit increment in either the horizontal or vertical direction 
corresponds to a distance of 1/10" on the printed output. The 
factor SCALE is a non-negative floating point number (if SCALE 
* 0.0, it is reset to 1.0) which assigns a physical distance to 
this unit increment. The physical distance may be either abstractly 
defined or in feet, miles, etc.; and 10.0*SCALE is the distance 
per inch of the printed output maps. This distance scale is 
printed at the bottom of each map if M^O. To obtain maps of dif- 
ferent regions within the same problem requires: the modification 
of the coordinates of all market centers and highways and/or 
SCALE leaving the remaining data the same. In this case physical 
distance must be invariant, i.e. modifying the value of SCALE or 
shifting the origin necessarily causes a consistent modification 
of the integer coordinates which will preserve actual distance. 
Thus, 

sd * c 

where 

c x the physical distance and is constant 
s = the scale factor 
d = distance on the map 
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Let us assume that data has been prepared for a problem 
which will produce output maps in a coordinate system (x,y) 
and that maps of the same problem in a coordinate system (x^y') 
are desired. Also, let s and s* represent the distance scales 
in the uKprimed and primed systems; and let (Xq^q) be the 
unprimed coordinates of the origin of the primed system. If 
(x,,y«) represents an arbitrary point in the unprimed system, 
its corresponding point in the primed system is (xj,yj) where, 

xj « s/s' (Xj- x Q ) 

y\ * s/s» (y r y Q ) 

All coordinates of market centers and highways must be transformed 

as above and rounded off to closest integer values-. Note that 

the new coordinate system is completely specified by x Q , y Q , and s* 

The following example illustrates the salient points of the above 
discussion. Assume that data has been prepared for a'problem 
having a market center located off the maps at ( x i> y i) s (120,-20) 
with a scale of 4.0. For the same problem it is aesired to obtain 
maps of the region surrounding the market center and having a 
scale of 3.0. Let (x 0 ,y 0 ) * (90,-40). Then 

xj « 4/3.(120 - 90) = 40 



* 4/3 (-2Q - (-40)) = 26 2/3 * 27 



and the new data would include the point (40,27) as the location 
of the market center. ^ 

The program contains an option which wi 1 I au* or a t i ca U t re-scale 
any problem. The problem as specified by \he daca immediately 
following the problem control card will bewailed the reference 
case. Each re-scale data card which follows a reference case 
contains the coordinates of a new origin and a new scale factor.* 
The coordinates of the market centers and highways in the reference 
case are- saved and new coordinates are computed relative to the 
reference problem for each re-scaling. Re-scaling may be done if 
M « 0, in which case the reference scale factor is assumed to be 
1.0. The number of the reference case is printed for each re- 
scaling along with the x and y coordinates of the new origin. The 
•cale ratio, i.e. the ratio of the reference scale factor to the 
new scale factor (s/s*) is .also printed. Under'the existing input 
data structure a problem may be *re-scaled as many as' ninety-nine 
times. The user is cautioned not to specify a scale^ factor on a 
re-scale card so large that a highway is effectively reduced to 
zero length, i.e. the initial^ and final coordinates coincide. This 
will cause execution to be terminated. 
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Deck Make-Up: A) Basic structure [CDC 6400] 

1) Systems cards 

2) EOR card r 

3) Program decks 

4) EOR card 

5) Data deck 

6) EOR card 

7) EOI card 

B) Data deck structure - The following sequence is repeated 
once for each problem: 

1) Problem control card 

2) Market center location card(s) 

3) Scaling and crop yield card(s) [if M > 0] 

4) Crop parameter card(s) [if M > 0] 

5) Highway parameter card(s) [if NROAD > 0] 

6) Re-scale card(s) [if NEWCASE > 0] 

Card Format: Problem Control Card - FORMAT (512, 7A10) 

Cols. 1-2 N = number <k centers [1 1 N 1 10 if password 
is used; 1 <\n f. 4 otherwise] 

Cols. 3-4 M = number of \r°j)S^f^*~ 12 if password is used; 
M <^4 otherwisej 

Cols. 5-6 NROAD = number of highways [NROAD 19 if password 
is used; NROAD <^2 otherwise] 

Cols. 7-8 ISYM < 0 if highways are to be disniay<?d on all 

maps; ISYM >_ 0 if highways .are to be displayed on 
the locat-on""of market centers map^only. |ISYM| k 1 
if highway symbol is to b<? an asterisk (*) and 
} ISYM j ■ 2 if the symbol is to be a plus (♦) » The 
highway symbol will be a dot(«) for any otner value 
of ISYM. 

Cols. 9-10 NEWCASE * number of re-scale cards for the current 
problem (if > 0). If NEKCASE <_0, the problem is 
not to be re-scaled. 

Cols. 11-70 Any title information 

Cols. 71 - 80 Any word with up to 10 symbols. If this word is 

the password, the limits on the number of centers, 
* crops and highways are increased. 



M arket Center Location Card(s) - FORMAT (8110) 

The X and Y coordinates are specified pairwise in integer format 

at field width increments of 10 in the following sequence: 
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Xj, Yj, the X and Y coordinates of center 1 
X~, Y 9 , the X and Y coordinates of center 2 



*N* t ^ ie X Y coordinates of center N 



Note that [—-5—-] ♦ 1 cards are required.* 



Scaling and Crop Yield Card(s) - FORMAT (8E10) t 
Cols. - 10 SCALE = the distance scale factor. Ten times 
this number is the scale printed below each 
output map. If SCALE = 0.0^ it will be reset to 1.0. 

Cols. 11 - (Y(I), I * 1, M) = the crop yields per unit 
area for crops 1 through M, respectively, at 
field width increments of 10. 

Note that [M/8] + 1 cards are required.* 



Crop Parameter Card (si - FORMAT (3E10)t 
Cols. 1-10 the unit market price 
Cols. 11 - 20 the unit production cost 

Cols. 21 - 30 the unit transportation cost per unit distance 

Note that there is one data card per crop per market center 
for a total of (N) (M) cards. The crop parameters for crops 
t through M with respect to market center 1 are read first; 
parameters of crops 1 through M with respect to market center 2 - 
are read in next; etc. 



Highway Parameter Card(s) - FORMAT (4110, F10) t 


CoU. 


1 - 


10 


X coordinate of starting point 


Cols. 


11 


- 20 


Y coordinate of starting point 


Cols. 


21 


- 30 


X coordinate of ending point 


Cols . 


31 


- 40 


Y coordinate of ending point 


Cols. 


41 


- 50 


ALPHA, the ratio of on-highway to 








off-highway movement charges. 
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Note that there is one data card per highway for a total 
of NROAD cards. 



Re-scale Card(s) - FORMAT (2I!0, E10)t 

Cols. 1-10 X coordinate of n$w origin in reference system 
Cols. 11-20 Y coordinate of new origin in reference system 
Cols. 21-30 SCALE = the new distance scale factor. Ten times 
this number is the scale printed below each 
output map if M > 0. If SCALE * 0.0, it will be 
reset to 1.0. 

* [X] means the integer part of X. x 

t The E10 FORMAT allcvs the user to specify input data either in 
For E form within a field width of 10 columns. The decimal 
point must be punched in any case. 

» 

Note: The program LANDUSE requires the subroutines 

PROJCT and REGION and needs 50000 g central memory locations 

to compile. The program only requires 33000 o locations to 

o 

load and go. It is written in FORTRAiN IV for the CDC 6400 
computer. Execution tine is on the order of frcn 5 to 20 
seconds for problems having up to 2 high'./ays, but increases 
exponentially with an increase in the mJiber of highways 
beyond 2. A strong factor affecting the amount of running 
time is the complexity, i.e. multiplicity of intersections, 
of the highway network. 



Programmers: Bruce Anderson and Duane F. Marble, August 1970. 
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APPENDIX B 



THE GEOGRAPHY PROGRAM EXCHANGE 



The Geography Program Exchange (GPE) has been established to assist 
universities and other non-profit organizations with the interchange of computer 
software which relates to problems of a geographic nature. The GPE operates on a 
self-sustaining basis with the assistance of a modest grant from the Mathematical 
Social Sciences Board Its services are made available to qualified users anywhere in 
the world. 

The GPE maintains a central file of specialized computer programs with 
associated documentation and test data sets, and makes copies of all or part of this 
material available to qualified users at cost. The programs held by the GPE 
represent specialized materials developed for use in geographic research and 
teaching, r 1 attempt is made to duplicate standard statistical programs. Distribu- 
tion of available programs is in the form of punched cards or BCD tapes, and 
associated documentation is also distributed in machine readable form. The GPE 
also produces a brief guide to its program holdings which is available without 
charge. Current program holdings are based mainly upon the specialized program 
libraries of the University of Iowa, Michigan, Michigan State and Northwestern 
University as well as individual programs from many other institutions. The GPE 
actively solicits programs apd. if the program submitted appears to be of enough 
general interest, it will be added to those distributed. All software handled by the 
GPE is obtained and distributed on a non-royalty basis. 

In addition, the GPE also publishes an irregular newsletter which, at present, is 
incorporated in Current Research Notes in Quantitative and Uieoretical Geography. 
This is distributed without charge by the International Geographical Union 
Commission on Quantitative Methods and requests for copies should be directed to. 

Professor Forrest R. Pitts * 
Social Science Research Institute 
University of Hawaii 
Honolulu, Hawaii 96822 U.S.A. 

Although the GPE operates with limited resources, it does attempt to provide 
assistance to those departments and individuals who are beginning to make use of 
the digital computer as a research and teaching tool. 

Overall management of the GPE is provided by a board consisting of Professors 
Duane F. Marble (Northwestern), Waldo R. Toblcr (Michigan), Robert I. Wittick 
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(Michigan State), and Forrest R. Pitts (Hawaii). Professor Witftck is in chaw ot 
operations and all inquiries relating to activities of the GPE should be directed to 
him at the following address. 

Professor Robert Wittick 

Computer Institute for Social Science Research 

Michigan State University 

East Lansing, Michigan 48823 U.S.A. 
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APPENDIX C 

KEYPUNCH INSTRUCTIONS 



Vogelback Corputing Center 
Nor thwos tern University 
Evans ton, III ino is 



IBM 2^ Printing Card Punch 



General Operation . 

Informa t ion is punched 
which punches letters, digit 
operator presses the "fued" 
from the hopper. When the 1 
Che card is positioned under 
left to right, one column at 
auttar^u ical ly positioned for 
the card automatically is te 
it passed colurrji by colur 
punching station. When the 
at the reading station moves 



in a card by opeiatmg a typewriter ^tylc kevboaid, 
s, or special characters into card column;., Ihe 
button on- the keyboard which roves a card down 
register" button is depressed/ the first column of 
the punching station. The cards ate punched tiotn 
a t iir.e . A f tcr each card is punched , the card 1 ? 
the ner. t column. After all SO colunn> j: t punched, 
glstered ur.dor the reading station, through uhwh 
si In phase with the card following it at the 
card at tho punching station is reload d, the card 
into the card stacker. 



Because passage through the reading station is in phase with pas>.n.o 
through the punching station, holes sensed at the read station can be duplicated 
into the following card parsing under the punch station. The ability to dupli- 
cate information from card to card reduces the amount of r.mual punching rueded. , 
Data common to a group of cards ray be punched into the first card of a group 
and duplicated into the remaining ones. Depressing the "duplicate" ke> vdn 
accomplish this, or it can be done au torr.at ica lly by the keypunch opeiatmg under 
preplanned control (drum or program card). 

Other functions that can be put under automatic control are printing, 
skipping colums, and ejecting and feeding cards. According to the requ 1 1 eir.ont^ 
of the card forrr.it, a pattern can be established so that the keypunch can >kip 
columns, duplicate fields, or release the card when punching is finished but 
before column 80 is reached. 

Instructing the keypunch to skip, duplicate, or release upon rejchinc. 
certain colunn locations is known as programing. The operator program, the' 
keypunch by punching a control card, wrapping the card around a eylindo: oi 
drum, and inserting the drum on a spindle in the machine.. The holes m the 
card control the operation of the machine, causing it to perform specified 
functions autorat ica Ily . 



Operating Features 
Card Hopper 

The card hopper is located on the upper right side of the machine. The card^ 
are placed in the hopper face forward, 9's down, and are fed front card fsr.t. 
A card is fed from the hopper to the card bed automatically or by iloptC"»s ion 
of a card feed key. "The first two cards to bo punched rrust ho fed by key de- 
pression, but all other card** in the hoppor may bo fed aotunu t u a My undo* 
switch control. 
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Punching Station 



Punching is performed at the fust of two stations in the card bed through which 
the cards pass. To start an operation, two c ards arc fed into the caid bed at 
the right of the punching station. As the second card is fed in, the first, card 
is automatically registered for punching -- that is, it is positioned at the 
punching station. 



Read in;; Station 

The reading station, where the cards arc read for duplicating, is located appro<i 
nately the distance of one card to the left of the punching station. Each card 
that has been punched passes through the read station as the next card is being 
punched. The two cards rove in synchronism, column by colurji, and information 
to be duplicated is transferred from the first card to the second. Reading fron 
one card to another can be controlled, field by field, so that cnl; the desired 
information is duplicated. 

This feature el mi nates card handling for the duplication of information from 
cards prepared durin 6 the punching operation. Card handling it nece .sary only 
when an operation requires duplication from prepunched master cards. 



Card Stacker 

The card stacker is located on the upper left side of the machine on a level 
with the hopper. After each card passes the id station, it is fed into the 
sticker automatically or by key depression, mien the cards are removed from 
the stacker, they are in their original sequence. 



backspace Key 

This key is located below the card bed between the rcafl and punch stations. As 
long as it is held depressed, the cards at the punch and read stations arc back- 
spaced continuously until column 1 is reached. At the same tine, the program 
card, which controls jk»pp;r.g and duplicating, is also backspaced. 



Program Unit 

Automatic skipping, automatic duplicating, and shifting from numerical to alpha- 
betic punching and vice ver&a are controlled by the program unit. F.AlTi ot these 
operation^ is designated by a specific code recorded in a program card. (See 
Program Card Codes J The program card is then fastened around a program drum 
and inserted in the machine, where it is read by a sensing mechanism. The drum 
revolves in step with the movement of the cards past the punch and read stations 
so that the program codes control the operations, column by colurn. 
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PrOftram Control Lever 



Operation of the program unit is controlled by the program control lever, 
located below the prOgra\n unit, When this lever is turned on, the program 
sensing mecl anisra is lowered so thait it rests on the program drum, and the 
codes punched in the program card control the various automatic operations. 



Column Indicator 



The indicator, located At the base of the program drum holder, indicate* to 
the operator the next column to be punched. Refer to this indicator us a gu-rde 
for spacing or backspacing to a particular column. 



Keyboard s - Number 3 - Combination alph ab etic and numerical keyboard, phis Li 
special characters . — 

The combination keyboards have the best features of both a typewriter and a 
numerical key punch. Ihe letter keys are arranged for operation by the stan- 
dard typewriter touch system, while the digit keys are placed so that a rapid 
three- finger touch system can be used. The usual numerical ke>b on a typewntei 
have been elimruttd, instead, a group of dual-purpose keys at the r;>Jit m»ivo> 
for digit as well as letter punching. This permits runvncal punch im; with the 
right hand. The touch s>stem for the ten numerical ke> > is. index finju-i t.u 
digits I, 4, and 7,^ middle finger for digits 2, 5, and S. and ring fuutW Uu 
digits 0, 3. >, and 9. The punching of a digit or a Utter with an> ot the 
comhmatior ..eys depends upon the shift of the keyboard. For example, depi le- 
sion of the 4 7 J key punches a A when the keyboard is in numerical shift, hut 
a J when in alphabetic shift. This shifting is similar to upper or towei ^.im' 
shifting on a standard typewriter and may be controlled au toma t ica 1 1\ by tlw 
program unit or manually by key depression. 

The fourth row of keys contains four special character keys at the Uti. llu^<; 
keys punch eight characters as shown on the key tops, four in iiumeru.U ami 
four in alphabetic hift. 



Program Card 

The program card controls the automatic operations for the corresponding tolurm 
of the cards being punched. The control ranching required m the program card 
depends on the functions to be controlled, e.g., skipping, duplicating, and 
alphahct ic punching. 



Field De f in it 'on (12 punch 

A 12 punch (+) should be punched -n every column except the first ^ left -hand 
position) of every field to be skipped, duplicated, or manually puncl ied , The so 
I2's serve to continue to the end of a field any i>kip or duplication started 
within that field. 
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Automatic Skin t \ I punch ^"J. 

An 11 punch in the first colurn of any field automatically start* a skip, 
which is continued over that field by the 1 2 1 s punched in the iei;»iinin£ columns 
of the field. If a single colur.iri is to be auto.ratic.ill> skipped, it is punched 
with an 11. This coding operate* in conjunction with the automatic skip- 
duplicatt switch, which must be on to sta^t the skipping automatically. 

Automat ic tp lication (0) 

A zero punched in the first column of any field au toruit ical ly starts duplica- 
, tion, which is continued over that field by the 1 2 ' s punched in the remaining 
columns o field. If a single colur.ji is to be automatically duplicated, 

it is pun^ ith a zero. 

Alphabeti c Shift (1) 

When the program card is in the machine, the combination keyboard is ncrnullx 
in numerical shift, and depiession of any one of the two-pu*pose keys cai'M> 
a figure to be punched, To punch a Utter, the combinat tor' keyboard must tv 
shifted for alphabetic punching. This shitting is performed au toma t u a I 1\ l«\ 
a I in the program card in each column of the alphabet c f*eld. 



Proy;^ h Card Codes 

Code Func t ion 

12 + Fie Id De f in it ion 

11 - Start Au tora tic Skip 

0 Sta r t Au toma t ic Duplication 

1 Alphabet ic Shift 

2 heft Zero Print 

3 Print Suppress ion 

Pro gram Drum 

Ihe program card is mounted on a program drum for invirlu«n in tin nutliuu. 
Ihe program drum has a clamping strip to hold the caid, and a handle on the lop 
to tighten or release the strip. To fasten a card around the ilium, laki the 
following steps. 

1. Hold the drum in a horizontal position with the handle to the ri^lit . Tut n 
the handle away (courite - -c 1 ockw ise) as far as it will go. This loosens the 
smooth edge of the clamping strip. 

2. Insert the column 80 edge of the card under the smooth edge of the clamping 
strio. Two alignment check holes in the clamping strip make it possible to 
see that the card is flush with the metal edge under the strip. Hie card 
should be positioned so that the 9's edge is against the rim of the drum 
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3. Turn the handle to the center posit tun. This tighten* the strooth edge of 
the clamping t>tnp and loosens the toothed edge. 

j * * 

. Wrap the card tightly arour.d tne drum and insert the column 1 edge under the 
toothed edge of the clanging strip. 

5. Turn the handle toward you (clockwise) as far as it wi'l go. Tins fastens 
the toothed ed£e of the clamping strip. The drum is >-ow ready to be fnserte 
in the machine. 

To remove a card from the drum, the above procedure should be followed m the 
reverse order. 

With the progran sensing mechanism raised, the drum is placed on the spindle 
under the center cover of tie punch, positioned so that the aligning pin falls 
in the aligning hole in the column indicator dial. The program c»> trol lever 
is turned on to lower the reading star wheels onto the program card, and the 
release key is djpressed^to engage the reading mechanism fully. Whencvei the 
drum is to be removed, the program control lever must first be turned off to 
r^ise the star wheels. 

When the punch is operated without progran co nCrol t the keyboard is in al oha - 
betic shift. When the program control lever is .turned on, the keyboard is in 
numerical <hift for all columns except those punched with the alphabetic field 
coding in the program card. If it is necessary to punch digits with the piojnam 
off or in an otherwise alphabetic field, the nurner nal shi f t key is lie Id down 
while the digit keys are depressed. Similarly, when it is necessaiy to puiuh 
some letters in an otherwise numerical field, the alphabetic shift ke y is. held 
down while the letter keys aie depressed. 

\ 

Skip Ke y 

Tie purpose of this key is to skip any field that is coded with I2's in the * 
program card. Depression of the key starts the skip and the 12 * s , in succeeding 
columns of the field, continue the skip for the rest of the field. 

Dash Ski, (or) Dash Key 
/ 

When .the keyboard is"*in numerical shift, depress >f thus key punched ju LI 
hole and causes skipping for* the rest of the fielc determined bv the miu^- 
slve columns punched 12 in the'prorgram card. When the ke>boaid .i b m .i Ifh.iK't u 
shift, depression of this Key punches an 11 hole but does not c.uise tkin«iu»,. 

.TO ^ 

Dupl icate Key * 

i 

Any reformation, manually, punched in one card and cordon to one or mule sutcet.- 
sive cards, can be duplicated by use of this^ key. Depression of this key at 
the beginning of a fiel<f starts duplication, and the 12 * s punched in the program 
card continue the duplication for the rest of the f ^ld. 



6 ' 67 

ERJC , 



00070 




Shift 
key 



C 



Spa 



Shift 
key 



Punch ing Keys 

Keys I- IS car. be depres^d op!> when the keyboard ib in alphabetic .shift, to 
punch the letters, ind n. at* 1 J". If ore of these keys n depressed while the » e> - 
boT<J is in ru .:r^t-al shirt, the ~achir,e lock . Operation can bu rusi^T. J ~> 
releasing the card, or b> ce ui.i^ fie alphabetic shift key which will cusc 
the letter to be punched. 

Combination keys 19-29 can be depressed when the ke> beard is in either purjrica! 
or alphabetic shift to punch the characters controlled by that key. 

Keys purch multiple digits for eight special characters. The> can bj de- 

pressed when the keyboard is in either numerical or alphabetic shift. 



Space 3a r 

The space bar can bo d^vessed at any_ time in a ranual field to cause spa*, inn 
over ere colurr. of the card. 



Functional Keys 

30 NUM (nurerical shift) shifts the combination k > beard into nunencal posi- 
tion as long a* it is held do'*n. 

31 AlPH (alphabet ic" shift) shifts the tombiration kt> board into a Iphaber im- 
position as long as it is held dc*Tu It is normally u;ed to permit the punching 
of letters in an otherwise numerical field. 

32 DUP (duplicate) causes puneh>in£ by reading from tKe preceding card. 

34 RBI. (relca.«0 causes the cards at the punch and r^ad stations ,to be advanced 
complete I > past those stations. Fields pro^ramr-ed for automatic duplication 
beyond the point of release are punched m the card before release is coi ,*h t^. 
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35 



FKrD (caid food) iLMubc a card fc*»! cycle 



1. Feeds a card from the hopper. 

2. Regis tets the card > at the punch ai:c read stations. 

3. Stacks the card frou the left'of th^ card be J. 

6 SKIP cause- skipping of the field for which it is depressed, 

J7 REG (card register) xs used primarily when inserting ca:d» nuncalH. 

44 MULT PCj. (n.ultiple punch) is held down to permit the pumhin£ ot note than 
one digit in a colurn. Hie ke>board is in nu^-ncal shift when tins key is 
depressed. 

(Note* An example *f the use of the m-ltipU punch key is 
| END OF RECORD 

$ * ; * * y y- * v * ★ * y< a * **. * * * * * * * y. * >. * y y. • * v- * * * ft * y y. 

Sap^l*" inf.^ra i „q;d fm Fortra n state- . nls . (Note that "A" punch defines an 
alphabetic lie Id length.) , 



f 



73 SO 



col 1-5 numeric (fitaterint number) 
6 numeric (continuation) 
7-72 alphabetic (to purch nureric mforntion, depress NUM key) 

73-80 autocratic skip and release of card 
★ * * * * * y- * * * * -v y * * ft * * v * * ft y ft y ft v, v- *' y **/-**■.. 



Sample proaran cu d for data ca rds. . 

I 6 11 21 31 61 m> 

col. 1-5 ID (miner J c) 

6-10 Skip 

11-20 Kuner u 

, 21-30 Skip 

31-60 Duplicate preceding card (nu-ienc) 

61-80 S^ip and release card 



* Inf Tritn«:i f.-i th: instruction >he. * ua > ^htai-uJ J ro .« tK H v ' p. 'mivUkm 

IbM JUvlerrru i ".run 1 (Vrd Pun. h_ *j I 1 rj nj |^__ r 2£ ' j* u jlJ' » r«'«'A2«-0"i 0, 

IV5°, .fmi s<.tu Ta* .a, Djua^ P_i oi »' *' i , Nor thw.» , t •■: n l , suut i .il\ Pi* -s , f 5 , 
pp. .8-51 . " 
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with National Science Foundation support. The ideas presented in these papers do 
not necessarily imply endorsement by the AAG Single copies are mailed free of 
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CHAPTER I 

DIGITAL COMPUTERS AND THE TEACHING 
OF COLLEGE GEOGRAPHY 



Although the large digital computer represents a relatively recent innovation in 
American universities, it is already established as a major tool for research 
geographers There is, in geographic research as in other disciplines, a definite 
interaction between tools and problems, Some problem areas may he dormant for a 
period of years or even decades because the necessary research tools are unavailable 4 . 
When a new research tool becomes available, either through primary innovation or 
through the transfer of knowledge from workers in other disciplines, there is 
frequently a surge of interest in topics which had previously been set aside, as well 
as development of new areas of investigation. 

An example of this phenomenon within geography can be seen in the impact of the 
high-speed, digital computer on problems of multivariate regionalizatton. The 
approach to an optimal multivariate classification of areas based upon a comparison 
of between region and within region variances had been developed, but empirical 
implementation was forced to wait until computing machinery capable of handling 
the calculations became available When it did, a strong upsurge of empirical 
interest in this topic followed These empirical investigations raised, in turn, 
additional theoretical questions and the cycle advanced through another phase. 

Parallel to the development of the computer as a research tool in many 
disciplines has been its use as an instructional aid. Larly applications of the digital 
computer in instruction involved the use of the machine to permit students to 
handle problems of a si/e and scale of complexity that had hitherto been impossible 
witfiin the standard classroom situation. In engineering, for example, students 
began to handle problems which, were considerably more realistic than the 
traditional, hand-computed examples This proved quite fruitful, and the intro- 
duction of the computer as a laboratory tool proceeded rapidly in many 
engineering disciplines and, in some schools, its use has been firmly established for 
well over a decade An illuminating discussion of this approach can be found in the 
report of the Ford Foundation project at the College of Engineering of the 
University "of Michigan (I7| Instructional use of the computer has lagged 
considerably in most of she social sciences, though there have been some recent 
attempts witfi in geography and other social science disciplines to provide programs 
of this nature. For example, see the work by Rushton (11 J as well as the more 
elaborate urban simulations such as that of Warneryd [ I9| . 

Recent technologic! developments have 1 permitted the central computer to 
handle a large number of users on a nearly simultaneous, real-time basis Efficient 

9 . , 
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time-sharing systems permit a number of users, equipped with remote terminals, to 
utilize a central computer through the intermediate efforts of a complex supervisory 
system which examines the priorities and needs of each user and allocates to him 
only that portion of the mam computer's capacity required for his work, and for 
only that period of time which is necessary to process the next step m his problem. 
Much of the capacity of these time-sharing systems arises out of the great disparity 
-between the internal computation speeds of the central computer and the speed at 
which information can be transmitted and received from the remote terminals. This 
disparity can be seen in an examination of the relative printing rates of most of the 
commercially available remote terminals. These are very slow when compared with 
the main computer's dedicated peripheral printers and the speed factor is 
frequently ' on the order of 250 or 300.1, if not more Even with this speed 
differential, the mam computer is frequently bound by the capacity of its "high 
speed" printers and will operate two or more of them simultaneously in order to 
attain greater output volumes. 

The introduction of ttmc-sharing facilities in many universities has prompted the 
development ot computer assisted instruction (CAI). There has been a great deal of 
interest m CAI on all educational levels ranging from preschool children to college 
graduates, and many ol the CAI developments in college geography arc discussed in 
previous technical papers of the Commission on College Geography. In the standard 
CAI environment, the time-sharing computer interacts with one or more students 
utilizing remote terminals. A scries of queries arc presented by the computer, based 
upon a program previously prepared by the instructor, with each query in the 
program designed to illuminate a particular point relating to the topic under 
discussion (e.g., map projections, water balance, etc.). After each query the student 
presents a response which the computer checks against a series of stored answers 
which have also been provided by the instructor, and then, depending upon the 
nature of the student's response, branches to another section of the program. 
Internal loops arc common in these programs so that the student can be given an 
opportunity to re-examine material that is giving him difficulty. This "match, 
branch and loop" format is basic to most CAI programs, In many cases, the CAI 
programs are not really computational ones since they utilize only a very small 
portion of the available facilities of the central computer and arc normally 
restrained to "match, branch and loop" operations. Although many CAI programs 
arc run on large, third generation digital computers which provide the necessary 
time-sharing environments the "match, branch and loop" mode of these programs 
does not really require the facilities of a computer of such power and some 
universities have responded by installing smaller computers which handle only CAI 
work. An informal survey conducted in 1970 indicated that only a handful of U.S. 
geography departments were currently utilizing CAI on cither an experimental or 
operational basis, but that many of the departments had a strong irtcrest in future 
uses of CAI. 

It is possible to utih/c the basic "match, branch and loop" format of the CAI 
program in combination with the major computational capacity of the machine. 
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This could be done, for example,, in a CAI unit dealing with problems of spatial 
optimj^aJion which, after providing instruction in the basic elements of linear 
programming and its application to spatial problems, would require the student to 
set up a specific problem in a, linear programming format. The CAI program would 
call upon a resident linear programming routine to solve the student's problem and 
then display the results of this computation* to the student together' with an 
appropriate commentary Because of the multi-processing nature of large modern 
computers, the CAI program could continue the foreground conversation with the 
student while the linear programming problem was running in the background. This 
represents a rather sophisticated use of CAI, but one that will be increasingly 
utilized. There are, however, several basic disadvantages to current CAI approaches 
which arise largely out of the computing environment as it stands at the present 
time. These include: 

1) The lack of comparability of the programs. A great variety of CAI 
programming languages are available, some of which can be utilized only on 
specific machines, and sometimes only at specific installations. Other 
languages which appear to have greater comparability, particularly the BASIC 
language, suffer from a lack of standardization which inhibits the easy 
transfer of CAf programs from one installation to another. 

2) The slow printing speed of myst remote terminals. An inexpensive terminal 
tends to print so slowly that *here is difficulty in retaining student attention, 
while terminals with higher printing speeds are, at present . unduly expensive. 

3) The difficulty of producing any form of graphic output at a remote terminal. 
This is an impediment which appears particularly critical folr geographic 
applications since many geographic CAI applications would benefit greatly 
from the output of simple graphics as well as maps. Available graphics 
terminals are currently both specialized and expensive, while printing a one 
page map at a Teletype terminal utilizing a printing area of, say 50 rows by 
70 columns would require the printing of 3,500 map elementsmharacteis) 
and would take five minutes or more. As terminals with higher printing 
speeds become available, this time will be cut to two minutes or undlh but 
even this is not adequate At present, then, geographically based exercises 
which require more than exceedingly modest amounts of graphic output, or 
which require map output, can be carried out only on expensive remote 
terminal installations or by utilizing the batch processing capability of the 
central computer 

It must bc„ emphasized that these problems will undoubtedly be overcome in the 
next few .years by rapid increases in the level of computer and communications 
technology It would appear that a profitable interim approach would be to develop 
a series of more complex programs based on highly standardized languages, such as 
FORTRAN, and oriented toward batch processing rather than time-sharing. 

As technology improves, these programs could be converted to the "foreground- 
background" CAI mode mentioned earlier. Experiments have shown that the 
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program described in this paper, LANDUSE, can be run effectively in such a 
configuration, but to do so at present requires the use of a high-speed terminal, a 
specialized communications link and a 300 line/minute printer. This type of 
installation has a rental rate significantly m excess of S 1,000/month. Within a 
period of four to six years, it should be possible to incorporate these programs into 
a CAI environment in which such equipment configurations are less expensive and 
much more common than they are today. 

The LANDUSE Program 

The program described in this paper, LANDUSE,, is designed as a laboratory tool 
for use in economic geography courses at both the beginning and intermediate 
levels. It is not a CAI program in the current sense of that term, since it does not 
utilize the basic "match, branch and k op" structure. It is, rather, directly related to 
the third type of instructional program termed "generative" by Huke, Fielding and 
Rummage (6). Programs of this nature permit students to enter parametric 
information into a process model and to note the fashion in which the model 
operates upon the information to produce the desired output. Students are 
expected to interpret the patterns produced with reference to the theoretical 
concepts contained in the process model. The impact of changes in the system can 
also be explored through alterations of the parameter values on successive runs of 
the program. In this manner, the program assists the students in attaining an 
understanding of some interrelationships present in complex models of spatial 
environments. A similar approach to industrial location is embodied in the 
computer model developed by Smith and Lee [13] . 

LANDUSE deals with a simple model of agricultural land use (the basic von 
Thiinen model), which is relatively rich and capable of successive complication The 
program was developed after an examination of an early version of the program 
prepared by Rushton [II]. An entirely new program was prepared* rather than 
modify the Rushton program, in order to provide a number of additional features 
in an efficient fashion. These include 

1) replication of the basic case consisting of a single market center and simple 
concentric rings, 

2) introduction of one or more space distorting transportation routes into the 
system in a completely flexible form, and ^ 

3) an ability to change the scale of analysis in a flexible but simple fashion. 

Chapter II provides a summary of the process theory which underlies the 
operation of the program. The third chapter describes the s,ructure of the program 
and some of the problems encountered in classroom use The final chapter presents 
six selected cases which demonstrate the flexibility and scope of the LANDUSE 
program. 
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CHAPTER II 

AGRICULTURAL LANE) USE: 
A BRIEF THEORETICAL STATEMENT 

The basic theory pertaining to the spatial structure of agricultural production 
originated a century and a half ago in the work of a German economist, Johann 
Heinrich von Thunen. Since then, other researchers have produced moic complete 
and complex statements of this theory but we are concerned here largely with the 
simple case as originally presented in von Thunen's classic book, The Isolated State 
[18). In this work von Thunen attempted to discover those laws which influence 
the prices of agricultural products and the way in which price variations are 
subsequently translated into patterns of agricultural land use. Basically, he was 
concerned with competition between various agricultural products and farming 
systems for the use of specific plots of land. A primary assumption of his model 
was that on each piece of land, the activity which yielded the highest net return 
from utilization of the parcel would be found and that competing activities would 
be displaced to locations where t|iey, in turn, would yield the highest net return. 

To facilitate discussion of thi$ problem, von ThUnen utilized a hypothetical case 
- his "isolated state" wherein he posited the existence of a market center located 
•at the center of a large, fertile plain which contained.no transportation routes. The 
undifferentiated plain in which transfer costs are a simple function of distance, the 
"isotropic plain," appears as a starting point in many geographic theories. Von 
Thunen further assumed that throughout the plain the soil was equal in quality and 
that the agricultural productivity of specific plots of land did not vary from one 
location to another. At some great distance from the town the plain turned into an 
uncultivated wilderness which effectively cut off ail communication between the 
"isolated state" and the outside world. In this isolated community, the central 
town would supply the surrounding rural area with ail manufactured products and, 
in return, would obtain all of its agricultural inputs from enterprises located on the 
plain. The problem as posed by von Thunen was what pattern of cultivation would 
take shape, given the stated conditions, and how would farming systems in different 
portions of the plain be affected by their distance from the central town? An 
explicit assumption in his work was that all agricultural entrepreneurs made 
decisions in an economically rational fashion. That is, entrepreneurs responded only 
to profit maximization as an optimizing goal and were assumed to be in possession 
of complete information pertaining all choices. 

Economic Rent and Location Rent 

As the basis of his analysis, von ThUnen introduced the concept of economic 
rent While he developed this notion independently, it had also been introduced a 
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few years earlier by an Lnghsh economist, David Ricardo, who was interested in the 
way certain factors conditioned agricultural production deci^pns. l J »e wa*s, however, 
concerned primarily with variations in productivity rather than with the impact of 
purely spatial factors. The Ricardian version of economic writ, as contrasted to that 
of von Tinmen, is clearly set forth by Michael Chisholm in his book. Rural 
Settlement and Land Use 1 2 1 . ^ 

The concept of economic rent is* central to manv questions of Isuid use 
competition and it is utilized in the vo'fi Tinmen model to provide a means whereby 
this competition is resolved in a rational fashion The observed gap between 
revenues and costs represents economic rent in a technical sense Von Tinmen 
observed that Ricardo had baseXl his discussion of the nature ot economic rent upon 
differences in the inherent productivity of \anous parcels of agricultural land, but 
that the same phenomena would anse*it productivity were held constant and the 
parcels were subjected to varying transportation charges based upon the cost of 
moving their production over different distances to the central market. In Ricardo's 
case, differential fertility of the plots influenced the level of gross returns, \whereas 
in von Tinmen's model,, gross returns remained constant, but the cost of ( 
^ transportation varied with distance to tiie-market. 

The two cases* are quite similar, buUrecetiHiterature has tended to apply the 
term, "location rents" to those economic rents arising solely out of differences in 
location. This serves to differentiate them from more general scarcity rents which 
may arise frou. factors such as simultaneous variations in productivity as well as 
location. It should be made clear, however, that location rents constitute a specially 
defined subclass of general economic orsjarcity rents and are identical to the latter 
only in the case whereinjocation 1$ the only factor differentiating the parcels of 
land. 

Given the assumption that we are dealing with a finite group of activities (which 
"in themselves may consist of fixed bundles of subactivities, e.^., crop rotation'' 
schemes) it is then possible to view competition for the use of agricultural land as 
equivalent to the following process, each entrepreneur washing to engage in an 
activity evaluates the utility of a given parcel for carrying out a specific activity and i 
determines the level of rent which can be expected. This fi^jrc icpresents the 
maximum amount that the operator wifl be willing to bid for use of the parcel of 
land since any bid higher than this estimate would require that the activity in 
question be operated at a loss. The bids for use t)f the parcel by various 
entrepreneurs desiring to undertake different activities may be ranked and the 
maximum identified The activity associated with this maximum bid will then, 
become the one which occ^)ies that particular unit of land. This notion may be, 
extended to*non-agricultural activities, of course, and is roughly equivalent to the 
"highest andV>est use 1 * notion which is frequently encountered in the urban land 
economics literature. In fact, a significant proportion of the current theory dealing 
with the internal spatial structure of the city rests upon the extension of this notion 
of location rents to form a basis for land use competition in the city, ^.g. f 
Alonso ( I ] . 
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A Mathematical Statement of the von Thiinen Model 

In his original work, von Thlkie^was concerned with specific cases, but he did 
produce a general statement relating fixation rent to distance from the market 
center. Typically the basic von Thiinen model can be stated as a simple linear 
equation- 

R» = Y I I> 1 Y,C, Y,T;D= YjIPj C, t\d| , (i = I,. .,n) (I) 

where. R, = Uv ' t fo. the \ th activity on a spcciPc parcel of land 
Y| s t .tion/umt of land for the \ th activity 

P, - uk. ^uce at the market for-the \ lh activity 
Q^f unit production cost for the i' /; activity 
T, = transport eharge/unit/kilometer for the \ th activity 
D = airline distance of the parcel of land from the market in km. 

This equation states that the bid rent for activity / at a specific parcel is equal to 
the gross returns on production from that parcel (Y,P,), less the total cost of 
production (Y,C,), le^s the cost of moving that production to the market (Y,P,D). 
The simplest situation^ obtained when we assume the existence of an isotiopic 
transport plain in which T, is"a linear function of distance and is assumed to have an 
intercept of zero, i.e.* there are no terminal charges in the system. It is also 
normally assumed that productivity is constant over the plain. That is, Y, does not 
vary locationally within ihc system and that a similar condition hoids with respect 
to unit production costs thus yielding a constant value, C„ for every parcel. The 

' unit price at the market is also assumed to be constant and the market will absorb, 
at that price, any amount which is offered. 

For a specific activity, these assumptions leave us with identical net returns (Y t P, 
Y,C,) excluding transport costs for e\ery parcel on our hypothetical plain. The 
only variations in rent are then introduced by variations in trapsfer costs (Y,T,D) 
which are based cn*ircJy upon differences in the value of D. That is, when the 
parcel is located arfhc market site the, third term in equation (i) is zero and 
location rent is equal to gross returns minus total cost of production As parcels 
further and further from the market are considered,, the si/.c of the third term 

\Y,P,D) in the xjuat&n increases. At some distance from the market transport 
charges will be equal to gross returns and the rent will. fall to zero. This distance 
marks the limit of feasible cultivation for the activity since production beyond this 
point would be unprofitable. The term relating to unit production costs, C,, is 
defined in a somewhat unique fashion. It includes. not only ttfe normal cost of 
carrying out the activity but also includes a specific term representing abnormal" 
profit for the operator This means, then, that location rents can be viewed as a 
potential "excess" profit for the entrepreneur. 

Since equation (1) is linear and continuous, it is easily represented in a graphic 
form. In Figure I, a bid rent curve for a specific activity has been plotted for 
parcels of land located at successively increasing distances from the market center 
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M distance — -> D 

l-igurc I . The Bid Rent Curve for a Single Activity 

at M. The rent R is the net return derived from a parcel located at the market place 
and is equal to YF YC since transport costs at that location are zero. The distance 
MD represents the maximum distance at which a parcel can economically engage in 
the activity, lor distances greater than MD either a subsidy must be provided or it 
must be assumed that the operator is not behaving in an economically rational 
fashion. The slope of the bid rentairve is obviously equal to the unit transfer cost in 
the system and becomes increasingly steep as unit transfer costs increase. Figure 2 
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Figure 2. Bid Rent Curves lor Two Competing Activities 

displays a competitive situation where there are two activities, A and B, whose 
different bid rent curves result from differences in unit market price, unit 
production costs, or unit trans^r costs between the two activities. In each case, 
these values remain constant for the activity at any location, but they are permitted 
„to vary between activities. As in the previous example. K A represents the rent at* 
the market for a parcel utilized by activity A, while K B represents the rent at the 
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market fW a parcel utih/ed by activity 8. The slopes of the two bid rent turves 
differ significantly reflecting higher unit transfer costs for A than for B If there 
were no competition present in thr system and each activity were considered singly , 
then D A would represent the maximum distance at which activity A could be 
undertaken with D b (Dh > D A ) representing the maximum limit of cultivation for 
activity B However, the two activities are in competition for use of the panels and 
within the distance MD* the bid rent of activity A always exceeds that of activity 
B, and hence activity A will then occupy all parcels of land between the market and 
location D*. However at distances greater than D* activity B produces a higher bid 
rent than activity A and activity B will occupy those parcels of land lying at 
distances from the market greater than D* and less than Dr. An entrepreneur 
located on a parcel of land exactly D* kilometers from the market would be 
indifferent between the two activities since both produce identical bid rents at that 
point and no economic basis exists for differentiating between them. 

The graphic analysis of the competitive situation shown in Figure 2 is only 
one-dimensional and total rent in this case is equivalent to the area under the curve 
Ka r Db. Wnen trans,ate( l t0 the two-dimensional case, two rings centered at M are 
formed. The first, with radius MD* encloses the are;: of cultivation of activity A; 
the second, starting at distance D* and extending to radius MDr encloses the area 
of cultivation of crop B (see Figure 3). The resulting spatial pattern is one which is 



commonly referred to as "lluincn Rings" I he total rent in the two-dimensional 
case is then given by the volume of the solid of revolution formed by rotating 
K^RDb about KaM as an axis (see Appendix A of Dunn |3| for a more detailed 
discussion). 

A great deal has been written on the von rhiinen mode' and for discussions of its 
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many complwationsSrfc! extensions, see the works by Dunn [3|, Isard [7j, and 
Alonso [1J. The wtility of even the simple model in explaining the observed 
patterns of land use has been noted b> a number of researchers and much of this 
latter work is admirably reviewed in Chisholm [2] . 

Extensions of the* Simple Model 

The basic von Thunen model outlined above displays a number of important 
factors relating to agricultural production location decisions, and does so in an 
exceedingly simple fashion. While it may be characterized as one of the simplest 
models of land use, it is exceedingly fertile since it is capable of Extension and 
modification in a variety of ways. Two of these, the expansion to'mfcltiple market 
centers and the introduction of space distorting transport routes will be discussed 
here. 

Multiple Market Centers 

Given the concept of land use competition, the extension of the basic von 
Thilnen model to the case of multiple urban centers is relatively straightforward. 
Basically. equation ( I) is recast into the following form: 

R,j = Y,P,j Y,C, - Y,T,jD, (i = l n; j = 1, . . .,m) (2) 

where R,j = rent produced by the \ th activity when it is oriented toward the j" 1 
"market center 
P,j = unit price for activity \ at the market center 
T,j = transport rate/umt/km for movement of output of activity i to 
market center j 

(the remaining terms are as previously defined) 

Utilization ol this formulation results in the establishment of multiple bid rent 
curves based upon the various market tenters (see Figure 4) The structures are as 
defined in the previous example except that, in this case, an evaluation of the bid 
rent curves shows that at certain locations, for example, L, three "different" 
activities are bidding for the use ol the land. At L we see positive bid rents derived 
from activity A orieiu. J toward market 1 and activity B oriented toward market 1, 
as well as a bid rent derived from activity B oriented toward market center 2. At L, 
activity A oriented toward market center I produces the highest bid rent and hence 
occupies ill at location. However, at distances greater than Mi? we see that 
orientations ol the two activities, A and B, toward the second market center 
produce higher bid rent levels than a continued orientation toward the original 
market center. Thus, at distance M,L* from the original market center, it can be 
seen that while activity B oriented toward market center 1 is still making a positive 
bid lor ii.* 1 ol the land, this bid is exceeded by those made by activity A and activity 
B when oriented toward the second market center. In this particular case, activity 
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Figure 4 Bid Rent Curves for Competition Between Two Market Centers 



A. with an orientation toward the second market center, occupies the location. 
Thus, along the line from M, toM 2 , we see a succession of land uses, first A then B 
and then A again. However, within the area devoted to the cultivation of land use B 
there is a breaking point, P, at which the land shifts from the tributary area of the 
first market center to the tributary area of the second center. 

An interesting problem involves the exact determination of the boundary 
between the tributary ateas of the two centers. The boundary is clearly defined in 
the one-dimensional case v shown in Figure 4, and can be easily inferred in the 
two-dimensional use where the parameter values for the two market centers are 
identical and an overlap between the bid rent curves exists, the boundary in this 
latter case will be a .>lraight line perpendicular to the line joining the two centers. 
However, as one of the centers becomes "dominant/' the boundary between the 
tributary areas of the two centers will gradually shift from a straight line to a 
hyperbola to a parabola and finally to a closed elliptical form around the second 
area with the second (smaller) market center occupying the foci of the ellipse 
nearest to the larger market center Demonstration of the development of the 
spatial structure in fhese cases is usually very difficult without resort to 
complicated mathematical derivations or an elaborate, arefully plotted graphical 
display. It is. nevertheless, an interesting extension of the basic model. 



System Wide Transport Improvements 

The introduction of system wide transport improvements, sudi as those that 
might occur when a new and more efficient type of vehicle is introduced into all 
parts ol the transportation system tan be handled within the model as given in 
equations (I) and (2). The net effect is simply to reduce the value of T, or T,,. 
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which, in turn, leads to a reduction in the slope of the bid rent curves. The benefits 
of any system-wide improvement in transportation ma> be differentially distributed 
between activities since, for example, the introduction of a specialized and more 
efficient bulk carrier might favor certain commodities over others. Another 
example of differential impact would be the introduction of refrigerated, milk 
transports which would obviously influence T, for the activity of milk production, 
but would have little direct impact upon the transport costs for commodities such 
as grain. 

! 

Route-Specific Transport Improvements 

A far more interesting, and theoretically complex, case is encountered when 
transport improvements in the system differentially influence locations. This is 
typically caused in the real world by the construction or improvement of an 
existing transportation link. For example, 3 new highway linking two nodes in the 
system may be built or an existing road may be improved through resurfacing, thus 
leading to lower vehicle operating costs. Within the context of the basic von 
Thdnen model, the introduction of specific transport routes can be viewed as a 
distortion of the uniform plain since unit transfer costs over the improved routes 
will generally be lower than transfer costs over the remainder of the area. Within 
the context of equation (1), this change can be handled in one of two ways: 
through variations in the values of the T, depending upon the type of route over 
which the travel takes place or, somewhat more easily, by assuming that variations 
occur in D itself. This latter approach arises out of a substitution of the notion of 
effective distance for that of airline distance which was previously utilized in the 
model. Effective distance is defined as follows, let D c be the minimum distance 
from the cell to the highway, D t the minimum distance from the market center to 
the highway, and D, the distance travelled along the highway. Then the effective 
distance D is given by the formula: 



with a in this case lepresenting the ratio of on-highway to off-highway movement 
charges within the range 0 < a. u 1. That is, for an a value of 0 5 travelling a unit 
distance along the highway would cost exactly half as much as travelling the same 
distance off the highway. The use of a as a scaling factor in the equation for 
effective distance compensates for the real differences in unit movement costs over 
the different routes. This proves to be particularly useful when the definition of 
effective distance is expanded to include travel over a potential system of routes 



D C = D c + a D, + D t 



(3) 



Di, = D c + a, D, + <x 2 D 2 + . . - + a n D n + D t 



(4) 



where there arc a series of n different highway segments in the route, each of which 
has associated with it an unique value of a. 
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The distorting effect of low cost transportation routes was noted by von ThUnen 
in his original analysis in a discussion of the distorting effects of a river winch 
would provide low cost transportation for agricultural commodities upon the limits 
of die area under cultivation. Derivation of ihe geometric impact ot highway 
improvement upon land u*e patterns has been undertaken by several authors, but 
the discussion that follows is adapted from that of Werner (20), Werner was 
conceme ' with the distortion of the sh^pe of the hinterland area that occurred 
after a new transportation route had been constructed. He defined hinterland as the 
area accessible from a given point with a specific level of transportation effort, 
where effort may be measured in terms of energy consumption, mileage, time or 
other criteria. The analysis assumes as gi\"n the isotropic transport plain, a point P 
corresponding to a market center and a straight transportation route running 
through the area with a cost per unit flow per kilometer of a, on the road and a 2 
off the road. This basic situation is illustrated in Figure 5. 

In Figure 5, for convenience,' the new transportation route is assumed to be 
identical with the x-axis and P is a point on the y-axis with coordinates (0,b). The 
situation as shown is a symmetric one so that the analysis need only be completed 
for points lying to the right of the y-axis. The minimum coSt connection from an 
arbitrary point Q 2 located sufficiently far from P and lying off the transportation 
route is given by a direct movement from Q 2 towards the road entering with an 
angle equal to arc cos a 2 /a,, then along the road itself to a point R, and then 
directly to point P. The angle at which the minimum cost path leaves the road is 
again equal to arc cos a 2 /a I . 

Werner demonstrates in ins analysis that the half-space under consideration can 
be divided into three regions, which are labeled I, II, and HI in Figure 5. Within 
region I all points have a minimum effective cost path which consists of a direct 




I'igure 5 Determination of the Sphere of lnlluetu,e of J Transport Route 
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movement to point P with no movement over the road. For ail points in region II, 
such as the arbitrary point Q 2 , al! minimum tost connections involve travel over the 
road in the fashion outlined above. The situation in region ill is more complicated 
since for points arbitrarily close to the road, minimum effective cost paths will 
involve a combination of direct travel and travel over the road. For points much 
further from the road, direct travel will be more efficient. The boundary separating 
the two parts of region III, is shown to be a parabola with the point P lying in one 
of its foci. The resulting situation is summarized in Figure 6 which displays the 
areas of access via the new road versus the areas which retain direct access to point 
P. The shaded portion in Figure 6 represents those points in the hinterland of P that 
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Figure 6. The Sphere of Influence of a Transport Route 
Source. Werner [20] 



have had their distance relationship changed by the introduction of the new 
transport facility , while points lying within the area of direct access retain the same 
distance relationships as before. For those points in the shaded area, the effective 
distance has been reduced by the introduction of the highway and this means that 
they are. in a purely economic sense, closer to P than they were prior to the 
introduction o* the highway By this transformation the areas of cultivation will be 
extended along the road, and in the areas immediately adjacent to it. in the fashion 
shown in f Figure 7 which displays a simple Thiinen ring structure for both the 
before and after cases Werner goes on to show how the patterns of changed access' 
are influenced when a variety ol preexisting transport networks arc assumed and 
the interested reader is referred to his original exposition tor further details (20). 
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nc-N hinterland boundary 



Figure 7 Ring Distortion Produced by a Transport Route 
Source: Werner (20L 



Summary 

Material presented in this chapter has provided a summary of basic aspects of the 
von'Thiinen model and some of its more direct extensions. The reader interested in 
a more detailed exploration of these problems is referred to the works by Alonso 
[ 1 J , Chisifolrh [2], Dunn [3], [sard [7] andLosch[9]. 



15 

00022 



CHAPTER III 



THE LANDUSE PROGRAM 



The objectives in using computer programs such as the one discussed here are set 
forth in the article by Huke, Fielding and Rushton [6] and in the report of the 
Engineering Education Project of the University of Michigan [ 17). The selection of 
this specific process model was based upon several factors: 

1) the von ThUnen model is perhaps the simplest generative model for spatial 
patterns of land use, but is rich tnough to permit a number of interesting 
extensions, 

2) the basic model is widely known and is discussed in most economic 
geography courses, and 

3) although the model is simple in structure, a number of operational 
difficulties frequently arise among students in visualizing ^the changes in 
spatial structure resulting from parameter shifts in the model. 




The program was developed to provide a means whereby the student could 
develop a stronger feeling for the' relationship between the abstract model and the 
spatial structure which it represents through manipulation of selected portions of 
the model and observation of the subsequent impact upon the spatial structure of 
land uses. The conventional "blackboard" approach utilizing one-dimensional 
cross-sections (such as those in Figures I, 2 and 4) has not been proven to be an 
effective solution to this problem. Secondary development objectives included a 
desire to acquaint the student with the complex spatial structures that can arise out 
of even simple linear models (e.g., case V in the following chapter), and to provide a 
large group of undergraduates with an introduction to the use of the high-speed, 
digital computer as an analysis tool. 

Program Structure 

The program presented here has been designed with flexibility of use in mind. It 
has the capability of handling exercises ranging from those that might be 
appropriate for an introductory course emphasizing simple concepts to those 
suitable for advanced courses dealing with more complex problems. Examples of 
the former include generation of simple Thunen rings and the demonstration of the 
development of a competitive boundary between two market tenters. The latter 
could involve any of the many complex 'cases of transport improvement and 
development, in all cases, the model relies upon the basic von Thunen model for 
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evaluation of competition among land uses. Competing market centers and crops 
are readily handled within this Iramework, and one or more distance distorting 
transport routes can be introduced through the concept of effective distance 
discussed in the previous chapter*. 

The program assumes the existence of an unbounded, isotropic plain upon which 
an arbitrary (x.y) coordinate system has been established. It then isolates a 
rectangular subrcgion from the plain and proceeds to develop the pattern of land 
use according to the basic rent equation. Unless otherwise instructed, the program 
assumes that the rectangular subrcgion is one whose lower left hand corner 
corresponds to the /.cro points of the coordinate system. For pragmatic reasons, 
related to both case of computation and case of output, the program assumes that 
this rectangular area is divided into a scries of small areal units which fonn the basis 
for all computations. There arc 110 of these units along the horizontal axis and 90 
of the units along the vertical or y-axis, producing a total of slightly under 10,000 
unit areas All coordinates which locate features on the plain, such as market 
centers, are specified as integer numbers in a read nght-and-up coordinate system. 
These integer numbers correspond to cell numbers. For example, the coordinates of 
a market center located approximately in the center of the rectangular region 
would be given as x = 55 and y = 45. Within the unbounded plain, the student Is 
permitted to locate up to a maximum ol ten market centers through the simple 
process of specifying their x and y coordinates. However, only those towns lying 
within the rectangular area defined by I < x < 1 1 0 and I < y < 90 will appear on 
the output maps. 

The student also specifics a >erics of parameters which define up to twelve 
abstract crops, with each crop specified m terms of the following parameters. 

1 ) Yield per unit area (assumed to remain constant for that crop over the entire 
plain); 

2) Unit price at caoli of the market centers, 

3) Unit cost of production (assumed also to be invariant over the plain). 

4) Unit cost of transport per unit distance (permitted to vary between market 
centers and between crops). 

With this information in hand, the program proceeds to evaluate the rent 
equation for each unit area in turn beginning with the unit in^ic upper left hand 
corner. The distance from this unit to each of the defined market centers is 
calculated and the renU according to the basic rent equation, that would be derived 
at this parcel from each crop and market center combination is obtained. For 
example, in a system that \u\d two market centers and two crops defined, each unit 
area would have four possible bid rents ba.sed upon crop A at market center 1, 
crop B at market center I , crop A at market center 2 and crop B at maiket center 
2. From this vector of tour bid rent values, the program selects the maximum, and 
thus identifies the particular land use which is assigned to the parcel and determines 
the assignment of the parcel to a specific tributary area. If all bid rents at that 
location are negative, then the use assigned is "idle" and the parcel is not assigned 
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to the tributary area of any urban tender. For positive bid rents, the maximum 
value is also stored so that subsequent summaries of total rent* and total areas 
, under cultivation can be provided by crop and by crop-market center combinations. 
Transport routes are defined by the simple expedient of providing the x and y 
coordinates of their points of origin and termination together with a value for a, 
the ratio of on-highway to off-highway movement charges. Program operation 
proceeds as described before except that the minimum effective distance is 
established for the computational cell prior to evaluation of the rent equation. The 
information on minimum effective distance, when compared with direct access 
distances, permits definition of the sphere of influence of the transportation route. 

Subroutines 

The program is written entirely in FORTRAN IV and consists of a main routine 
and two specialised subroutines The mam routine. LANDUSE, is responsible for all 
input and computational phases. The only exception occurs when highways are 
present and the subroutine PROJCT is called to determine the shortest distance 
between the analysis cell urider consideration and a specific hjghway. It will be 
recalled that Werner [20] demonstrated that entry into the transport route should 
take place at an angle which is a function of the ratio of the two costs. PROJCT, 
for reasons of computational simplicity, assumes that entry always takes p! ce at 90 
degrees This introduces some distortion, but it is felt that this can be safely 
disregarded in this particular case. 

The subroutine REGION is responsible for all map output. It is designed to store 
and output up tc three one-page maps using the smallest possible amount of core 
memory. In order to fall into the most favorable level of the priority system, the 
program is written to load and go using only 14,000 words of core. Sto .ge of the 
three 54 x 110 maps would normally require nearly 18,000 words of computer 
memory, but REGION reduces this storage requirement to under 1,700 words 
through the use of masking statements. The amount of storage required b.y 
REGION ts dependent upon the computer word length and since word length for 
the CDC 6000 series machines is significantly greater than for comparable 
computers, such as the IBM 360, the CDC 1 version of REGION will not cfpfejaje on t/ ^ 
other machines. Versions of REGION have been prepared, however, for use on the TV 
IBM 360 as well as on CDC 3000 series machines, both of which have different" 
word lengths than the CDC 6000 series. These modified routines produce output 
maps of slightly different sizes than those described here (the 1 1 0 units in the 
x-dimension of the analysis area represents eleven CDC words each holding ten 
characters), but the general structure is identical. / 

Production of maps on the high-speed printer requires that all output be scaled 
by a factor of 3/5 to produce identical vertical and horizontal scales. This results 
from the si'e of the characters on the print chain (10 per inch) and the spacing 
between lines (6 per inch) This is why LANDUSE uses a 90 x J 10 analysis region 
and produces 54 \ 1 10 maps II this scaling were not done, the Thiinen rings would 
be significantly distorted. 
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Program Limits 

The program allows maxima of ten market centers, twelve abstracfcrops and 
nine transport routes. The number of computations required in the maximum case 
can be quite large^ For example^ fn a problem with ten market centers and twelve 
crops, the program must evaluate 120 rent valu' . for each of nearly 10,000 cells. 
The problem of determining patterns of land use in a system involving a relatively 
large number of rrnrket centers and crops is not nearly as severe, however, as the 
case where there are a large number of transportation routes with multiple 
intersections. Since each transportation route is permitted to have an individual 
value of a the determination of the minimum effective distance becomes quite 
complex. It is possible, therefore, within the defined maxima of the prograr*. to 
establish problems which cannot be economically solved even with a large 
high-speed computer. These rather high limits have been established to provide the 
program with the maximum, amount of flexibility and not with the view that all 
will be reached at the same time. From a pragmatic standpoint, it is desirable to 
have a lower set of limits for student # use. Within the program, these lower limits are 
established at four marRet-centers, four crops and two highways and this is the 
configuration which is utilized in the normal classroom situation. Access to the 
larger capabilities of the program may be obtained by punching a single decoding 
"password" on one of the problem control cards. If this password is absent, then 
the lower set of maximas ate in operation; if present, access is permitted to the full 
capabilities of the program. 

Data Input 

* Data input to the program takes place via a sequence of control cards, the first 
of which is a general problem control card specifying the number of market centers, 
the number of crops, the number of highways and the number of re-scale runs that 
are to be made. This card also carries basic title information which is reproduced to 
identify the output, as well as the password if one is entered. The re-scale option is 
discussed in detail later, but basically a problem may be re-scaled through the 
specification of a new origin and a new scale factor for the map. If a problem is to 
be re-scaled, an indication of intent must be given on the problem control card. 
Following the problem control card, ene or more data cards giving, the market 
center coordinates are read and are followed by a third group of cards containing 
information on the desired scale of the map (given in terms of the number of 
distance units per inch) as well as yield levels fi : each of the crops. The fourth 
group of cards is called crop parameter cards, and specifies for each crop and market 
combination, the unit market price, the unit production cost, and the unit 
transportation cost per unit of distance. If transportation routes are specified in the 
problem, then a filth .group of cards follows specilying the starting and ending 
point of each transportation route as well as its associated a value II the problem is 
to be re-scaled one orjnore times, then the final group of cards consists of re-scale 
cards specifying x and y coordinates of the new qngm as well as the new distance 
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scale factor. ^Any number of problems can be run in sequence by repeating the 
above deck structure- 
Detailed instructions on the preparation of input cards are given in Appendix A. 

" * 

Program Output 

The output of the program consists of a series of maps as well as basic summary 
statistics. The first output stage consists of an echo check of data inputs and 
repeats, in tabular form,, the information pertaining to market centers, crops and 
^"highways that has been read from the data cards (see Figures 8 and 9). This permits 
the student to check quickly for errors in his control cards. Following the echo 
t l heck, from one to four maps are produced depending upon the problem control 
card parameters. The firsLmap is always a general location map wh'ch displays the 
location* of all market centers and highways in the problem. (see Figure 10). All 
maps are identical in format with market center locations shown as "X M \ and all 
carry small crosses at the corners of the rectangular analysis area to serve as 
registration points should the student desire to overlay one or more maps for 
comparative purposes. If highways have been defined in the system, then the 
secon J map is one showing the location of the market centers as well as those cells 
which fall within the "sphere of influence" of each of the highways (see Figure 1 1). 

Highways ^re identified by the program by reference to the sequence in which 
the highway control, cards are read. That is, the first card is identified by the 
program, as being associated with highway number I, the second card with highway 
number 2, and so forth. A cell whose effective distance to any of the market 

^. centers is reduced through the existence of highway I has a "1" pnnted in it. 
Similarly, any cell which is influenced by highway 2 has a "2" printed in it, etc. 
Those cells through which direct access still represents the minimum effective 

t distance are left blank. Thus, the area of the map which is filled with "1° *s 

represents the sphere of influence of highway 1, the area filled with "2" *s 

represents the sphere of influence of highway 2. and so forth. 

v 
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Figure 1 1. Output Sample - Sphere of Influence of Highway 




The next map produced is a general land use map (see Figure 12). This map 
shows the location of each of the market centers as well as \ic pattern of land use. 
The first control card read is associated with crop A, the second with B, etc. Any 
cell for which crop A has the highest bid rent has an A printed in it on this map, 
any cell with crop B having the highest bid rent has a B printed in it, etc. Cells for 
which ail bid rents are negative are left blank. For those cells for which two or more 
bid rents are identical an sign is printed indicating that a tie exists between 
two or more land uses. The last map identifies the tributary areas of each market 
center (see Figure 13). The first set of coordinates read is assumed to define market 
center !, the second set of coordinates market center 2, etc. A cell is assigned to 
center 1 if the highest bid rent for that cell is one which is associated with sales at 
market center 1, it is assigned to center, 2 if the sales are associated with market 
center 2. e|c. Ties are again identified by printing of an sign in the cell. 

The final output section consists of a set of summary statistics relating to rent 
and area under cultivation for each market center and crop combination (see Figure 
14). Total rent and total area under cultivation are also given for each crop and the 
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total area with zero rent level, that is, idle land, is aiso presented. The summary 
statistics are useful in providing quantitative measures of impacts due to parametric 
changes. 

Introduction of Transr station Routes 

Transportation routes are defined within the study area through specification of 
the coordinates of their points uf origin and termination as well as the associated a 
value. The a value represents the ratio of on-highway to off-highway movement 
charges and an a value of one implies that these are identical. Thus, an a value of 
0.5 implies that travel on the transportation route is onjy one-half as costly as travel 
off the highwa>. Based upon empincar investigations of vehicle operating costs, it is 
possible to select a values which .are at least roughly in accordance with reality. A 
number of studies have examined the effect of transport improvements on vehicle 
operating costs and if the assumption is made that an a value of 1 is roughly 
equivalent to the crudest ol dirt roads, then relative values of a may be established 
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l igurc 13 Output Sample Tributary Areas 



for other common road surfaces as well as lor different types of vehicles. Table I 
presents some suggested values for a which are based upon empirical investigations. 
These values are only suggestive since many other factors influence the level of 
operating costs. However the fact that they arc based upon "real-workT data tends 
to make them more viable in the mind of the student. 

The sphere of influence of a highway is defined as the set of all cells for which 
thj average effective distance from the cell to all N market centers is less than the 
average airline distance to the towns, or the average effective distance to the towns 
using all other highways A cell winch falls within the sphere of influence of a given 
5 highway is identified on the map with that highway's number, while those cells 
whose minimum average distance is the average airline distance are left blank. If the 
average effective distances are equal lor any two highways with' respect to a given 
cell, an ,<= " sign is then printed in the cell. 

A transport route may be of any length and may run in any direction, but care 
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Figure 14. Output Sample Summary Statistics 



must be taken that, rather through initial specification or subsequent re-scaling, the 
effective length of the highway is not reduced to zero. If this should occur then an 
error message is issued and the program terminates. 

Changes in Origin and Scale 

All points on the isotropic plain are assumed to be specified in terms of integer 
coordinates so that a unit increment in either the horizontal or vertical direction 
corresponds to a distance of one-tenth of an inch on the printed map. The factor 
SCALE is a non-negative, floating point number which assigns a physical distance to 
this unit increment. The physical distance may be either abstractly defined or in 
feet, miles, kilometers, etc., and ten times the value of SCALE is the distance per 
unit of the printed output maps. The distance scale which is in effect for a 
particular problem is printed at the bottom of each map. In order to obtain maps 
covering different regions within the plain, the coordinates of the defined market 
centers and highways must be modified and/or the value of SCALE must be 
changed leaving the remaining data the same. In this case, physical distances must 
be regarded as invariant, that is modifications of the value ofSCALF or shifting of 
the origin must necessarily cause a consistent modification of the originally 
specified integer coordinates which wril preserve the "actual" distances in a spatial 
system. Appendix A contains information on how these scale shifts and shifts of 
the origin are computed by the program. 

The program contains an option which will automatically re-scale any specified 
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Table I 



OPERATING COSTS AND ALPHA VALUES FOR ROAD SURFACE 
AND EQUIPMENT CLASSES IN A DEVELOPING NATION* 



Equipment 
Class 


Road Surface 


Earth 


Gravel 


Paved 


Small passenger car 


3.26 




I .70 




(1.0) 


\V.fO) 


\v.o I ) 


Large passenger car 


5.32 


■4.13 


3.23 


• 


(1.0) 


(0.78) 


(0.61) 


1-ton truck 


6.79 


4.78 


3.56 




(1.0) 


(0.70) 


(0.52) 


3.5-ton truck 


11.79 


7.47 


5.13 




(1.0) 


(0.63) 


(0.44) 


15-ton tractor 


20.53 


12.19 


7.92 


semi-trailer 










(1.0) 


(0.59) 


(0.38) 



♦Figutes arc cents/ vehicle-kilometer, bracketed numbers are alpha values assuming earth 
surface to be equivalent to a = 1.0 



Source Adapted from figures given in United Nations. Transport Modes and Technologies for 
Development, 1970 

problem, A problem which is specified by a data set immediately following a 
problem control card is called the reference case. Each re-scale data card which 
follows a reference case may contain the coordinates of a new origin as well as a 
new scale factor. The coordinates of a new origin on the re-scale card are assumed 
to be given in terms of the original coordinate system. That is, the system may be 
viewed by the student as taking the original rectangular analysis area (which the 
program assumes to be located with its lower left hand corner at die zero point of 
the coordinate system) and moving its lower left hand corner to the new 
coordinates specified on the re-scale card. If a scale change is also specified, the 
relative size of the region and the analysis areas will be changed as well. The original 
coordinates of the market centers and highways, as specified in the reference case, 
are saved and new coordinates are computed relative to this reference problem for 
each subsequent re-scaling. Changes of origin may or may not be accompanied by 



subsequent changes in scale. In many complex problems, it has been found that 
changes in origin together with simultaneous changes in bcale produce useful results 
in terms of clarifying the analysis. For example, a complex case may be developed 
and the student may desire to examine in greater detail some limited poition of it. 
He may then request that the scale of the map be shifted to provide moie detail and 
that the origin be translated to such a location that the area of greatest interest is 
centered in the new map. 

In the case of re-scale operations, the number of analysis cells defined in the 
rectangular region remains the same while their area is effectively changed In the 
present version of the program, this means that the summary statistics in the 
reference case and the re-scaled case will not be comparable. This problem will be 
eliminated : n a future version of the program. 

Using LANDUSE 

The LANDUSE program was develooed at the Department of Geography, 
Northwestern University, over a three year period. The program was initially used in 
an upper division course in economic geography having an enrollment of about 25 
btuoents, most of whom were graduate students from other disciplines. Problems 
were assigned on an individual basis and were fairly complex in their structure since 
the student group possessed a substantial background in felated areas such as 
transportation engineering. The use of the program with this small group of mature 
students posed a few problems. 

During the next year, with the cooperation of Dr. Ian Hamilton (Department of 
Geography, London School of Economics) who was a visiting professor at that 
time, the use of the program was extended to a larger, lower division course in 
introductory economic geography. This course had a registration of .about 70 
students and because of the larger class si/.e, the students were "teamed" in groups 
of two or three in order to reduce operating costs and the load on rather limited 
keypunch facilities. Teai-.ing in groups greater than three did not prove successful 
because of the increasing remoteness of the individual student from the exercise. 

At the present rime, the program is used in both the lower division and upper 
division economic geography courses. The types of problems assigned differ 
significantly with the more complex problems in the upper division course 
reflecting the more substantial background of these students, and their specific 
interest in questions of transportation and its relationship to the spatial structure of 
land use . 

General experience with the program indicates that \i has a high level of student 
acceptance and that it does increase their level of knowledge m the desired areas. 



Necessary Student Background' 

Aside trom a theoretical grounding in the notion of economic rents and location 
rents, and in the general structure of the simple von ThUnen model, the student 
must know how to prepare the data cards for the program, how to submit the 
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program to the computer and how to retrieve his output. In most cases, particularly 
when dealing with lower division students, the instructor cannot assume student 
familiarity with either computers or the operation of the keypunch. Our experience 
has shown that basic keypunch instruction can be given in about 15 or 20 minutes 
and that an introduction to the program, card layout and computer center 
operation take^ about an hour. In many universities, basic keypunch instruction is 
provided by the computer center through either direct demonstration or by student 
access to video recordings. A sample of the keypunch reference material given to 
North western'students can be seen in Appendix C. 

The first exercise winch is normaUy assigned at both class levels is a "set piece" 
wherein all parameter values are given This is designed to test the student's 
familiarity with card preparation and system operation. Case I in the following 
cfiapter is usually used for this purpose. 




Operating Costs 

The cost of operating the program is a function of the charge level established by , 
the local computer center, as well as of the number and complexity of the problems 
which are assigned by the mst/lictor. In the introductory course at Northwestern, 
where only a few simple problems are assigned to teams of students, operating costs 
are approximately $ 14 00 to S 15.00 per student. More complex problems, although 
not a greater number, are assigned in the upper division course on an individual 
basis and here the cost per student rises to between $50.00 and $70.00. In many 
cases these charges, particularly those for the lower division class, reflect a level 
which is related only to the number of job submissions made, since at Northwestern 
a minimum charge of $1.00 per job is set by the computing center. In dealing with 
simple problems, one of the major factors influencing total class cost is the number 
of submissions per student. This can be minimued by making use of the provision 
within the program for stacking problems so that a number of different cases can be 
run on the same submission When used, this tends to reduce significantly the cost 
per student. , 

The cost per student is> also increased when students substitute repetitive tests on 
the computer lor thinking about the operation of the model. Prior to being given 
die exercises, the students are instructed to produce trial solutions utilizing 
standard graphic methods, such as those shown in Figures 2 and 4. However, the 
student finds that it is less work to change merely one or two numbers in his input 
deck and resubmit the problem. This substantially increases the number of runs and 
the total cost. Some computer installation^ have monitor routines available that will 
check the number of submissions made by an individual student and will either 
prohibit access after a certain number ot tries or will call the situation to the 
attention of the instructor Such a program is not available at Northwestern, but a 
comprehensive monitoring prugiani is being developed which should alleviate this 
problem and permit the instructor to identify tiiose students who are having 
difficulties with the operation of the program. 
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While the dollar figures cited here refer to a $7.00/minutc computer charge, it 
must be kept in mind that at many schools the assignment of dollar values to 
classroom use of the computer does not result in any actual transfer of funds. The 
"charge" in these cases provides only a means ot allocating a fixed resource (wh.Ji 
is paid for in the school's annual budget) among competing users. Charge policies 
vary widely among schools nd the proscpctive user should make local inquiries. 

Program Availability 

Because LANDUSE is a large, complex routine which is undergoing continuing 
development, a listing of the source deck is not included in this technical paper. 
However, machine-readable versions of the current distribution version may be 
obtained from the Geography Program Exchange. The activities of the Exchange 
and methods of obtaining programs arc described in Appendix B. The GPE will 
provrde either CDC or IBM versions of the program and persons requesting 
machine-readable copies should specify the machine which will be used. 

Persons using the program are requested !o forward their comments to the 
author. 
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CHAPTER IV 



SOME APPLICATIONS OF THE LANDUSE PROGRAM 



The following six cases have been chosen to demonstrate the versatility and 
flexibility of the program. Parameter values were selected to reflect this goal. 
Exercises featuring more "realistic" parameters can be developed easily within the 
context of a particular course. The cases presented here begin with a simple 
situation which reproduces the basic ThUnen ring structure and then proceed 
through sucessive complications to a complex case of transport impact combined 
with a demonstration of the use of the re-scale option to provide a finer grain 
examination of the results. 

For case of interpretation, boundaries and market center locations have been 
emphasized through the addition of hand inserted symbols in most illustrations in 
this chapter. 



Case I - Simple Thunen Rings 

This case demonstrates the ability of the model to replicate the formation of 
simple Thunen rings around a single market center, it is normally utilized asa 4t set 
piece" in order to familiarize the student with card preparation, program 
submission and interpretation of the output. The parameter values for generation of 
this case are shown in Table 2, and Figures 15 and 16 display the resulting land use 
pattern and market area structure. 
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Table 2 

CASE I PARAMETER VALUES 



Market Center Locations 



Market Center 
- 



X Coordinate 
5? 

Crop Yields 



Crop A 
Crop B 



30.00 Units 
30,00 Units 



Y Coordinat 




Market/Crop Parameters 



Market 
Center 

1 
1 



Crop 

A 
B 



Market 
Price 

35,00 
30.00 



Production 
Cost 

5,00 
5.00 



Transportation 
Cost 

1.50 
LOO 



BIBBBBBBB8B 
BBBBBIlBBBBBBBBBBBBBB 
BBBIBMIIIBBBBBBBBBBBBBBttBe 
BBBBBBBBBSIBBBBHBBBBBMBRlBBBiB 
BBBBSHBBBBBBBBBBMftBBBBBBBBlMBBBBBB 
BBBBMBilBBBWlBBBBBBBBBB*6IB**BBB86BB 
BBB8BCTMBBBBBBBBBBBBM8BBB8BMBBBBBBBBBB 
BBBMBBBBBBBBtBBBIIBBBBBBBBMBtMBBHBBBBBBB 
■BBBBBMBBBBBBBBBBB BBBBBIHB BBBBBBBMBBBBBBBBB 
BBBBMBBBBftMBBWlf^ntn%|MBBBMMBBBftBBBBH 
•BBBMIBMMBB»B#Tu****A«lXkeSBe..9BBBBBB5BB 
BBeBMBBBBB#BMBjir*«««*****AA**aB6BBIBBIBBBIMSB 
HWWIUIIMlflAMUU kkkMkkk AAMBBBitBBIBBOBBlB 
BBBBBBBBBBB8BBIP AAA A A AAA A A AAA AAAAAbHBBBBBBBBBBIBB 
tlBBBBBBBBBBBBl* AAAAAA AAAAAAA AAAABBBBBBMBBBBBBB 
BBBBBBBBBBBBBBBtAAAAAAAAflAAAAAAAAAfcBBBBBBBBBBMBB 
BB8BBB8lB8B8B8^> A A AAAA A A AAAAAAAAAAaBBBBBBBBBBBBRW 
BBBBBBMBBBBBB^AAAAAA* AAAAAA AAAMBBBBBBBBBIBBBB 
|BBBBiBiBBBBBBB%AA»»»*A»»"»»A*#BBBBllMBBB»lBB 
BMtBMMBBMBBtm#>********A#fBBBBBBBBBB»B»0«B 
BB«miBBB«mBBBt«4M44"#MBMBBBlBBB6BBBBft 

bbbbbisbbbbbbbbbbbimiUIIIbbbbbbbbbbbbbiibbb 

BBBBBBf 8BB8BBBM8BBBBBB88BBB8BBBBBBBIBBBBBB 
BMBBMMfMBBBBBBBBBBMBBBBBlBBBBBBBBBBB 
BBBBBBMMBBBBBBBS6BBBBBBB6BBMBBBB8BBB 
BBMBBBBB*BBBBBB8B68BBBB6B88BBBBBBBBB 
BBBBBBBBBBBBIBBBBMBBBBBBBBBBBBlB 
8BB8BBB»BBBBBBB8vJB6B8BBMRBB 
••BBBMIBBBBIIBBBBBBBBB 
*BSf«BBBBM8BB a 

Figure l^. Case 1-Land Use 
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Figure 16. Case I Tributary Area 
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Case II - Competition Between Market Centos^ 

The case of competition between two market centers can be dealt with easily 
within the context ol the program. Case H represents a somewhat arbitrary 
selection of parameters leading to a situation where the first tenter is ''dominant" 
with respect to the second. The land use pattern which results is displayed in Figure 
17. An examination of this map shows that Thiinen rings develop around both 
centers and that land use B merges between them. A more interesting pattern is 
revealed when the tributary areas map, Figure 18, is examined. Here the dominance 
of the first center is clearly evident and the boundary between the trade areas of 
the two centers bends in the direction of the second center which occupies an 
off-center location within its trade area. Through suitable adjustment of the 
parameters, the boundary between the two trade areas may be bent even further to 
the right and m the Extreme, caused to close on the far right thus embedding the 
trade area of center 2 entirely within the trade area of center 1. Table 3 gives the 
parameter values which will reproduce this case. 

Table 3 

r 

CASH II PARAMETER VALUES 



Market Center Locations 



Market Center X Coordinate Y Coordinate 

1 36 45 

2 74- 45 



Crop Yields 

Cmp A 30.00 Units 

Crop B 30 00 Units 



Market/Crop Parameters 

Market Market Production Transportation 

Center Crop Price Cost Cost 



I A t 45.00 5 00 1.00 

1 B 40 00 5.00 0.75 

2 A 35.00 5.00 I 50 
2 B 30.00 5 00 1 00 
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8BBBB8BBBBBBBB8B»BB$«BBBBBBB8 
BBBBBBBBBBBBBBBBBtfBBBBBBBBBBBBBBBSBBB 
BBBB8B8B8BBBBBB8BBBBBB8BBBBBBB.BBBBBBB8BB8BBSB 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBBBBBBBBbBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBSBBBBBBBBBBBB 
BBBBBBB«8BBBBBBBBB8BBBBBBB8BBBBBBBBBBBBBBBBB8BBBBBBBBBBBB8B 
8BBBBBBBBB988BB9b8BBBBB8BBBBBB8BB8BBBBBB8BBBBB88BBBBBBBBBBBBBBB 
8BBBBBBB8BBBBBBBBBBBBBBBBBBBBBBBBB8BBBBBB8BBBBBBBBB9BBBBBBBBBB8BBBB 
BBBBBBBBBBBBBBBBiSBiBBBBBBSBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBBBBBBBBBBBBBBBBMBBBBBBBBBBBBBBBBBBBBBBBB8BBBBBBB8B8BBBBBB8BBBB8^BBB8BB 
BBB6BBBBBB8B8BBBBMBBBBBBBBB8BBBBBBBBBBBBBBBBBRB8B88B8BBBBBBBBBBBBBBBBHBBB 
BBBBBB9BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB8BBBBBBBBBBBBBBBBBBB8B8BBBBB 
BBBBBBBBBBBBBBBBBHBBMBBBBBBBBBBBBBBBBBBBBBBBBBBBBB1I8BBBBBBBBBBBBBBBBBBBBRBBBBBB 

bbbbbubbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbBbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 
bbbb8bbbbbbbbbbbbbbbbbbb9bbbbbbjwm&bb88^bbbbbbbbbbmbbbbbbbbb8flb8bbbbbbb8bbbbbbbbbbbbb 

B»BBMBBlBBB98MB88tfBBM8»BBB 8jBBt UUUAU %BBW BBBBB8BJBBBBBBBBaB8BBBBBBBBB8BBBBBBtfB8BBB8B 
P BBBB^9BBBB6BB6BBHBBBBBBJUffnnAAAAAAAAAAAHT%^ 
BBBBBBBBBeBBBBMBBBBBBBjRAAAAAAAA»AAAAAAAAAAlA^B8BBlWBBBBBBB8BBBBBBB8BPBB8BBBBBBB 
RBBBBB9BBBBBBB8BBBB8BJP*****AAAAAAAAAAAAAAAAAAA^bBBBBflBBBBBBBB8BBBBBBBBBBBBBBBBBBBBBB§BBBBBB8 
flBB6HBBBB^BBBBBBBBBBjKAAAAAAAAAAAAAAAAAAAAAAAAAAlfcBBBBBB8BBBBBBBBBBBBBBBBflBB8BBBBBB8BBBBBBBBB 
BBBBBBBBBBBBBBBBBbBjKaAAAA#*AAAAAAAAAAAAAA AAAAAAAAaSMBB8BB8B8B8BBBBBBBMBJBW8 

RBBBBBBBBBBBBHBBBJTAAAMAAAAAAAAAAAAAAAAAAAAAAAAAAAAMBBBBBBBBBBBBBBJBiAAAAA^B^BBBBBBBBBBBBBBBB 
BBB8iBa8BB8BBBBBBfAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^BBBBBBBBBB#lAAAAAAAAAA^BBBB^^ 
RBBBBBHBBRBBBBBI ^AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAa! JbbBBBBBBJ^AAAAAAAAAAAA^BMBBMB^ 
RBB8BMBB8BBBBBBI kAAAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAi iBBBBBBBffA***AAAAAAAAAAAAA«BBBBBBBBBBBBBB 
RRBBBBBBBBBBBBif UA AAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAAI IBRBBBBrAaAAAAAAAAAAAAAAAA^BBBBBBBBBBBBBB 
BBBBBBBBBBBBBBBf IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* IBBBBBBB*AAAAAAAA * ""+ B »» B ^""«^ 
flBBBBBBBBBBBBBBt lAAAAAAAAAAAAAAAAAAftAAAAAAAAAAAAAAAAAA* IBBBBBBSbVaAAAAAAA^AAAAAAA JBBBBBBBJBBBBBB 
BBBB88BBBBB8BBBI UAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAI »WBBBRBB>AAAAAAAAAAA« 

RBBBriBBBBBBBBBBI UAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* IBBBBBBBfjMAAAAAAAAAAAAAAAAJBBBBBBBBBBBBBBB 
RBBBBBBBBBBBBBBf kAAAAAAAAAAAAA AAAAAA AAAAAA AAAAAAAAAAAAJ iBBBBBRRfl^ftAA************AjB"BBBBBBBBBBBBB5 
BBBBBBBBBBBBBBBBV******************^ 

RB88B8BBBBBBBB8BWAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^BBBBBBBRRBflBBBl|||iai||#BBBBbBBBBBB8BBBBBB 
BBBBBBBBBBBBBBSBBBilMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAlBBBBBBB^BBRBBBBBBBBBBBM 

«BBBBBBBBBBBBBBBBBB^AAAAAAAAAAAAAAAAAAAAAAAA*WAJKbBBBBBBBBRRBBBBBBBBR88BBBB8BBBBBBBBBBB8BBB 
BRBBBB88BBBBB8BM8BB» 

r8bb8bbb--bbbbbbbbbb5Kaaaaaaa^ 

RAHHHBBBBBBBBHBB8tfBBBBeeB%AAAAAAA6AAAAAAAAAJP6BBBBB8BBBBBRB8RRBBBBBBBBBR6BBdBBBBBBBBBBB8BBB 

HBdBBBBBBBWBBBBBBBBbBBBBBlftBiAJiAA^ 

HB8BBBBBBBB62aBB8bBBBRB8BB8BBlBTOlTn(CTTO 

RRBBBBBBBBBBBBBBBBBBBHBeBBBBBBB8BBBBBBBBBBBBBRBBBB88BBBBBHBBBRBBBBBBBBBBBBBBBBBBBBBBB 

RBBBHBBBBBBBBBBBBBBBBAB6BBBBBBBRBBBBBBBBBBBB8BB8BBBB8BBRBBB9BBBBBBRBBBBBBBBBBBBBB 

BBBBBBBBBBB8BBriBBBBBbRBB8BBBBBBBB8BBBBBBBBBBBRBBBBBBBBBBBBBBBBBBBBBBB8R8BB 

BRBBBBBBBBBBBtJB^ 

RBBBBBBBBBBBBHBBBBBBBBBBBBBBBBBBBBB8BBBBBBRBBRBBBBBBBBBBB8BBBRRRBBBBBBB 
BBBBBBBBBBBBBBBBBBBRBBBBBBBBBRBBB8BbBBB8BBBRBBBB8BBBBBBRHRBBRRRBRBB 
BBBB8BBBBBBBBBBBBRBBBBBBBRBRBBBBBBBBBBBBBRBRBBBRBBBBBBBBBHBRHBR 
BBBBB8B88BBBB88«8888BBBBRBBBBBBBBBBBBBBBRB8BB8BB8BB8BRBBRRRBB 
BBBBBBHBBBbBB«B8BBBBBBBBBBBB8BBBB8BB8BBB8BBBB8BH8BMBBB8BB 
BHBBBBBBBBBBHBB?J8BBBBBBBBOBHBBR8BBBRBBB8B8BBBB8H«B8 
BaBBBBB8HBBBu8BB9BBB86BB8BBBBBBBBB9BBBBBBRBBH 
8»BBBHdB8dBBB8BBB<lBBBBtfB8BBB8RBHB9BBBBf» 
tt888BHBB8RBBBBBBBBBBB8BBB8BBB8B 



Figure 17. C*:sc II Lind Use 
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Figure 18 Cjs- i\ Tributary Areas 
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Case III - Transport Route Bypassing a Marke* Center 

This case is a replication of case I with the introduction of an* east-west 
transportation route running slightly tu the south of the market center as shown in 
Figure 19 Figure 20 gives the sphere of influence of the transportation route with 
those cells containing a "1" having a minimum effective distance which involves ase 
of the transportation route The thin vertical line of I's immediately below the 
market center results from an anomalous condition wherein crossing of the 
transport route is counted as use of the route. This will be removed in a future 
version of the program. Figure 21 demonstrates the resultant patterns of land use 
and should be compared directly with Figure 15. Although there is an observed 
extension of land use B along the transportation route, there is no impact upon the 
pattern of cultivation of land use A, such as one might have expected from an 
examination of Figuie 7. The reason for the lack of distortion of the A rings can be 
seen through a comparison of Figures 20 and 2 1 . Superimposition of Figure 20 over 
Figure 21 reveals that the sphere of influence of the highway does not extend into 
the A ring and hence no change is recorded for those cells. Figure 22 displays the 
resulting distortion of the trade area along the transport route and Table 4 gives the 
parameter values which will reproduce this case. 



Figure 19 Case III Market Center and Highway Locatio. 
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I igure 20 Case HI Sphere of Influence ot Highway 
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BBMBBBBB6B 
tHIHUHl BB8BBB 6BB6 
IMMMBIINIIIUHIMIUM 
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HIMMHNMHIIIMWItUMIWMlllKM 

bbbbbbmbbbbbbbbbbbbbbbbbbbbbbbbbbbbsbbbb 

m»mih|iimhimhnm»bh«mi8nimhhh 

6M»«BbBMMBBB6BjrMAAAAAAATX%6ftBBBBBaBBB6BS88 
* BABBBWemmBaWAAAAAAAAAAAA *%B6»MB06BBB6BB8B 

888BB8BBBBBBBBl»AAA»»»»AA»»»»»A»»»B8B0a8BB8B80Be 
BBBBBMBBBBMBBjUAAAAAAAAAAAAAAAAJIBtfBBtfBBBBBBBtfBB 
bS8SBIfiMBlM»eiBaBB****AA******A*A**AN9BMBBBBBBBBBBBB 
•BBBBBBBBBBMMBMBMBBBAAAjUAAAafUAAAAAAJBBBBBBBftdBBBBBBBBBBBB 
b»««»»»^6BB0t#B¥bBb0B8*B|^»A»AAAAAAAAAAAAAAk0RB0B8BBe8«Bee8BB6BBeB8l 
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Figure 21. Case III -Land Use 
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Figure 22. Case IH Tnbutary Area 
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Table 4 

CASE III PARAMETER VALUES 



Market Center 
I 



Market 

C' 

1 
1 



Crop 

A 
B 



Market Center Locations 



X Coordinate 



55 

Crop Yields 



Crop A 
Crop B 



30 00 Units 



30.00 Units 
Market/Crop Paramete. 



Market 
Price 

35.00 
30.00 



PFoducuon 
/ Cost 

5.00 
5 00 



Highway Data 



Y Coordinate 
45 



Transportation 
Cost 

1 50 
1.00 



Starting Point 



Ending Point 



Highway X Coordinate Y Coordinate X Coordinate Y Coordinate Alpha 
i 1 39 108 39 0,45 
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Case IV - Market Center Competition and Transport Impact 

This case is a replication of case II with the addition of a transport route 
bypassing the two competing centers a short distance to the south. The spatial 
positioning of the two market centers and the transport route is shown in Figure 
23. Because of the existence of twecenters, a more complex pattern emerges when 
the sphere of influence of the road is examined (see Figure 24). The impact of u * 
transport route produces an extension of cultivation of both crops A and B to a 
point where the new band of - Jtivation of crop A along the highway nearly join* 
the two centers (see Figure 25). There is, however, still a small island of 
production remaining at the edge of the A-B ring around the second center. / .n, 
many of the more interesting results of the insertion of the transportation ro , can 
be seen in the examination of Figure 26 which displa* the configura** of the 
two market areas. Some very significant changes can be observed h jmparing 
Figure 26 with Figure 18. The introduction of the transport route h<u> extended the 
trade area of center I far to the east of center 2 and the lower boundary between 
the two trade areas has now b$en wrapped up and around center 2, so that its trade 
area has become a parabola oriented toward tne northeast. The parameter inputs 
which replicate this case are given in Table 5. 



Figure 23. Case lV-Location of Market Centers and Highway 
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Figure 24, Case IV Sphere of Influence of Highway 
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Table 5 



CASE IV PARAMETER VALUES 



Market Center Locations 



Market Center X Coordinate Y Coordinate 

1 36 45 

2 74 45 



Crop Yields 

Crop A 30 00 Units 

Crop B 30.00 Units 



Market/Crop Parameters 

Market Market Production Transportation 

Center Crop Price Cost Cost 



1 A 45 00 5 00 1 .00 

1 B 40.00 5.00 0.75 

2 A 35.00 5.00 I 50 
2 B 30 00 5.00 L00 



Highway Data 
Starting Point Ending Point 



Highway X Coordinate Y Coordinate X Coordinate Y Coordinate Alpha 
1 I 37 108 37 0.30 
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Figure 25 Case IV Land Use 



ERIC 



42 

00047 



mumuuiiuumimumiimit 
nmmn mum mm mini mm mum 
mmumimmm mtmmmmmimmtm 
iiiiiiiiiiiiiiiiiiiiiiiimiimiiiipiiiii i mm um 
mmmmimmmmmmmimmmmmmim 
mmmtnmmmmmimmmmmmmmmim 
.L , M!J l MI} mmu,mmmmimimuim »»»»»»»iiiiim 
uiiuiiiiuiiiiiuiiiiiniiiuiiuiiiiiiiiiiiiiiiiiniiiiiiii 
immmmmmmmmmmmmimmmmmmmi 
immimmmmmummmimmmmmmmimm 
iiiimiuuiuuiuiiiuiiiiuiiiiuiiiiiiuiiiiiiiiiiiiiiii 
mmmmmmmumimmnmmimmmmmmm 
nnmummimmnmimmimmmmnmiimmm 
miiiiuimiuiuiiiuuiuuiniuimimimmiumiuu 
mini u um i mm mini imiuu minimum mm nm 

lUllllllllUllllllllillllllllllllllllllllliilliiiiUllIlllllil 

uiuimuiuuuuuiiuiiiiiuiimiuuiuimimiumm 
mniimmmmmmmmuiimmmmimimm nm 

IlllltlllliUlllllllltllllllllllllllllliliilillllllllllllllllll 

mimmimmmmimmmmmmmmmmmmm 
uuiiumMuumumuiuimuiiu minimi in iimnu 
numuimumuiiuiummiuuimuiuimmun nm 
mmmmmmiimmiimmimimm iiimmiimmi 
miimmmmmm mm 1 1 1 nmi limn mmmi n unii 
minimum minimi in n iiimmmi im immmm 
uuuuiuumiiuiiiu muu mini minium mm mm 
tmmimmmimiimmmmimmmiiim mmmm 
immmmimim immm mm i imm iiimmmi mi i 
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiiiiiii 
mumm mmimmmmmnm mim mmmim mil 
tmumuuimmmimimmiii imumummiii mm 
imimuiumiimmmiiMmiuimiiiimmmiimim 
mmmimmiimmmimmim imuiiiiiiimiiimiii 
imiimmmmmiimiuiumm mmiiium mil mm 
immmuimuiuimiummimmiiii mini iimiiui i 
in milium mmmm nm niiiii niiimimmimuiui 
immmummmim nmumuim^u ii nm nun mm 
mini mi mini mini mumum i inn mini mm mm 
mm muu mu i mm mum u ii mum itimiiiiiiium 
mumm mumm mum mum mm muu nuimmi 
miumuumnmumnui nun imuuiunn mumm 
mmumm im mum mm uu milium um uuimiu 
mmumummiimiiiiumumi Hum iiumiimiim 
ii im muu i m mm in in 1 mum mm mum iiniium 
mmmmuim in im mu muu i uiim mi 1 1 1 in umiii 
mmumimummiuiutiimiiimuiiuimi nunum 
immmuuiiii mini mim iiim nun u i im mm uu n 
mmumumuimmuiiiniiuu/muuiiiiuum nun 
nmnmmmum mm mmi/mmmuumiii i nm 
mu m mini mum im/iunu mm mini in 
1 1 u i m m 1 1 u m m m/i i n i m u u m m in 
m ii mi muu mini mu mini i in 
• mmuimmnuimiimui 




______ mumm 

1 1 1 1 ilimillm 1 1 m n i ui 
i mmmi iimmmmu 
im mu m imuiuu uu 
immimum mu i u m 
uu imm um mini in i 
mmuim luiuinu i m 
n i uu mi milium mu 
i minium ii n mini m 
iimiiuu iiiuiium m i 
u ii mum uu nun mu 

i in i um mmmm ii m 

ii m muu i im mm u 1 1 
in m 1 1 mi i n i nun mu 
i in i minim i ii um uu 
in m iiuiiiuiii i in 1 1 

n 1 1 mm mu i u 
in 1 1 mm 1 1 
i u i m 
1 1 



imi 
ifcftt 
uu i 



Figure 26 Case IV Tributary Areas 
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Case V - Complex Transport Impact 

This case is designed to demonstrate that the application of the simple linear 
model can produce relatively complex spatial patterns m the presence of both 
competing market centers and off center transportation routes. Figure 27 displays 
the spatial situation of market ^enters with respect to the two transport routes, and 
Figure 28 defines the relative spheres of influence of the two transport routes. The 
operation of the model produces a very complex pattern of land use, as shown in 
Figure 29, which nevertheless bears a resemblance to that seen in cases 11 and IV. 
Three Thunen rings form around the first market center while an incipient ring 
structure seems to be forming based upon the second and third market centers. The 
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border of cultivation is obviously distorted along each of the transport routes. 
Turning to Figure 30. the tributary areas map. we see that the particular parameter 
structure adopted here has generated a discontinuous market area rbr center 2. A 
northeast extension of the market area of center 1 occurs as a result of the effect of 
the second highway and divides the market area of the second center inuwo parts. 
Hie market area of the third center displays some of the characteristics seen in the 
maiket area of the second center in the previous case. Table 6 gives the parameter 
inputs for replication of this example. 
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Figure 29 Case V Land Use 
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Figure 30. Case V Tributary Areas 
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Tabic 6 

CASH V PARAMETER VALUES 



Market Center Locations 



Market Center 

i ~ 

2 

3 



X Coordinate 

^~ 46 

55 
64 

Crop Yields 



Crop A 
Crop B 
CropC 



10.00 Units 
10.00 Units 
10.00 Units 



Y Coordinate 

50 " 

50 

50 



Market/Crop Parameters 



Market 


Crop 


Market 


Production Transportation 


Center 


Price 


Cost 


Cost 


1 


A 


47.00 


1960 


1.10 


I 


B 


47.00 


19.80 


1 03 


! 


C 


47.00 


20.00 


1.00 


2 


A 


39.00 


19.60 


1.10 


2 


B 


39.00 


19.30 


i.03 


2 


C 


39 00 


20.00 


.1.00 


3 


A 


31.00 


19 60 


1.10 


3 


B 


31 00 


19.80 


1.03 


3 


C 


3 1 .00 


20.00 


1.00 






Highway 


Data 






Starting Point 


Ending Point 




Highway 


X Coordinate 


Y Coordinate 


X Coordinate Y Coordinate Alpha 


1 


62 


I 


50 90 


0.40 






1 


<)4 90 


0.30 



ERIC 
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Case VI ~ Change of Scale and Origin 

This final case is included to illustrate the manner m which changes of scale and 
origin may be used to illuminate complex spatial patterns. Figure 29 displays a 
relatively complex spatial land use pattern whose interpretation is greatly hindered 
by the size of the printer elements. The insertion of a re-scale card at the end ot the 
prev.ous example produces a more detailed map of a portion of the original area. 
The origin was moved up and to the right, and the base scale of the map was 
changed from 10 units to the inch to 2.5 units to the inch, producing a iouriold 
enlargement of the subarea. The resultant map shows all three market centers as well 
as portions of the two roads, Figure 31, and the appropnate portions oi the sphere 
of influence map. Figure 32. The re-scale operation produces, in Figure 33, a 



ligute 31 Case VI Location of Maikct Centers and Highways 
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significantly enlarged version ol a portion ol the onginal land use map. Here the 
relationship between the various trops and market centers can be seen more clearly, 
as can the impact of the two transport routes. Figure 34 shows the market area 
boundaries in the subarea u> a more detailed fashion than did Figure 30 Thus, the 
change in scale and location of the origin leads to increased ease of interpretation 
by the student. Table 7 shows the factors which were entered on the re-scale cards 
to produce this case. 
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CCCCCvUCCtUtCCCcCccCCCctCcCC 
CCCCCCCCCCCCCCCCCCUCC CCCCCCCC 
CCCCCCCCCCCCCCCCCcCCC CCCCCCCCCC 
CCCtCCCCCClCCCCCCcCCCCCCCCCCCCCC 
CCCcCCCCCCCCCCCCCcCCcCCCC CCCCCCC 

ccccccccccccccccccccccccccccccccc 



CCCcCCCCCCCCUCCCcCCCCCCCCCCC CCCCC 

ccccccc:ccc:ccccccccccccccccctccc< 
ccccctccccccccccccccccc;cccc:ccccc 

CCCC UCCCCCCCCcCCc CCC CCCCCC CCCCC 
cCcctCCUCCCUCtCc CCCCCCCC CCCCCCCCCC' 



tccuct.tcccccrtZ^ccccccccccccccccrccccctcuctcccccccctaccccccccctc 
^cccc.cccccccccccccccccccccccccccc:ccccccccccccccccrcccccccccccccccc 

^CCCCCCCCCC CCCCCCC vCCCCUCcCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCC CCCCCCCCC -CCCC CCCCCCC CCCCCCC CC CCCCCCCCCC CCCCC CCCCCC CCCCCCC 

iccccccccccccctcccccccccccccccccccccccccc«.cccccccccccrctcccccccccci 
to»eB»8B8M:ccccccccccccccccccccccccccccccccccccccccccccccccccccccctcct tcctt c 




CCCCCCCCCCCCCCCCCCCCCCCCC CCC CCCCCCCCCCCCCCCCCCCCCCCCC CCCCC CCC CCC 
CCCCCCCCLCCCCCCcCCCCCCCLCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCC 

xccccccccccccccccccccccccccccccccccccccccccccctccccccctccccccccccccc 

CCCCCCCCCCCCCCCCCCCCCCCCC CCCCCC CCCCCCCCCCCCCCCCCCCCcCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCLCI,CCC r CCCCCCCCCtCCCCCCCCtC(.CrCCCCCCCCCCCCCCCCCCCCCCCCCC 



CCCCCCCCCCCCCCCCCvCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCcCCCCCCCCCCCCCCCCC«.t.Cl.CCCCCCCCl.CCCl.CCCCtCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCtCCCCCCtCCtCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

ctcccccccccccccctcccccccccccccccctttccccccccccccccccccccctcccccccccccc(.cccccccccccccccccccccccccccccci 

KgiBdttdBlltC CC L L '-CCCCLSCCCCCCCCCCCCCCCCCCLCCCCCCCt.CCCCCCCCCCCCCCCCCCCCCCCCC.CCCCt.CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
HBMHHwwMMUBBBt ^f CC C C LC C C C C C C CCC C C r 1C CL C C C LC c cc c cC ccc c Ccccc c c c c c c c c cc cc c c c c c c tc c c c c c c c c c c c c c c r c C C C C C C C>* uC cc c c c 
»«»BBMBB»BBBl»(>B8e^tCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCtCCCCCCCCCCcCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCtCC 
48«BBBbB8»bBB»BBiiiiBB\CCCCCCC.CCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtC r CCCCtCCCCC£CCCCCCCCC''CCCCCCCCCCtCCCCCCCCC' 
MBbiHBH#b8»»»B8BBesJ*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
BBB8BBHB8bbbbBBBBB8BBB*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC£CCCCCCCCtCCCCCCL-CCCCCLCCtCCCCCCCCCCCCCCCCCCCCCCC(.C 
BiBBB>MBBBBBBBB»BB<t»3BreEcCCCCCC r CCCCCCCCCCC CCCCCC CCC cffMlfcc CCCC CvCCCCC CCCCCCCC CCCCCC CCCCCCC CCCCCCCC CCCCC C 
BBBBBBBBSBBKBS^bllBBB»«bb«aCCCCCCCUCCCC£££CCCCC^ 

^•BBBBBBHBBOBBBaBBBBlfcCCC CCCCCC C UUA ^C C CC CCg BHBBB 8 BBBBBlfc C C C CCC C C C C C C C C C C C C C C c cc c c c c cc c c c c c cc c c C C C C C C C 
tftB«B8BBBBBBBBBBbB&»aCCCCCCCCCCC 1»aaa1*TObBB8BBB«»bBBBBR«CCCCC C CCCCCCCC CCCC CCCCCCC CCCC CCCCCCC CCCCCCCC 
bBBBBBBbBBtfBBBBBBUKCCCC CCCCC lAAAAA^^BBBBBBBBHBHBBftR^CCCCCcAcC C CC CCC C C C CCC CCCCC CCCCCCCCCCCC C 
AAA A AAAAA^BBBHbBtBSBtfbBBBlIC CCCCCCC CC lAAAAAAAfc89nBB8B»BflB»*n«*CCCj^RCC CCC CCCCCCCCCC r t CCC CC CCCCCCCCCC CCCC 
|HBBtfDBBBBWH(*BB9KCCCCCCCCC kAAAAAA AAW9RHBBB8bH8BRH»BU9BHH^CCCC CCC C CCCCCC CCCCC r~ CCCCCCCCCC CCCC 
\»BBB8BBBBbBBb*»B< CCCCCCCCC I AAAAAAAAA»HHHBUBb8BBBB8BbBtiSBR 



AAAAAAAAAAAAAl 

aaaaaaaaaaaaaaS 



fccccccccccccccccccccccrccccccccccccc 



bBBBBB»tWBB8BB{ CCCCCC 
BBe8BBBbHB4B88E CCCCCC* 
AAAaAAAAAAAAAA«bD8BMBBHB*9dBt CCCCC 

jBBBBSBBBbBbBBBC :CCCC , 
AAAAAA.^,AAAAA^abBbBB9bbHBBBB« CCCCCCC _ 
AAAAAAAAAAAAApBBBBdBBHBtfBntB%E CCCCCCtCC 1 
AAAAAAAAAAAjflB8dBBBB<lBbdBtfB(CCCCC CCCCC 
AAAAAAAAAAj<^BBtft)»BBBB>»BBBBBi>t£CCCC CCCCC 
A AAA AAAjff V>»BBBBBBpBBBBBBt>BBtfcCCCC»wCCC r 
IBBBBH&BBBBHBBB^ J»BB Bi+bHBB^CC C C C C C CC C C C C 
«(MJ«tf9BB8BBBBBBBbBBHBBB«BiCCCCC CCCCCCCC 
BBBBBr»»BeHBBBBBBBb3BB*HBI*CCC CCCCCCCCCC 
BRBeBHBBb' .1MBBBB»BB»BHBj£CCCCCC CCCCCCCC CC]»AA/ 
ARBbBdWbBrfBHBBHBBBSBJJ^CC CCCCCCCCCC CCfCC CC( 
Mne&BBBb&HbBbBBHBBBJ^CC CCCCCCC CCCC CCCCC r« U*f * 
BB8B6b»Be-BrfBt:BrfJ**CCCCC CCCCCC CCCCCCCC CC |U>»Befj 
hBHBBBBBB*BHB»b1Tc:C CCC CCCCCC CCCCCCCCCC 4BBRUH8bSr*Bb 



8BBJBBHBBB8BBBBHBB*CCCC 
BBBBBBHe8HB«bBB»BBltcc': 
BRHBBBRBBBBHBBBHBtfTccCCl 
BbBBBbBBBBBBBBBBBfCCCCC 
HB8bBBB8BBHBBHbBHjC CCCCCC' 
BBBBBHBBBbBBBBRBBjCCCCCCCCCC 
BBMBBHBbRBBBBflBBHJCCCC CCCCCCCC 



.CCCCCCCC CCCCC cCCCCCC CCCCCCC 
iCC CCCCCCCCCC CCC CCCCCCCC 

;cccccccccccccccc 

'CCCCCC CCrCCCC 
ibBRflSBBBBBBBBBBJCCCCCc C' 
^BBBBBBBBB^BeBbJCCCCCCC 
jBUBBBeaflBBBBboBW: CCCCCCC 



BBfiHBUBBHBBBBBHHBttCCCC CCCCCCC CCCCC^B8BBBBHBBHB»8b- CCCCCCC 
&4RBtfB»BbBBBBRftl|j/cCCC CCCCCCCC CCCC CCCCT^BJBBBHflBBBKC CCCCCCC 
BBB^tJBSHBBBBfBBBtCCCUCCCC CCCCC CCCCCCC CCC CC^flflfiHBfi: CC CCCCCCC 
BBBBBBBHBBBBBBBBf CCCCCCC CCCCCT CCCCC ICCCCC^F AhhA cc/cCCCCCC 
9Bfc6ftBbb88Rft«E^ CCCCCCCC CCCC CCC CCCCCCCCCC CCCCCCCCCC CCCCCCC 1 
^HBflBBfCCCCCCCCC CCCC CCC CCC CCCCCC CCCCCCCC CCCCCCCC CCCC i 
BBBBBJC CCCCCCCC CCCCCC cCCCi-CCCCCC CCCCC CCCCCCCCC CCCC -C 
CCCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCC< CCCCCCCCCCCCCCCCCCCCCC 
'CCCC CCCC CCCC CCCC«XttC CCC CCCCCCCC CCCCCCCCCCCC 



KBBB BBBBBKjC cTTcCCCCwCl.CLCCl.ee CCCCCCCC CCCC 9BBBB88B99BBBBeHHj:CCC CCCCCCC CCCCC CCC CCCCCCC CCCCCCCC CCC CCC CCCCCCCCC 
^^CCccCCCl CCC CCCC ww u CCCC CCCCCC CC^B(iiflBBe8H98BBtirl8BrfcCcr CCCCCCC CCCCCCCC CCCCCC CCCC Cr CCCC '.CCC CCCCCCCC 
f CCCCCC CCCCCCCCCCCCCCCCCCCCCC CCCCCCC CCCBBBa By8888HB8t*>*BtCcCCCCgCCCCCCCCC CCCCCCCCCCCC CCCCCCJCCC CCCCCC 
C CCCC CCCCCCCC CCC CCCCCCC CCC CCCCC CCCC CJ^BBBBB&BtJbBBBHbdBtjCCCCCCC CCCCCCC CCCCCCCCCLCcCCCCCCCTCCCCCC CCCCCCCC CC 
cCCcCCCCCCcCCCCCCCCCCCvCCCC CCCCCCC^pBBBBBHBBflBB^BbBSHjfcC CCCCC CCCCCC C CCCCCCCCCCCCCCCCCCCCCCtCCCCCCCC 
C C CCC CCC r CCCCCC CCCC Cc CCCCCCcCCCCCCcJPb8BB88BBBB8B9«B!*H*CCCCCC'"C CCCCCC CtCCcCCCCCCrC CCCC >. r CCCCCCCC CCCCC 
fCCCCCCCC ' <"CCCCC.CcCCCCCCCCCt <*ct*BBBBbB«rtBBBB^^&-CCtCCCCCCCCCCCCCCCCCCCCCi CcCO C C C CCl. C C C C C CCCC CCCCC 

cccccCC(..c.cccccccccccccccuCCcCccc!*f«ttBBBBH8BBHt;dB9dt»dKj:ccc ccccc ;cc cc c cccccc cc c r r .c^cccc crcccccccr 

CCC' i" f ^' CCC CCCCCCCC CCCC?BBae*&*H«HfcrtMB^HHi*'CvCCC< fC" CCCCCCC CCCCCCCC .Ci d C'CCCf ' c etc CCC ccc 
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1 1 1 1 111 U 1 111 11 I Hi it VR*222222*22ii 222^2222222*22*222222222 222222*2***22^}*}** 1333333 3 333 33 
111 I Hill I 11 11 II 1 1 II i $22222222222 2222* 2222222* 2*i 2222 1 2222 22222*222 2222222(l_*********i***l* 
H I I ill 11 1 1 11 1 il illl 1 V22222222?222?22 ( 2222 222222222 222222222222222222222*222*X****** 33 333 3 
H 1. Hi 11 1 t. "Ui 1 1 IK 2222 2222222 222 *2 *2* 22222 22222 2 2222222222 2222 22222 222222^0**** 3 3 1 3 1 ] 
11 1 1 111 1 I li 111 1 11 1 1 11 1 1 >2222222222*22**2*22222*22**222222**22*22**2222t222222>22l%lYi))*M 
111111 11 111 1 1 1 11 HI 1 U222222?22*2*2*2 22222222**2** 2*222222*22* 22222222 2222222222X****** 
lllllll inillllllllll *22?2222222222?*2222?<222**222222222222222*22222222222222222222Tti3* 
I 1 111 1 1 1 1 I i 1 1 1 I III 1 1 1 '222222222* 2222*222 22 2*22**22222 2222 22222*2222222 12222222 112222221 _ 
111111111111111 till I JJ12* 2 222 2222 22222 222222c ' 222222*2 2222 222222 222*2222222222 22222222*227 
I lllilllllUHHlllf 222222*2222*22 *2222222< 2*22222222222222222222222222222222*22222222 

I 1 1 ill 111 11 111 1 11 Hf 222222222222222222222222*2t%:-??2>222222222222222222222222??22222222? 
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1 1 11 111 1 11 HI 1 iiJ222222222222222*222t>2222222222222222222222c22*222£22222222222222222*22 
I U Hill 1 I II 1 1 222 222 22 22222222 2222 2 222222 2 2 2 2222222222 222222222 2222 222 ?2i2 2222 2*2 2 ? 2 

lilllllll H 1 Hill 1X2222222222222*222222*2 222222222* 22*2 222222 22i 2 2*2 2 2*22222 222222 22222 222 
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Conclusion 

I ho examples presented here demonstrate the variety of Ways in which 
LANDUSL tan provide a b jsis lot the development ol computer-based laboratory 
exercises lor economic geography mmsus I he flexibility ol the program permits it 
to be utilized in couises nn neaily any level and the type ul exeiuses developed car 
be tailoied to any gioup ol students 

LANDl'Sl in still undergoing development (e g , it will soon be possible to stud) 
the ellects ul s\slenuliL t areal variations m productivity) and commenb and 
suggestions (mm useis aie welcomed * 
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APPENDIX A 



DETAILED OPERATING INSTRUCTIONS / 

Version 2.0 /Department of Geography 

/ Northwestern University 

LANDUSE 

Purpose: Development of land use and market area patterns according to 

the basic Thunen model. 

Description: Within a rectangular area tMs program develops the pattern 

of land use according to the basic Thunen model :< R = Y(P - C - TD) , 
where R is the economic rent, Y is the yield per unit area, P is 
the unit market price, C is the unit production cost, T is the 
unit transportation cost per unit distance and D is the distance 
to the market, The use producing the highest rent level is 
assumed to occupy the unit of land. In addition to N market 
centers and M crops, NROAD highways may be defined for a given , 
problem. Crop yields (Y) are assumed constant for each crop over 
the total area. * 

In general, the addition of a highway will oxpand the amount of 
usable land by reducing the effective distance from a cell to a 
market center and, thereby, reducing the transportation cost. 
The restriction of a, the ratio of on-highway to off-highway 
movement charges, to 0 < a < 1 implies that the effective distance 
along the highway is reduced by the factor <* below the sane 
distance off the highway. As it is possible for a problem to have 
multiple highways, it is necessary to consider all possible paths 
, including the direct path from a given cell to each market center 
and to select from these the minimum effective distance. Let 
i^, i^» ... i r be subscripts which identify r highways to be 

traveled in sequence from aicell to a market center, and assume 
further that this is the minimum effective distance path. Let 
D^ be the minimum d 1 stance from the cell to highway ij., the 

minimum distance from the market center to highway i p , and 

D ij> D i2» •**» D i tne distances traveled along the r highways. 

Then the effective distance D is given by the formula: 

D = D r ♦ a. D ♦ D ♦ ... ♦ a D. ♦ 0. 
C ij ij i 2 i 2 i r i T T 

The rectangular area relevant to program computations is between 
1 and 110 units horizontally and 1 and 90 unit*, vertically. This 
corresponds to the printed maps which are 54 Hnes at 110 
characters per line. There are 10 characters per inch horizon- 
tally but oni? 6 lines per inch vertically sy that the computer 
output must be scaled by a factor of 3/S to accommodate the printer 
and provide the same scale vertically and horizontally. All 
coordinates locating towns and roads are specified as arbitrary 
integer numbers, however, only the towns and portions of highways 
lying within the rectangular 'area defined by 1 <^ X < 110 and 
} < Y 1 90 will -appear on the output maps. 



ERIC 



53 

00057 



The program currently allows maxima of 10 market centers, 12 
crops, and 9 highways. Since the minimum distance must be 
computed from each cell to each town, problems could be devised 
using these maxima which might require a month or more of 
computer time. The lcrfer set of maxima of 4 market centers, 
4 crops, and 2 highways has been included in the program for 
student usage. Access to the expanded set of limits is pro- 
vided by punching a decoding "password" of up to 10 characters 
in Columns 71 - 80 of the problem control^ card . 

Data input to the program is an arbitrary sequence of problems 
each one beginning with a problem control card. The problem 
control card contains', among other things, N, the number of market 
centers, M, the number of crops, NROAD, the number of highways, 
and NEWCASE, the number of re-scale runs to be made. A problem 
is re-scaled by specifying a new origin and a new scale factor 
for the maps. The original (reference) problem is then recomputed 
with updated coordinates for the market centers and highways. 
As N must be greater than zero, a data card giving the market 
center coordinates always follows the problem control card. If 
•M > 0, the scale/crop-yield cards followed by the crop parameter 
cards appear next in sequence. If NROAP > 0, the highway data 
cards (one per highway) are next. Finally, if NEWCASE > 0, 
the re-scaling cards (one per re-s-ale for a total of NEWCASE) 
are included. It should be noted that program execution will 
be terminated and an error message will be printed out under the 
following conditions of the input data:- 

1) N < 0 or N > MAXN, where MAXN is cither 4 or 10; 

2) M > MAXM, where MAXM is either 4 or 12; 

3) NROAD > MAXR, where MAXR is either 2 or 9 ; 

4) a highway of zero length is specified, i.e. the initial 
and final coordinates are identical. 

Tie program output consists of A] an echo check of data inputs, 
Bj from one to four maps (depending upon the problem control 
card parameters), and C] sumary tables giving areas in cultiva- 
tion and levels of economic rent (if M > 0). The maps are:- 
1) a location of market centers nap; 2) a sphere of influence 
of roads map; 3) a land use map; and 4) a market areas nap. 
The first map always appears in the output, and the second 
appears if NROAD > 0. Maps 3) and 4) are printed if M > 0. 

The market center location map contains an X at the posit i6n 

of each town appearing within the previously defined rectangular 

area. Each highway or portion of a highway (if any) within this 
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rectangular region is also displayed with identifying numbers 
at the initial and final coordinates. An option exists for 
selecting the highway symbol. It may be either a dot (• ) > 
an asterisk (*) , or a plus (♦)". The default option is a dot. 

The Sphere of influence of a highway is the set of cells for 
which the average effective distance from the cell to all N 
market centers is less than the average airline distance to 
the towns or the average effective distance to the towns using 
all other highways. A cell within the sphere of ^influence 
of a highway is identified on the map with the highway number. 
Those cells whose minimum average distance is the average 
airline distance are left blank. If the average effective 
distances are equal for two highways with respect to a given 
cell, an equal (=) sign is printed in the cell. 

The land use map displays the regions of crop cultivations 
with a letter symbol printed in each cell associated with a 
crop. The market areas map indicates with a number symbol over 
the same region the market where the crops are sold. An equal (=) 
sign is printed in those cells where it is equally profitable 
to grow either of two crops or sell the product at either of 
two markets, respectively. Regions for which the rent is i 0 
are left blank. 

If M>0, distances on the map may be scaled by a multiplying factor, 
SCALE, This factor has meaning only if M>0 as a result of the 
effect of distance on the rent computations. All points on or 
off the maps are specified in terms of integer coordinates so 
that a unit increment in either the horizontal or vertical direction 
corresponds to a distance of 1/10" on the printed output. The 
factor SCALE is a non-negative floating point number (if SCALE 
* 0.0, it is reset to 1.0) which assigns a physical distance to 
dus unit increment. The physical distance may be either abstractly 
defined or in feet, miles, etc.; and 10.0*SCALE is the distance 
per inch of the printed output maps. This distance scale is 
printed at the bottom of each map if M^O. To obtain maps of dif- 
ferent regions within the sane problem requires the modification 
of the coordinates of all market centers and highways and/or 
SCALE leaving the remaining data the same. In this case physical 
distance must be invariant, i.e. modifying the value of SCALE or 
shifting the origin necessarily causes a consistent modification 
of the integer coordinates which will preserve actual distance. 
Thus, 



where 



3d = c 

c * the physical distance and is constant 
s * the scale factor 
d = distance on the map 
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Let us assume that data has been prepared for a problem 
which will produce output maps in a coordinate system (x,y) 
and that maps of the same problem in a coordinate system (y l ,y') 
are desired. Also, let s and s ' represent the distance scales 
in the uHprimed and prined systems; and let ( x q» v q) be the 
unprimed coordinates of the origin of the primed system. If 
( x l>yj) represents an arbitrary point in the unprined system, 
its corresponding point in the primed system is (xj,y*) where, 

xj = s/s' (Xj- x Q ) 

y[ s s/s' (yj- y Q ) 

All coordinates of market centers and highways must be transformed 

as above and rounded off to closest integer values. Note that 

the new coordinate system is completely specified by x Q , y Q , and s*. 

The following example illustrates the salient points of the above 
discussion. Assume that data has been prepared for a" problem 
having a market center located off the maps at ( x l >y 1 ) - (120,-20) 
with a scale of 4.0. For the same problem it is aesired to obtain 
maps of the region surrounding the m ket center and having a 
scale of 3.0. Let (x Q ,y 0 ) = (90,-40). Then 

xj = 4/3 (120 - 90) = 40 

yj « 4/3 (-2Q - (-40)) = 26 2/3 = 27 

and the new data would include the point (40,27) as the location 
of the market center. 

The program contains an option which will autoTaticall y re- scale 
any problem. The problem as specified by the data immediately 
following the problem control card will be called the reference 
case. Each re-scale data card which follows a reference case 
contains the coordinates of a new origin and a new scale factor.* 
The coordinates of the market centers and highways in the reference 
case are saved and new coordinates are computed relative to the 
reference problem for each re-scaling. Re-scaling may be done if 
M = 0, in wnich case the reference scale factor is assumed to be 
1.0. The number of the reference case is printed for each re- , 
scaling along with the x ar.d y coordinates of the new origin. The 
scale rano, i.e. the ratio of the reference scale factor to the 
new scale factor (s/s*) is also printed. Under *the existing input 
data structure a problem may be Te-s<aled as many as ninety-nine 
times. The user is cautioned not to specify a scale^ factor on a 
re-scale card so large that a highway is effectively reduced to 
zero length, i.e. the initial^ and final coordinates coincide. This 
will cause execution to be terminated. 
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D2ck Make-Up: A) Basic structure [CDC 6400] 

1) Systems cards 

2) EOR card i 

3) Program decks 

4) EOR card 

5) Data deck 

6) EOR card 

7) EOI card 

B) Data deck structure - The following sequence is repeated 
once for each problem: 

1) Problem control card 

2) Market center location card(s) 

3) Scaling and crop yield card(s) [if M > 0] 

4) Crop parameter card(s) [if M > 0] 

5) Highway parameter card(s) [if NROAD > 0] 

6) Re-scale card(s) [if NEWCASE > O] 

Coni Format: Probler, Control Card - FOP.VAT (512, 7A10) 

Cols. 1-2 N = number 6£ centers [1 IN f. 10 if password 
is used; 1 o^S <, 4 otherwise] 



Cols. 3-4 M« number of sxoj>vfM*i 12 if password is u^ed; 
M <_ 4 otherwise]" 

Cols. 5-6 NROAD = number of highways (NROAD 19 if password 
is used; NROAD ±2 otherwise] 

Cols. 7-8 ISYM < 0 )f highways are to be displayed on all 

map,; ISYM >^ 0 if highways are to be displayed on 
the l^cation~of market centers map^oUy. |lSifM| = 1 
if highway symbol is to be an asterisk (*) and 
| ISYM | = 2 if the symbol is to be u plus (♦) . The 
highway syiibol will be a dot(«) for any otner value 
of ISYM. 

Cols. 9-10 NEWCASE = number of re-scale cards for the current 
problem (if > 0) . If NEWCASE <^0, the problem is 
not to be re-scaled. 

Cols. 11 - 70 Any title information 

Cols. 71-80 Any word uith up to 10 symbols. If this word is 
the password, the limits on the number of centers, 
' crops and highways are increased. 



M arket Center Location Card(s) - FOP-MAT (8110) 

The X and Y coordinates are specified pairwise in integer format 

at field width increments of 10 in the following sequence. 
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Xj, Yj, the X and Y coordinates of center 1 
X-, Y~, the X and Y coordinates of center 2 



X N* t * ie X anc * Y coordinates of center N 
Note that 4 *] ♦ 1 cards are required.* 



Scaling and Crop Yield Cord(s) - FORMAT (8E10) + 
Cols. 1 - 10 SCALH = the distance scale factor. Ten times 
this number is the scale printed below each 
output map. If SCALE = 0.0, it will be reset to 1,0. 

Cols. 11 - (Y(I), I = 1, M) = the crop yields per unit 
area for crops 1 through M, respectively, at 
field width increments of 10. 



Note that [M/8] + 1 cards are required. 



Crop Parameter Card(s) - FORMAT (3E10)t 
Cols. 1-10 the unit market price 
Cols. 11 - 20 the unit production cost 

Co2i>, 21 - 30 the unit tiansportation cost per unit distance 

Note that there is one data card per crop per market center 
for a total of (N) (M) cards. The crop parameters for crops 
1 through M with respect to market center 1 are read first; 
parameters of crops 1 through M with respect to market center 2 - 
are read in next; etc. 



Highway Parameter Card(s ) - F0RMT (4110, E10) A 



CoK. 


1 - 


10 


X coordinate of starting point 


Cols. 


11 


- 20 


Y coordinate ot starting point 


Cols. 


21 


- 30 


X coordinate of ending point 


Cols. 


31 


- 40 


Y coordinate of ending point 


Cols. 


41 


- 50 


ALPHA , the ratio of on-highway to 








off-highway movement charges. 
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Note that there is one data card per highway for a total 
of NROAO cards. 



Re-scale Card(s) - FORMAT (2110, E10)t 

Cols. 1-10 X coordinate of npw origin in reference system 
Cols. 11-20 Y coordinate of new origin in reference system 
Cols. 21-30 SCALE = the new distance scale factor. Ten times 
this number is the scale printed below each 
output map if M > 0. If SCALE = 0.0, it will be 
reset to 1.0. 

* [X] means the integer part of X. * 

+ The E10 FORMAT allows the user to specify input data either in 
F or E form within a field width of 10 columns. The decimal 
point must be punched in any case. 

Note: The program LANDUSE requires the subroutines 

PR0JCT and REGION and needs 50000 8 central memory locations 
*.o compile. The program only requires 33000 g locations to 

load and go. It is written jn FORTRAN IV for the CDC 6^.00 
computer. Execution tine is on the order of from 5 to 20 
seconds for problems having up to 2 highways, but increases 
exponentially with an increase in the nunber of highways 
beyond 2. A strong factor affecting the amount of running 
time is the complexity, i.e. multiplicity of intersections, 
of the highway network. 



Programmer j : Bruce Anderson and Duane F. Marble, August 1970. 
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APPENDIX B 



THE GEOGRAPHY PROGRAM EXCHANGE 



The Geography Program Exchange (G^E) has been established to assist 
universities and other non-profit organizations with the interchange of computer 
software which relates to problems of a geographic nature. The GPE operates on a 
self-sustaining basis with the assistance of a modest grant from the Mathematical 
Social Sconces Board Its services arc made available to qualified users anywhere in 
the world 

The GPE maintains a central file of specialised computer programs with 
associated documentation and test data sets, and makes copies of all or part of this 
material available to qualified users at cost. The programs held by the GPE 
represent specialized materials developed for use in geographic research and 
teaching, no attemp* is made to duplicate standard statistical programs. Distribu- 
tion of available programs is in the torm of punched cards or BCD tapes, and 
aooouatcd documentation is also distributed in machine readable form. The GPE 
aJso produces a brief guide to its program holdings which is available without 
charge. Current program holdings are based mainly upon the specialized program 
libraries of the University of Iowa, Michigan. Michigan State and Northwestern 
University as well as individual programs from many other institutions. The GPE 
actively solicits programs and. if the program submitted appears to be of enough 
general interest, it will be added to those distributed. All software handled by the 
GPE is obtained and distributed on a non-royalty basis. 

In addition, the GPE also publishes an irregular newsletter which, at present, is 
mcorpon ted in Current Research Notes in Quantitative and Theoretical Geography. 
This is distributed without charge by the International Geographical Union 
Commission on Quantitative Methods and requests for copies should be directed to 

Professor Forrest R. Pitts 1 
Social Science Research Institute 
University of Hawaii 
Honolulu. Hawaii 96822 U.S.A. 

Although the GPE operates with limited resources, it does attempt to provide 
assistance to those departments and individuals who are beginning to make use of 
the digital computer as a research and teaching tool. 

Overall management of Hie GPE is provided by a board consisting of Professors 
Duane F Marble (Northwestern). Waldo R. Tobler (Michigan). Robert I. W;Uick 
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(Michigan State), and Forrest R. Pitts (Hawaii). Professor Witfick is in charge of 
operations and all inquiries relating to activities of the GPE should be directed to 
him at the following address 

Professor Robert Wit tick 

Computer Institute for Social Science Research 

Michigan State University 

Bast Lansing. Michigan 48823 U.S.A 
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APPENDIX C 

KEYPUNCH INSTRUCTIONS 




» Vogelback Ce ,'JUng OnUr 

Nor chwo> toi :i I'tuvor «ity 
Evanston, I i 1 moii 

IBM 2^ Printing Card Puruj] 
Ccncra 1 Ope ra t ion . 

Inform t ton is punched in a card by opera ting a typewriter t> lc ko\boaid, 
which punches letters, digits, or special character into ^ard column.. The 
operator pnssi ? the "feud" button on the keyboard which rvves a card dour, 
from the hopper. When the "register" button is depr es;>ed , N tht* first co!%rn of 
the card is positioned unde: the punch mj station. The card.s are puncru d Lion 
left to right, one column at a t l~e . After each card is punched, the car*' i« 
audo^itical ly position d for the ne^t colur. w. Attjr all SO colu-rv* »il punched ( 
the card autcrvu ical ly is registered urder the reauir.g station, thrvu t h •. hi^h 
it passed column by colur n in phise with the card following it at the 
punching station. Wi._n the cad at the punch in.; s'Tticn Us re leased, the card 
at th^ reidmg station neves Into the card stacker. 

Because p.is &.»£.■■ through the rtidirg station is in phase with p.is^n.e 
through the punching static", holes s.-rwd at the read stition can be Juplieauv' 
into the following card pacing v.nder the punch station. The abiiitv to dupli- 
cate infoircdtitoi fron card to card reduces the an our. t of r*inua 1 punching r.tded. 
Data corrwn to a group of cards ruy be punched into the first card ot a grou,« 
and duplicated into the rer-iimng ones. Depressing t\ e "dunlin ate" U\ jljsi 
accor-plish this, or it ca be done au torv 1 1 ca 1 ly by the keypunch opL'.atuu' order 
preplanned control (drun or program card). 

Other functions that cc" be put urd^-r au tor-.it lc control axe printing, 
skipping colurns, and ejoctir , .' J and feeding cards. According to the re t u r ei ont- 
of the card for--i., a pattern can be established so that th^ keypunch i an *A:p 
columns, duplicate field-. , or release the card when punching is fin. -hud but 
before colu'nn 80 is leached. 

Instructing the keypunch to skip, duplicate, or release upon rc.iehim. 
certain colurn locations js known as programing. The operator progij-*^ the 
keypunch by punching a control card, wrapping tho card arornd a ^vlinde: e: 
drun, and inserting the dru^ on a spindle in the nachir>«. The holes in tJ e 
card control the operation of the machine, causing it to perform specified 
functions autor-jt ' ca 1 ly. 

Opera t ing Features 
Cftrd Hopper 

The card hopper is located on the upper right side of the machine. The cjrd> 
are placed in the hopper face forward, 9's ifoun, and are fed front card hr^t. 
A card Is fed fron the hopper to the <.ard bed au tnrr i 1 ua 1 I y or by deploys ion 
of a c,»rd feed key. The iir-st two cards to be punched ru^t be fed bv ^ > ih - 
press ion, but all other cards in the hopp' x ruy b< ftd au tuna t ua I 1 y . umU i 
switch control. 
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Punching JtHt_ioti 



Punching is perfotru»d at the fu.st of two stations in the cud bed through which 
Che cards pass. To start an operation, luo c ards arc fed into the cud bed at 
the right of the punching station. A-> the second card is fed in, the first Laid 
is automatically registered for punching that is, it is positioned at the 
punching station. 

p cad .nt, St it inn 

The reading station, where thi c.v.ds are read for clupl u atmg, is located approp- 
riately the distance of one card to the Left of the punching station. Each eatd 
that has b«"n punched p.i,,^ through tru read station as the next catd is being 
punched. The two cards- rove in synchrunisn, column by colum, and information 
to bo duplicated Lj transferred from the first Card to the second. Kead uig fiom 
one card to aaolhei car. be concr^lled, field by field, so that only the desired 
mfot ration is duplicated. 

This feature eliminate, cud handling fei the duplication of inf,. rr.it u«n from 
cards pn pared during th< punch ir. tl operation. Card handlir is r, ( v. • ar> only 
wh»*n an operation r<c.uues duplication ft or prepunch*'*' ma.ter card-.. 



Card Staetf r 

The card stacker is located on thv aprer iei't side of the machine on a level 
with the hopper. After each card passes the read stall"*, it is fed into the 
sticter autit-vit lea Uy or by key depression. When the ds are rencvod from 
thf stacker, they arc in thc:r original sequence. 



This key is located below the card bed between 
long as it is held dep: e..sed,> the cards at the 
spaced con tipao'is ly until column ! is reached, 
card, vp icr. comroli .k»ppn: 6 at 4 -! duplicating, 

Program I'nit 

Automatic skipping, au terra tic duplicating, and shifting fron numerical to alpha- 
betic punching and vice ver^a ate controlled by the program unit. f\u h oL Chc-e 
operation* is designated by a specific code r corded in a program card. (See 
Program Card Codes) The progran eird is then fastened around a program drum 
and inserted in the machine, uh*>re it Is read by a sensing mechanism. The drum 
revolves tn step with the movc-'uit of the caids past the punch and reui stations 
so that the program cod»*s control the operations, colurn by column. 



the reafl and punch stations. As 
punch and n ad station^ ai e bitk- 
At the same t ire , the progran 
is also back -.paced. 
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Program Control Lever 



Operation of the program unit is control led by the progiam cutitrol lever, 
located below the program unit. When this lever is turned on, the program 
sensing mechanism is lowered so that it rests on the program drum, and the 
codes punched m the program card control the various automatic opera tionb 



Column Indi c ator 

The indicator, located at the base of the program drum holder, indicates to 
the operator the next column to be punched. Refer to this uidi^atoi «s a guide 
for spacing or backspacing to a particular column. 



jL^A u _llds_ ._~__Nurher 3 - Combination alpha b etic 'indjnuj^cri cal keyboard , plus I I 
sj>ec_wi l_ char ac ter s . " ' 

The cunbinition keyboards have the be^t features of both a typeuiiter and a 
numerical key punch. The letter keys are arranyd for operation by tru- stan- 
dard typewriter touch system, while the digit keys are placed so that a rapid 
three-finger tu'^i :i syste- crni be used. The usual nurerical keys on a t>pt~iitei 
have bten e I i - r- 1 t»d, instead, a group of dual-puipose keys at the rivht seiw> 
for digit »■» wrjil as letter punching. This permits :.u-eri«.a! punchm., uith :lu 
ri^it hand. The touch bwtt'n for the ten nur.eiical key * is. indes f ir^cei Uu 
digits I, U, a-d 7, middle finger for digits 2, 5, and S. and ring finyVi loi 
digits 0, 3, 6, ?rd r The punching of a digit oi a Utter with any ,ot the 
combination keys de f _uds upon the shift of the keyboard. For ovanple, depus- 
sion of the 4- J key punches a uhen the keyboard is m numerical sh»lt, but 
a J when m alphabetic shift. This shifting is similar to uppei or !ewei ubo 
shifting un a standard typewriter and may be controlled autorat ica 1 U H tlu 
program unit or manually by key depression. 

Tru fourth row of keys contains four special character keys at the Kit, Iho-o 
keys punch ei^ht characters as shown on the key tops, four in uuierual and 
four in alphabetic shift. 



Program Card 

The program card controls the automatic operations for th*_ cox i e-pond i n v, toluru 
of the cards being punched. The control punching : -,jir.! ir trie prjffram card 
depends on th<« functions to be controlled, e.g., skipping, ii*r>licatin#?, and 
alphabetic punching. 



Field De f init' on, ( it punch 

A 12 punch (+) should be punched in every column except the first (lift-hand 
position) of every field to be skipped, duplicated, or rami* 1 1) puiivhed. These 
12's. serve to cuntmue to the end of a field any skip or duplication started 
within that field. 
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An 11 punch in the first coiu-'i of any field au "uitically starts a skip, 
which is continued over chat field by the 12's punched in the iei"-iiniri£ c 
of the field. If a sir^le tolm-i 1=. to b«» auto, a i u a 1 \ > sripp»d, it is p 
w^than li. 1hi3 coding op°ratet> in conjunction with thy auttnitu skip- 
duplicate switch, which must be on to start the skipping au tora t ica Uy . 



Autor at'.c Pup I icu t 'n (0) 

A zero punchcJ in tho first eolurn of any field auto^-at leal ly starts duplica- 
tion, which is continued over that field by the 12' s punched in the rorjinuiy, 
coiur-ns of the f i« I 1 'fa single culur, n is to bf au tona t ically duplicated,* 
it is punched with 



When the pio^i.n iiid la in the machine, the comb l nat io* keyboard is norna!l\ 
in nupcrvcal shi't, and depie ?siun of an; one of the two-pu; po-»t kevs caus»t > 
a f i£iir»- to be punched. To punch a letter, the combination ke)hoird trust ho 
shifted for alphabetic punch irn;. This sh.ftim 1 , is p r furred aut on at ti a I 1\ h\ 
a I in thr program card in each colunn of the alphabetic field. 



Prq/ra-i Card Cod'-s 



Code 


Funct ion 


12 + 


Field Def mi t ion 


li • 


Start Autorvitie Skip 


0 


Start Autorvitie Duplicatiot 


1 


Alphabetic Shift 


2 


left ?ero Pr int 


3 


Pr int Suppress ion 



Pro gram _ Q ru _^ 

J h «• program card is mounted on a program drum lor insertion iu llu luuhuu . 
Die program drum has a cla-ipm^ strip to hold tlu tatd, and a handU yu the 
to tighten or release the strip. To fasten a card around the dmri, laki lit 
fol lowing s tcps . 



top 



1. Hold the drum in a horizontal position with the handle to the ri^hf lutn 
the handle away (counter -clockwise) as far as it will go. This loosen- the 
smooth edge of the clamping str r . 

2. Insert the column 80 edy,e of Ch* card under the smooth edge of the clamping 
strip Wo alignment check holes in the clapping strip make it possible to 
see that the card is flush with the metal ed&e under the sttip. Hie card 
should be positioned so that the 9's edye is .if a ins t the run of the drum 
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3. Turn the handle to the u'nUr posit iun. This tightens the smooth edge of 
the clamping strip a-.J loosens the toothed edge. 

_j * * 

4. Wrap the card tigntl) around the dru'n and inset t the cotut..n I edge under the 
toothed edge of the clanging strip. 

5. Turn the handle toward you (clockwise) as far as it will go. This fastens 
the toothed edge of the clamping strip. The drum is. now ready ry be inserted 
in the machine. 

To rtmove a card from tha drum, the above procedure should De followed in the 
reverse order . 

With the p-ogran sensing mechanism raised, the drum is placed on the spindle 
undei the center cover of the punch, positioned so that trie aligning, pin falls 
in the aligning hole in the column indicator dial. The program control leve 
is turned on to lower the reading star wheels onto the program card, and the 
release key is depressed.^ to engage the reading mechanism fully. Whenevei the 
drum is to be removed, the program control leaver irust first be turned off to 
raise the star wheels. 

When the pu nc h is operated without pro^r^n control, the keyboard _is in_a_lpha - 
betic sh ift. When the program control lever is .turned on, the keyboard is in 
numerical shift for all columns except those punched with the alphabetic field 
coding in the program card. If it is necessary to punch digits with the juov;iam 
off or in an otherwise alphabetic field/ the x .u m er i c a 1 s h i f t__kev is lie 1 J doun 
while the digit keys are depressed. 'Similarly, when it is necessaty to putu h 
some letters in an otherwise numerical field, the alphabet ic shjJJ ko\ is held 
down while the letter keys are depressed. 

\ 

SV. ip K ey 

The purpose of this key is to skip any field that is coded with 1 2 * s in the** 
program card.. Depression of the key starts the skip and the 12's, in succeeding 
columns of the field, continue the skip for the rest of the field. 

Dash Skip (or) Das h Key 

When the keyboard is-un numerical shift, depression of thi y punched an 11 
hole and causes skipping for* the rest of the field as dote* .d bv the miuo- 
»lvc columns punched 12 in the*pr*gram card. When the keyhoaid js m alj''>iKt-u 
shift, depression of thus key punches fin 11 hole but doo*, not cause sktpiMiu,. 

Duplicate Key * ^ 

Any -tnToriracion, manually, punched in one card and coT,on to one or mote succes- 
sive cards, can be duplicated by use of th lfi key. Depression of this key at 
the beginning of a fiel<3 starts duplication, and the 1 2 * s punched in the program 
card continue the duplication for the r*st of the field. 
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t 



Key bo .ad S 




Shift 
key 



C 



Space Bir 



(Aux j f A I pi-1 

Shift 
key 



Pu ri c h ii '.', i\ 0 s 

Keyb 1-18 car, oe iVprvs Jt 'J orly uhen the keyboard is m alphabetic shift Co 
punch the lettur* mdi^at^J, If ore cf tlwse r j> b i' depr^bed vnllo the r e> - 
board is ir. ru ?r«*cal shift, the rachi«ie locks. Operation can be resu- .1 by 
releasing th j card, or b\ uepressi:.; t alphabetic shiit ke> v*hich will ciu «e 
the letter to be punched, 

Conbmation ke>s 19-29 can be dtpre*»sej '-hen the keyboard is in either i ur.-rual 
or alphabetic sh »ft to purch the characters controlled by that key. 

Kevs ~C-^3 par^h multiple digits for eight sp^ial characters. The> can de- 
pressed when tho keyboard is m either r.iirvr icaL or alphabetic shift. 



Space Bar 

The space bar car be depressed a: any^ tire in a rranual field to clause bpa*. in£ 
over era colurv. of tha card. 



Functional kevs 



30 NLM (numerical shift) shifts the combination ke> board into nu .ct ual posi- 
tion as long as it is h*ld down. 

31 ALPH (alphibetic shift) shifts the rorcbina t ion ke/ board into alphabetic 
position as long as it is held down. It is noma 11} u*ed to permit the pushing 
of letters in an othoi-'ise rurorical field. 

32 DUP (duplicate) causes punching by reidirg fro*n tWe preceding card. 

3* REI. (rele:t*e) causes the cards at the purch and read station* t ti* be idwiued 
conpletely past those nation*. Fields pro^ra-u ed ir: auter.iCic dup 1 ic i r i.'n 
beyond the point of relea>o aru punched tti the card before rel«-Js>e is u 1 ' pjaiji. 
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35 



hKKD (card feed) jraases *i card feed cycle 



1. Feeds a card fror.j the hopper. 

2. RegLstets the card', at the punch and read stations. 

3. Stacks th.j car< f fr^n the left* of the card bed. 

36 SKI use^ skipping of the field for which it is depleted* 

i7 Riv, i^ard r-.^ster) l* used primarily when in-*ertir.£ caids nunuall\. 

44 MULT PCH (nultipl-' punch) is held down to pt_ i m 1 1 the punching ot note than 
one digit in a colurn. The ke> board is m nuMuncal si lft when this key is 
depiessed. ' 

(Note* An exanple of the use of the inuitiplo punch Ley is* 

| END OF RECOKD ' 

Jt A i A A i A A v A * A * A it « A * * * * A * A * * V A .' A * A A w * * t. * V 

Satrpl** [iM.'.rj i v-aid f"t Fortra n stnt'j - jnts . (Note that "A M punch defines an 
alphabetic field length.) 



f 



73 SO 



col 1-5 nuner c (j> tat er.cn i nunber) 
6 numeric (continuation) 
7-72 alphabetic (to punch nureric information, depress NUM key) 

73-80 aucornatic skip and release of card 

A * * * A * #• A A A A " *„' A * A X * A x A A A A A A V A * # x * * > A ^ * 

Sample progra-i card for d ata c rds . 

1 6 11 2j 31 6j , SO 

'col. 1-5 ID (nui.er ic) ' 

6-10 Skip 

U-20 Numoru 

21-30 Skip 

31-60 Duplicate preceding card <nuneric) 

61-80 S*ip and release card 



Informal »or fi>r tin* instruction iheot was ohUi'mJ f:o:n the IP" pi hi it Us on 
IBM Kcfe reia i n..j \ 2 : '« Ijrd Pun. h_ ?»j Pri n_' - fj r^i Pum.li , Ft'i"i A 2 - 0 "< : , 
193°, and .\».nr t.i Ja: 7a , Da t .f JProio > c n,.. >t Northwe , t ern I'ulwisUs I'u , f°f> 5 , 
pp. 'iS-Sl . ' 
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